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Preface 


This book describes in 2 parts experimental data with simple explanations (Part I) 
and itinerant electron theories (Part II) about magnetism and its related properties 
of 3d-intermetallic compounds. Unlike 3d-metal alloys and oxides, theoretical as 
well as experimental studies on 3d-intermetallic compounds such as 3d-pnictides 
and chalcogenides, on which we focus in this book, seem unfortunately delayed. 
The objective of this book is to motivate active studies in this field in the future. 

We discuss in detail magnetic and related properties of the 3d-transition-metal 
pnictides and chalcogenides, which include the intermetallic compounds expressed 
as MX and M>X, and their mixed compounds M,_,M’,X, MX, _yXy and M>_,M{X, 
where M (M’) is a 3d element and X (X’) a pnicogen (P, As, Sb, and Bi) or a 
chalcogen (S, Se, and Te). Most of the MX-type compounds crystallize either in 
the hexagonal NiAs-type structure or in the orthorhombic MnP-type structure which 
is regarded as a distorted NiAs-type structure. Crystallographic phase transition bet- 
ween the NiAs- and the MnP-types occurs in some of MX-type compounds when the 
temperature changes. The M,X-type compounds crystallize in the tetragonal Cuz Sb- 
type structure. As described in detail in this book, many of the compounds mentioned 
above exhibit very interesting magnetic and crystallographic phase transitions caused 
by various means such as change of temperature, applications of external magnetic 
field or pressure, and change of the composition x in the case of mixed compounds. 
This variety of phase transitions seems to have motivated many researchers to in- 
vestigate the group of compounds mentioned above. We, one experimentalist and 
three theorists, had been collaborating for a long time to clarify the mechanism of 
phase transitions as well as the peculiar magnetic properties appearing in the group of 
compounds mentioned above. Considerable progress in understanding the properties 
of the compounds mentioned above has been achieved; however, there still remain 
various phenomena not clarified yet. Besides the fundamental problems, those of po- 
tential applications are also considered for the phenomena observed in some of the 
compounds mentioned above. 

In Part I the crystallographic and magnetic properties are described together with 
simple phenomenological analyses the consideration of applications is given in a few 
sections. In Chap. 1 we show the method to prepare a polycrystalline sample, crystal 
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structures and the phase diagram of a typical MX-type compound. In Chap. 2 we sur- 
vey basic magnetic data such as magnetic transition temperatures, magnetic moments 
and magnetic orders, etc. for the compounds of MX- and M2X-types. Chapter 3, the 
main part in Part I, presents detailed descriptions of magnetic and related properties 
for respective compounds. In Sects. 3.1—3.6 properties of the compounds containing 
Mn or Cr, which are the most interesting from both the fundamental and applica- 
tion view points, are described in detail, together with phenomenological theory to 
explain magnetic properties, and application of MnAs as a magnetic refrigerant is 
also included in these sections. Properties of the M2 X-type compounds are shown in 
Sect. 3.7. Many pages of this section are used for the Cr-modified Mn2Sb. 

Part I, where theoretical results are shown, is divided further in to two chap- 
ters. In Chap. 4 theoretical results of the NiAs-(MnP)-type compounds are shown. 
The magnetism of materials arises from 3d or 4f electrons. There are two differ- 
ent models to describe the magnetic electrons: the localized model and the itinerant 
electron model. Theories shown in Part II are based on the itinerant electron pic- 
ture. In Sects.4.1 and 4.2 we show and discuss in detail electronic band structures 
calculated by the first principle method for paramagnetic, ferromagnetic or antiferro- 
magnetic states of the MX-type compounds in the 3d-pnictides and chalcogenides. 
In Sects. 4.3-4.7 we theoretically investigate various experimental data on the ba- 
sis of the calculated electronic band structures. In Sect. 4.3 we clarify origin of the 
anomalous behavior of paramagnetic susceptibility and thermal expansion of MnAs 
by taking the spin fluctuation into account. Correspondence between Fermi surfaces 
and de Haas-van Alphen oscillations in NiAs and NiSb is discussed in Sect. 4.4. 
Pressure induced magnetic transition in CrTe is discussed in Sect. 4.5. In Sect. 4.6 
we present a theory to search for the most stable magnetic ordering in FeAs, where 
instability of the paramagnetic state is investigated from the wavevector-dependent 
susceptibility tensor. In the final Sect. 4.7 structural transformation from the NiAs- 
type to the MnP-type structures is discussed. On the basis of the calculated band 
structures, we investigate the instability of the NiAs-type structure against MnP-type 
lattice distortion by taking the electron-lattice interaction into account. 

In Chap. 5 we show the theoretical results of the Cu2Sb-type compounds. Crystal 
structure and magnetic ordering are shown in Sect. 5.1. In Sect. 5.2 band structures 
of the Cu2Sb-type compounds calculated by the first principle method are shown, 
followed by a comparison between experimental and theoretical results in magnetic 
ordering and optical spectra. 

This book is the result of a longlasting close collaboration between experimen- 
talists and theorists, which is also the tradition of the study of magnetism in Japan. 
It would be our great pleasure if this book proves helpful to graduate students, 
researchers and engineers in the field of science and technology of magnetism. 


Miyagi, Osaka Kazuko Motizuki 
October 2009 Hideaki Ido 
Tadaei Itoh 


Masato Morifuji 
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Part I 


Experimental 


1 


Basic Properties of 3d-Pnictides 


1.1 Main Compounds of 3d-Pnictide 


Among various intermetallic compounds between 3d transition metal M and 
pnicogen X (X = P, As, Sb, or Bi), compounds expressed as MX and M2X are 
of special interest due to their rich variety of magnetic properties. The MX com- 
pounds crystallize basically in the hexagonal NiAs-type (B81) structure, however, 
some of them show distortion into the orthorhombic MnP-type (B31) structure. 
The M2X compounds take the tetragonal Cu2Sb-type structure. The MX and M2X 
compounds are listed in Table 1.1. 

In this table, the MX compounds denoted by boldface letters take the NiAs- or 
MnP-type structure. The M)X compounds written in boldface take the CuzSb-type 
structure. It should be pointed out that the ternary compounds MnAIGe and MnGaGe 
also have the CuzSb-type structure. Some compounds take nonstoichiometric com- 
position such as M;,X (x > 0) although they are written as MX in Table 1.1 [1]. For 
example, MnSb actually forms a compound Mn;+,Sb with 0 < x < 0.2. Some other 
NiAs-type compounds also have similar tendency [1]. As the excessive Mn atoms, 
which occupy interstitial positions, have a strong effect on magnetic properties, a 
special attention should be paid for preparation of stoichiometric MnSb compound 
[2]. For the NiAs-type structure, we can regard that M atoms with smaller size than 
that of X atoms occupy interstitial positions of the hexagonal close-packed struc- 
ture formed first by the X atoms. Therefore, this may permit existence of excessive 
M atoms. Among the compounds shown in Table 1.1, only Cr-, Mn-, and Fe- com- 
pounds take magnetic ordering. 


1.2 Crystal Growth by Sintering Method 


Crystal growth is the starting point of experimental studies of solid state physics. 
Although there has been rapid progress in equipment and technique, we show here 
a conventional method to prepare a polycrystalline sample. This method is called 
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Table 1.1. Compounds between 3d metal and pnicogen 


P As Sb Bi 
Ti TiP TiAs TiSb - 
Vv VP VAs VSb - 
Cr CrP CrAs, Cr2As CrSb - 


Mn MnP,MnoP MnAs,Mn2,As MnSb,Mn2Sb MnBi 
Fe  FeP, FeoP FeAs, Fe, As FeSb - 


Co CoP, Co2P CoAs CoSb - 
Ni NiP NiAs NiSb NiBi 
Cu - - Cu,Sb - 


The compounds written in boldface have the NiAs-, MnP-, or 
Cuz Sb-type structure 


a sintering method, in which mixture of powders of constituent elements settled in 
an evacuated silica tube is sintered in a furnace. This method is applicable for most 
intermetallic compounds. 

Usually particle size of metal powder is about 100 um, but larger or smaller par- 
ticle sizes can be acceptable. Mn can be ground easily in a mortar of agate; however, 
other 3d metals are very difficult to pulverize. It is recommended to obtain the pow- 
dered metal from specialty store. Pnicogen is ground easily, but special care is nec- 
essary for phosphorus and arsenic. As for phosphorus, we use red phosphorus (not 
white phosphorus), which may explode while grinding. One can purchase powder of 
red phosphorus. Since As is toxic (As203 is deadly poisonous), one must wear mask 
not to inhale the fine particles and be careful so that particles do not scatter. Tidy- 
ing up after experiments is also necessary. Purity of all ingredient should be more 
than 99.9%. 

We show procedure of sintering for the case of MnAs. First, we wash electrolytic 
manganese in dilute nitric acid and then in water. Next we grind it in a mortar of 
agate. High purity bulk of As taken out of a vacuum glass container is crushed into 
powder. Powders of Mn and As are weighted in desired proportion and mixed in 
a mortar carefully so that the powder does not scatter. These processes should be 
finished as quickly as possible to avoid oxidation of the powdered metals. Then, the 
mixture of Mn and As (total weight of 5—8 g) is sealed in an evacuated quartz tube. 
An example of shape of the quartz tube is shown in Fig. 1.1. 

The neck of the quartz tube should be carefully prepared. If this part is too thin, 
keeping the tube in high vacuum may fail, and ingredient powder cannot go inside 
through the neck with too small internal diameter. 

The mixture sealed in the evacuated quartz tube is heated in an electric furnace. 
Since As sublimates at 613°C and has large vapor pressure at high temperature, fur- 
nace temperature must be raised slowly spending a few days up to 500°C in order to 
avoid explosion. Then the temperature is kept constant for another few days so that 
rough reaction of the mixture finishes. In the final stage, the temperature is raised 
up to about 800°C (you can heat up rapidly this time) and kept for a week. We can 
take out the reacted product after cooling down. As the sintered crystal thus prepared 
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100mm 


20mm 70mm 
+ | +> 


evacuation 
—> 


i, 


ingredient [ outer 6 4~5mm 
inner 6 2~3mm 


Fig. 1.1. A sample tube of quartz. After evacuating up to about 10~> Torr, the neck of the 
quartz tube is fused and cut by a gas burner. The dimensions shown in this figure are typi- 
cal ones 


may be inhomogeneous, the sintering at 800°C for a week is desirable to be repeated. 
Namely, the reaction product is pulverized and sealed again in an evacuated quartz 
tube. One can heat the sample rapidly this time. Longer the period of sintering, better 
the quality of the crystal. 

As phosphorus also has risk of explosion, the procedure mentioned earlier must 
be taken for compounds containing P. Other compounds also can be prepared in the 
similar way. The diagram of binary alloys [3] will be helpful for determining the 
sintering temperature. 


1.3 Crystal Structure of NiAs-, MnP-, and Cu,Sb-Type Structure 


1.3.1 NiAs-Type Structure 


As the MnP-type structure is regarded as distorted NiAs-type structure with small 
shifts of atomic positions, these two structures are similar to each other. We explain 
below these crystal structures connecting them with each other. 

Figure 1.2 shows the hexagonal NiAs-type structure. Pnicogen atoms form a lat- 
tice like the hexagonal close-packed structure, and 3d metal atoms a simple hexag- 
onal lattice. Table 1.2 shows lattice constants of the MnX compounds at room 
temperature. 

As seen in Table 1.2, with change of pnicogen from Bi to P, the lattice constants 
decrease. In particular, decrease of the a-axis is remarkable. It is noted that a-axis of 
MoP is 73% of that of MnBi, corresponding to the fact that radius of P atom is 75% 
of that of Bi. Hence we see that the a-axis is almost proportional to radius of the 
pnicogen atom. Such tendency is also found in the compounds CrP, CrAs, and CrSb. 


1.3.2 MnP-Type Structure 


The MnP-type structure is obtained by introducing an orthorhombic distortion as 
well as slight shifts of atom positions for the hexagonal NiAs-type structure. As 
distortion and shifts are small, these two structures are basically similar. However, 
even the slight difference affects magnetic properties through change of electronic 
structure. 
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r 3d metal (M) 


(0, 0, 0), (0, 0, 1/2) 


C) pnicogon (X) 
(2/8, 1/8, 1/4), (1/8, 2/8, 3/4) 


Fig. 1.2. The NiAs-type structure (B81) 


Table 1.2. Lattice constants a and c (room temperature) of MnX [4] and radius of the X atom 


(MnP) MnAs MnSb-— MnBi 


a (A) 3.17 3.68 4.15 4.34 
c(A) 5.26 5.72 5.78 5.98 
Atom radius of X (A) 1.28 1.39 1.59 1.70 


Note that MnP takes the orthorhombic MnP-type structure 
(a and c are lattice constants corresponding to those of the 
hexagonal NiAs-type structure) 


tub 


ba f 
@@) M: (0,0, 0); (0, 0, 1/2) 
OO, X: (2/3, 1/3, 1/4), (1/8, 2/8, 3/4) 


Fig. 1.3. Relation between the NiAs-type and the MnP-type structure projected onto the 
c-plane. (a) The NiAs-type structure (The thick lines denote the unit cell. See also Fig. 1.2.) 
(b) Distortion into the orthorhombic lattice and displacement of M atoms along the b-axis are 
shown. The displacement is written as ub using a displacement parameter u. (c) Displacement 
of X atoms in the b-direction. The displacement is written as wb (w: a displacement param- 
eter) whose origin is 1/3 or 1/6 of the b-axis. Note that the displacements and distortion are 
enlarged in order for the illustration to be legible. See also Fig. 4.38 


We explain the difference between the MnP-type and the NiAs-type structure. 
Figure 1.3 shows atomic positions projected onto the c-plane for the NiAs-type and 
MnP-type structures, where displacements of atoms in the c-plane are also illustrated. 
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In Fig. 1.3, crystallographic parameters u and w indicate atomic displacement 
along the b-axis. The parameter w is one order smaller than u. Distortion from the 
NiAs-type to MnP-type structure is expressed with a parameter 6, which denotes 
deviation of the b-axis from J 3a-axis : 


5 = (b— V3a)/vV3a. 


Atoms also shift along the c-axis as shown in Fig. 1.4. 

In summary, the MnP-type structure is characterized by the five crystallographic 
parameters: the displacement parameters of metal (u and x) and anion (v and w) and 
the distortion parameter (6). The NiAs-type structure is thus expressed by setting 
u=v=x=w= 0 =0. The crystallographic parameters are listed in Table 1.3 for 
MoP and CrAs, which take the MnP-type structure. 

Other MnP-type compounds have similar values of the parameters. Since w 
and x are generally smaller than wu and v, only u, v, and 6 can be determined in 


Fig. 1.4. Unit cell of the MnP-type structure. Displacements of M atoms and X atoms in the 
c(A)-direction are given by xc and vc, respectively. The parameter x is smaller than v by one 
order. In the expression of lattice parameters of the MnP-type structure, (A,B,C) is often used 
in addition to (a,b,c). Note that displacements are enlarged, and displacement of X atom in 
the b-direction (w) is not shown to avoid confusion (see Fig. 1.3) 


Table 1.3. Crystallographic parameters of MnP and CrAs in the MnP-type structure (room 
temperature) [4] 


6 u v x w 


MnP 0.077 0.05 0.06 0.005 0.01 
CrAs_ 0.036 0.05 0.05 0.007 0.006 
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Table 1.4. Positions of atoms in the MnP-type structure expressed in the coordinates based on 
A-, B-, and C-axes 


M (—x,0,u), (x, 1/2,1/2—u), (1/2—x,0,—u), (1/2 +x, 1/2,1/2+u) 
X (1/4— 1,1/2,5/6—w), (1/4+v,0,1/3—w), 
(3/4—v,1/2,1/6+w), (3/4+v,0,2/3+w) 


@ Cu(D: (0, 0, 0), (1/2, 1/2, 0) 

© Cu CD: (1/2, 0, u), (, 1/2, 1—-w) 
O Sb: (0, 1/2, v), (1/2, 0, 1-v) 
u=0.295 

v=0.280 


a (4.08 A) 


Fig. 1.5. Unit cell of the Cu2Sb-type structure. 3d transition metal occupies the Cu(I)- and 
Cu(I])-sites. In MnAlGe, Al atoms occupy the Cu(II)-sites and Ge occupies the position of Sb. 
The crystallographic parameters u and v are these of Mn2Sb 


reliable accuracy by X-ray diffraction for powdered samples. Table 1.4 shows posi- 
tion of atoms (see Figs. 1.3 and 1.4). We note that the coordinates are based on the 
A-, B-, and C-axes (c-, a-, and b-axes) shown in Fig. 1.4. 

In some literatures, origin of coordinate for the atomic position is shifted to 
(1/4, 1/4) in the ab (BC)-plane shown in Fig. 1.3 (or 1.4). In this case, the positions 
shown in Table 1.4 must also be shifted. There have been several (about four) man- 
ners of expression for the lattice parameters of the MnP-type structure. The (a, b,c)- 
system used in Figs. 1.3 and 1.4 is useful because the a- and c-axes of the NiAs-type 
lattice correspond to the a- and c-axes of the MnP-type structure. We use mainly 
(a,b,c) system and sometimes (A,B,C) system in this book. 


1.3.3 Cu2Sb-Type Structure 


MnAIGe, MnGaGe, and MX denoted by boldface letters in Table 1.1 take this crys- 
tal structure, whose unit cell is shown in Fig. 1.5. In this structure, distance between 
the Mn(II) atoms (denoted by the thick dotted line) is large. For example, the distance 
is 3.94A for MnoSb. Nearest distance between Mn(1)—Mn(II) and Mn(1)—Mn (J) is 
2.81 and 2.88 A, respectively. In the case of MnAIGe, which is obtained by substi- 
tuting Mn(II) and Sb of Mn2Sb by Al and Ge, respectively, Mn—Mn distance along 
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the c-axis is very large and equal to c(= 5.93 A). This is also the case of MnGaGe. 
Crystallographic characteristics mentioned above affects magnetic properties, which 
we discuss in Chap. 3. 


1.4 Nonstoichiometric Compounds 


One of typical nonstoichiometric MX compounds is Mn; +,Sb. Figure 1.6 is a phase 
diagram of Mn-—Sb binary system [3]. The NiAs-type structure exists at room tem- 
perature for Sb composition between 44% and 49%. In other words, Mn;+,Sb with 
0.04 < x < 0.25 has the NiAs-type structure. Curie temperature Tc of Mnj,+,Sb 
decreases linearly from 314°C (51 at% Mn) to 90°C (56 at% Mn) with increasing 
Mn composition. This diagram shows that stoichiometric (i.e., 50 at% Mn) MnSb 
does not exist. MnSb is considered to have about 1% excessive Mn atoms occu- 
pying partially the interstitial sites ((2/3, 1/3, 1/4) and (1/3, 2/3, 3/4)) shown in 
Fig. 1.7. It has been reported that an excessive Mn atom at the interstitial site is 
surrounded by six nearest Mn atoms and does not have magnetic moment, while 
an Sb atom polarizes with small moment 0.2 antiparallel to the total moment [5]. 
For Mn, +,Sb (0.04 < x < 0.25), Tc and spontaneous magnetization rapidly decrease 
with increasing x. Magnetic anisotropy constants also vary with x [6]. 

Some other MX compounds that forms the NiAs- (or MnP-) type structure in 
nonstoichiometric composition region are shown in Table 1.5. 
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Fig. 1.6. A phase diagram of the Mn-Sb binary system. The label 2P denotes two phases 
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Fig. 1.7. Interstitial positions of the NiAs-type structure at which excessive Mn atoms of 
Mn;,Sb locate are denoted by the filled circles 


Table 1.5. Composition region where MX-type compound exists (room temperature) 


VSb Single phase only for x = 0.4 [7] 

CrSb 0 <x < 0.04. Upper limit of x increases with temp [1,3] 

MnSb _ See Fig. 1.6 

FeSb 0.2 <x < 0.3 [8], 0.1 <x < 0.3 [3] Region incr. with temp [8] 
Exists only for excessive Fe 

CoSb O<x< 0.04. Region incr. with temp [1]. —0.04 < x < 0.15 [3] 

NiSb —0.06 < x < 0.02. Region incr. with temp [1]. —0.08 < x < 0.10 [3] 

VAs 1:1 compound 

VP 1:1 compound 

CrP 1:1 compound 

CrAs 1:1 compound 

MnP 1:1 compound 

MnAs _ 1:1 compound 

MnBi MnBi for T < 340°C. Mn, ogBi for T > 340°C 
Phase separation into Mn and Mn-Bi liquid for T > 446°C [3, 8] 


The letter x in this table means composition in M;,,X. Generally, the region 
expands as temperature rises 
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Overview of Magnetic Properties of NiAs-Type 
(MnP-Type) and Cu,Sb-Type Compounds 


2.1 Compounds That Have Magnetic Ordering Phase 


In this section, we look over structural and magnetic properties of the MX and M2.X 
compounds. In Table 2.1, we list the NiAs- (or MnP-) type and the Cu2Sb-type com- 
pounds that show magnetic ordering. Basic magnetic properties of typical MX and 
M>2X compounds are shown in Tables 2.2—2.6. In the tables, only representative data 
are shown when there are several data. Refer to Chap. 3 for more detailed magnetic 
properties. 

As shown in Tables 2.1—2.4, magnetic ordering occurs only in compounds con- 
taining Cr, Mn or Fe. The compounds with the MnP-type structure show the double 
helical magnetic ordering as shown in Table 2.6. 

As a typical example of the double helical structure, we explain the magnetic 
structure of CrAs in Fig. 2.1. 

We summarize the double helical magnetic orderings of the MnP-type com- 
pounds in Table 2.6. MnP has a long wavelength and thus is close to a ferromag- 
net, and is known to show metamagnetic magnetization process. Magnetization of 
MoP saturates at magnetic field about 4kOe when the applied field is along the 
a-axis. There are no magnetization measurements for the other compounds listed 
in Table 2.6. There is only an unpublished data by us for powdered sample of CrAs 
measured in pulsed fields up to 35 T at liquid nitrogen temperature. According to 
this report, magnetization shows no saturation and reveals anomaly at 10 T. CrAs is 
an interesting compound with various unsettled problems. We will discuss some of 
these problems in Chap. 3. 

MnxX (X = As, Sb, or Bi) are ferromagnets. CrAs, CrSb, MnAs, and MnBi show 
a first order phase transition at Tc (or Ty). See Chap. 3 for the detail of these phase 
transitions. 

The CuySb (M2X)-type compounds have rich variety of magnetic structures 
as shown in Fig. 2.2 [1]. In particular, ferrimagnetic Mn»Sb is interesting; Slight 
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Table 2.1. The shape of the signs indicates crystal structure at room temperature 


x \M Ti Vv Cr 


P MP O QO 


Fe Co Ni 
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Compounds that show magnetic ordering are denoted by the gray sign 


Table 2.2. Magnetic properties of CrX [1] 


Mx CrP CrAs CrSb 
Magnetic structure yhasapeak~ 200K Doublehelix Antiferro 
Ty(K) 250 (Ist) 718 (Ist) 
DA (Ltp/Cr) x(RT ~ 3) x 10-6 1.67 3.0 


Unit of magnetic susceptibility 7 is emu g—!. See also Chap. 3. Ty and 
Pa are Néel temperature and magnetic moment per atom. (lst) means 
first order transition 


Table 2.3. Magnetic properties of MnX [1] 


Mx MnP MnAs MnSb- MnBi 

Magnetic structure Double helix (T < 47 K) Ferro Ferro Ferro 
Ferro (47K < T < 291K) 

Tc (K) 291 318 (ist) 587 628 (Ist) 

Da (Ltp/Mn) 1.3 3.4 3.6 3.8 


See also Chap. 3. Tc is Curie temperature. (1st) means first order transition 


substitution of Cr etc. for Mn induces a first order phase transition from an anti- 
ferromagnetic to a ferrimagnetic state with increasing temperature. We will discuss 


this phenomenon in Chap. 3. The various magnetic ordering shown in Fig. 2.2 is also 
discussed in Part II of this book. 
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Table 2.4. Magnetic properties of FeX [1] 


MX FeP FeAs Fe,,,Sb 

Magnetic structure Double helix Doublehelix Triangular in c-plane 
Tn(K) 125 77 105 ~ 211 (depends on x) 
Pa (Up/Fe) 0.41 0.51 0.88 


Table 2.5. Magnetic properties of M2X (See [1] and Chap. 3) 


MoX Cro As Mn)As Fe2As Mn2Sb 
Magnetic structure Antiferro  Antiferro Antiferro Ferri 
Tn(K) 393 573 353 550 


Pa (Up/M),.M(D) 0.40, 1.34 2.24.1 —-1.28,2.05 2.13, 3.87 


Table 2.6. Double helical structure 


MnP CrAs FeP FeAs 


O13 16° 126° 169° 140° 
2n/k(=2) 8.2b 2.67b 5.0b  2.67b 


See Fig. 2.1 for 0,2. b is length of the b-axis. 
Data are taken at 4K (12 K for FeAs) (see [1] 
and Chap. 3) 


ae 


A =27/k=2.67 b 


Fig. 2.1. Double helical magnetic ordering of CrAs. (a) Cr atoms projected onto the ab-plane 
(see Fig. 1.3b). The spins are parallel within the plane 1, 2, 3,... normal to the b-axis. (b) As 
the plane changes 1, 3, 5,... (or 2, 4, 6,...), the direction of the spin rotates by 67.4° with 
the relative angle 126° between adjacent (i.e., between 1 and 2) planes. The wavelength is 
A = 2.67b and Uc; = 1.67 Up [2]. Slightly different data are reported in [3]. See also Fig. 4.33 
in Part I 


2.2 Compounds without Magnetic Ordering 


The compounds denoted by unshaded symbols in Table 2.1, which do not include 
Cr, Mn, or Fe, are paramagnetic or diamagnetic. We summarize properties of these 
compounds that show no magnetic ordering in Table 2.7. 


14 2 Overview of Magnetic Properties of NiAs-Type and Cu2Sb-Type Compounds 


O@M(® 


O O 
wa e o ee Oo @ Oo 
® MUD 


OQ g 


*®@ <O +@ <0 e © 
iv Lieely 
® © 
unit 
cell 
i Zoe Zz eo 8 © & o 
Cr.As Mn.As Fe.As Mn.Sb 
(antiferro) (antiferro) (antiferro)  (ferr1) 


Fig. 2.2. Magnetic orders of the Cu2Sb-type compounds. M(I) and M(ID (see Fig. 1.5) have 
different moments. As (Sb) is not shown [1] 


Table 2.7. Properties MX compounds that have no magnetic ordering 


Mx Struct. Magnetism (unit of v is emu g~!) 


TiP TiP 
TiAs  TiP Pauli para. y ~ 1.5 x 10-6 (T =4.2K) [4] 
TiSb NiAs Pauli para [5] 


VP NiAs x slightly decreases with temp. 7 ~ 2 x 107° (RT) [6] 
VAs MnP x slightly increases with temp. 7 ~ 1.5 x 10-6 (RT) [6] 
VSb NiAs Single phase in V; 4Sb. Pauli para. y ~ 1.5 x 10-6 [7] 
CrP MnP x has a broad valley at 200K. y ~ 3 x 10-6 (RT) [8] 


CoP MnP Pauli para. y¥ ~ 1.5 x 10-6 (T > 100K) [9] 
CoAs_ NiAs x has a broad peak at 200K. y ~ 2 x 10-6 (RT) [10] 
CoSb_ NiAs Pauli para. y ~ 1 x 10-6 [4] 
NiP NiP 

NiAs_ NiAs Pauli para. y ~ 2.5 x 10-7 (T > 300K) [11] 
NiSb  NiAs Diamag [12] 
NiBi  NiAs Superconducting, T¢ = 4.25 K [13] 


MX compounds that have no magnetic ordering are of little interest in mag- 
netism. However, we can investigate genuine electronic structures that are helpful 
to understand properties of other MX compounds. In Table 2.7, some of such com- 
pounds with weak temperature dependence of susceptibility are denoted by “Pauli 
para”. The susceptibility v of the order of 10~° to 10-7 emu g! may be considered 
by not only 3d states but also diamagnetism of anion, and possibly by some kind of 
magnetic impurities. Some compounds such as VSb has a nonstoichiometric phase 
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in which the composition deviates from 1:1. Magnetic properties are often affected 
by nonstoichiometry. As already described in Sect. 1.4, ferromagnetic Mn) +,Sb is a 
typical example of nonstoichiometric compounds [7, 14, 15]. 


References 


[1] 


[2] 
[3] 


[4] 
[5] 
[6] 
[7] 
[8] 


[9] 
[10] 
[11] 


[12] 
[13] 
[14] 
[15] 


Landolt-Bornstein III/27a, Magnetic Properties of Pnictides and Chalcogenides, ed. by 
K. Adachi, S. Ogawa (Springer Berlin, 1989), p. 70 and references therein 

G.P. Felcher, F.A. Smith, D. Bellavance, A. Wold, Phys. Rev. B 9, 3046 (1971) 

K. Selte, A. Kjekshus, W.A. Jamison, A.F. Andresen, J.E. Engebresen, Acta Chem. 
Scand. 35, 1042 (1971) 

H. Ido, J.Appl. Phys. part ITA 3247 (1985) 

K. Adachi, J. Phys. Soc. Jpn. 16, 2187 (1961) 

K. Selte, A. Kjekshus, A.F. Andresen, Acta Chem. Scand. 26, 4057 (1972) 

J. Bouma, C.F. van Bruggen, C. Haas, J. Solid State Chem. 7, 255 (1973) 

K. Selte, H. Hjersing, A. Kjekshus, A.F. Andresen, P. Fischer, Acta Chem. Scand. A29, 
695 (1975) 

K. Selte, L. Birkeland, A. Kjekshus, Acta Chem. Scand. A32, 731 (1978) 

K. Selte, A. Kjekshus, Acta Chem. Scand. 25, 3277 (1971) 

LL.A. Delphin, K. Selte, A. Kjekshus, A.F. Andresen, Acta Chem. Scand. A32, 179 
(1978) 

H. Schmit, Cobalt 7, 26 (1960) 

N.E. Alekseevskii, N.B. Brandt, T.I. Kostina, Izu. Akad. Nauk SSSR 16, 233 (1952) 

Y. Yamaguchi, H. Watanabe, T. Suzuki, J. Phys. Soc. Jpn. 45, 846 (1978) 

T. Okita, Y. Makino, J. Phys. Soc. Jpn. 25, 120 (1968) 


3 


Properties of the Compounds with NiAs-Type 
(MnP-Type) and Cu,Sb-Type Structures 


3.1 MnP and Related Compounds 


A magnetic phase diagram of MnP with respect to magnetic field and tempera- 
ture is shown in Fig. 3.1. This diagram was obtained from magnetization measure- 
ments for a single crystal MnP [1]. It has been known that in the temperature region 
T < 47K, MnP shows a double helical ordering with relatively long wavelength 
of 8.2b propagating along the b-axis [6]. Above 47 K, MnP is a ferromagnet with 
Tc = 291K. For 420K < T < 660K, susceptibility obeys the Curie—Weiss law with 
Mert = 2.36hg (2S = 1.56) and 6p = 344K. This value of Herr is close to saturation 
moment |.3g observed at T = 0K [7]. The ah curve shows upward convex in 
the range 660K < T < 1400K [8]. Studies on mixed compounds with the type of 
Mn,_,M,P (M = other 3d metal) also have been reported [9-12]. 


3.2 MnAs and Related Compounds 


3.2.1 Magnetic Transition of MnAs and the Bean—Rodbell Theory 


MnAs is especially interesting among the MX compounds. In this section we show 
magnetic and crystallographic properties of MnAs and give a simple explanation for 
them on the basis of the phenomenological theory developed by Bean and Rodbell. 
We also mention some interesting properties of mixed compounds between MnAs 
and MnP (and MnSb).! 

MnAs has long been of interest since the study carried out by Guillaud et al. 
[13]. Temperature dependence of magnetization and inverse susceptibility are shown 
schematically in Fig. 3.2 [14]. With rising temperature, spontaneous magnetization 
disappears discontinuously at Tc(up) = 318K. This first order transition is accom- 
panied with structural transformation from the NiAs-type to the MnP-type structure. 
Susceptibility shows a peak at 7, = 398 K and obeys the Curie—Weiss law for T > T. 


' See Part II of this book for the itinerant theory on magnetism of MnAs. 
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Fig. 3.1. Magnetic phase diagram of MnP [1-4]. See also [5] for more detailed diagram 
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Fig. 3.2. Spontaneous magnetization o, and inverse susceptibility y~! of MnAs plotted as 
functions of temperature 


As we will see later, 7; is a structural transition temperature from the MnP-type to 
the NiAs-type structure. Anomalous temperature dependence of paramagnetic sus- 
ceptibility in the temperature range 7c < T < T, is explained from the viewpoint of 
itinerant electrons (see Sect. 4.3). 

Figure 3.3 shows temperature dependence of the crystallographic parameters 
[15]. The a-axis and the unit cell volume show remarkable discontinuous changes 
of 1.2% and 2.1%, respectively, at Tc(up), while the c-axis does not show visible 
change at Tc(up). These facts suggest that there is a correlation between the a-axis 
and ferromagnetism of MnAs. Thermal expansion coefficient has a large value in 
the MnP-type phase just above 7c. This is due to the rapid structural transformation 
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Fig. 3.3. (a) Temperature dependence of lattice constants of MnAs. (a, b, c are lattice constants 
of the hexagonal lattice.) @ is the thermal expansion coefficient. (b) Temperature dependence 
of the parameters u, v, and 6. Ymax means the temperature where the v —!_T curve shows the 
minimum as shown in Fig. 3.2 (Reproduced from [15]) 


from the MnP-type to the NiAs-type in the small range of temperature (see Part II for 
detail). The crystallographic parameters u, v, and 6 (see Table 1.3) take the maximum 
value just above Tc(up), approaching zero at 7; on further increasing temperature as 
indicated in Fig. 3.3b. 

MnAs shows a metamagnetic magnetization process in the temperature range 
Tc < T <T,, corresponding to a paramagnetic—ferromagnetic transition accompanied 
with a crystallographic change from the MnP-type to the NiAs-type structures. As 
the magnetization process is similar to that of MnAs;_,P,, which we will discuss 
later, we show in this section only temperature change of magnetization in H = 2 T 
and H = 40T in Fig. 3.4 [16]. The curve in H = 40T suggests that, if the NiAs-type 
structure below 7c were maintained up to the temperature region just above 7c of 
the spontaneous magnetization, Tc could be about 450 K. It is noted that this roughly 
extrapolated Curie temperature Tc ~ 450 K is considerably higher than 6, ~ 270 K 
estimated in the high-temperature NiAs-type region, which is characteristic of MnAs 
as Tc ~ @p is satisfied for MnSb and MnP. 

Various characteristic properties of MnAs mentioned above are well explained 
by a phenomenological theory based on an assumption of strong volume dependence 
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Fig. 3.4. Temperature dependence of magnetization of MnAs in two magnetic fields H = 2 
and 40 T. Tc means Tc(up) (Reproduced from [16]) 


of exchange interaction between Mn atoms [17,18]. Namely, the first order transition 
of MnAs (such as the discontinuous change of crystal volume and spontaneous mag- 
netization, and their hysteresis for temperature variation), metamagnetic phase tran- 
sition in the temperature region just above 7c (e.g., magnetization process and tem- 
perature dependence of hysteresis), and pressure dependence of 7c can be generally 
explained. We have to note, however, that 3d electrons in MnAs from which magne- 
tization arises have itinerant character. Therefore, it is necessary to treat 3d electrons 
in MnAs as itinerant electrons as we will show in Part II of this book. In this section 
we describe the phenomenological theory to show how the various peculiar proper- 
ties of MnAs are explained. 

First, we assume that the Curie temperature depends on crystal volume (or lattice 
constants) [17]. This assumption is reasonable as Tc depends on applied pressure 
in most ferromagnets. The crystal volume is affected by various factors, for exam- 
ple, thermal expansion, pressure, applied magnetic field, and magnetovolume effect 
(magnetostriction). Structural phase transition also changes crystal volume. Consid- 
ering such factors, Bean and Rodbell [17] set a relation 


Tc(V) = To[1 + B(V —Vo)/Vol, (3.1) 


where 7c is assumed to be a linear function of crystal volume V. This may be justified 
as long as volume change is small. In (3.1), Vo is the volume at T = 0 K in the absence 
of magnetovolume effects. Jj is an expected Curie temperature when V is assumed 
to be kept in the constant Vo. Volume dependence of Tc is given by the coefficient B, 
which takes a large value for MnAs. 

We consider a localized model where a magnetic material consists of N magnetic 
atoms with spin j. Under an applied magnetic field H and pressure P, free energy of 
this material is written in the form 
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G =—(3/2)[i/(i + 1)|NkpTco? — HM,o + (1/2K)[(V —Vo)/Vo]? 
—TkgNS + P[(V —Vo)/Vol, (3.2) 


where o(= M,(T)/M,(0)) is a magnetization normalized by Ms(0) specific magne- 
tization at OK. The first term denotes exchange energy written by using 7c.” The 
second and the third terms are Zeeman energy and elastic energy (K is the com- 
pressibility), respectively. The fourth term is the entropy one and kgNS (kg is the 
Boltzmann constant) is the entropy, where S is summation of magnetic contribution 
So(o,j) and lattice contribution S,. The last term denotes energy due to applied 
pressure. 

When o is smaller enough than 1, Sg of the magnetic entropy term kgNSg(o, /) 
can be expanded in the following form: 


So(6,j) =So(0, j) +ao* + bot +co°+..., (3.3) 
where the coefficients a, b,... are given by 
a= —(3/2)[//(i+1)], b= —(9/20)[(2+ 1)*- Y/RG+ DF. (3.4) 


See Appendix for the detailed expression of (3.4). This expansion is useful to in- 
vestigate magnetic phase transition. We note that the exact expression of Sg(o, j) is 
given by 


So(0,j) —So(6, jf) = ao — [ B(a) der, (3.5) 


where B (0) is the Brillouin function with angular momentum j. For the discussion 
in this section, we use (3.3) and (3.6), which is derived from (3.5) as given by 


OSo/00 = —@. (3.6) 


Substituting (3.1) into (3.2), and using a condition of minimum free energy 
dG/OV = 0, we have a relation 


(V —Vo)/Vo = (3/2) [77 (J + 1)|NkpKToB 0? — PK. (3.7) 


This equation means that magnetovolume effect (or exchange striction) is propor- 
tional to 0? if P=0. 

First, we consider a case of P = 0. Inserting (3.7) into (3.2), we obtain an expres- 
sion of G, which has been minimized with respect to V. For this free energy, we use 
a condition dG/do = 0 and (3.6), and then we have a relation among spontaneous 
magnetization o(= M,(T)/M,(0)), temperature 7, and a magnetic field H as 


(g.jMe /kp)(H/T)) = 2ao[1 + (2b/a)no*] + (T/To)o +, (3.8) 


? Exchange energy is written as —(1/2)AM2 with the coefficient A of molecular field. The 
first term in (3.2) is derived from this expression. 
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where [lp is the Bohr magneton and g is the g factor (=2 in the present case). In this 
equation 77 is defined by 


n = —(a’/2b)NkgKTB°. (3.9) 


We note that 7 is proportional to B*. As we will show below, 7) is an important 
parameter contained in the coefficient of the o* term of the free energy expansion. 
For 7) > 1, a first order transition occurs, while order of transition is second for n < 1 
(see also Appendix). The coefficients a and b are given in (3.4), and q@ is a function 
of 6 (i.e., 6 = Bj(@)). 

Using (3.8) with a given value for j and setting H = 0, we can calculate tem- 
perature variation of spontaneous relative magnetization o = o(T/Tp) for various 
values of 7. We can also calculate magnetization curve by fixing temperature and 
changing H. By considering that magnetization of MnAs at T = 0K is 3.4Up, we set 
j =3/2 (with g = 2). Figure 3.5 shows temperature dependence of o calculated from 
(3.8) for various values of 1. Magnetization curves calculated for 7 = 2 at various 
temperatures are plotted in Fig. 3.6. 

We see from Fig. 3.5 that 7 = 1 is the critical value between the first and the sec- 
ond order transition. For 1 = 2, relative spontaneous magnetization o is 0.81 att = 1 
and is 0.6 at t = 1.08 (i.e. just below 7c). These results agree well with experimental 
values shown in Fig. 3.4. We also see that magnitude of hysteresis of spontaneous 
magnetization is A(T /Ty) = 0.08. This value is also close to the experimental value 
A(T /Tp) = 0.06 (i.e., Tc(up) = 318 K and Tc (down) = 7p = 301 K). Therefore, tem- 
perature dependence of spontaneous magnetization of MnAs is well reproduced with 
the parameter 7 = 2. 


1.5 
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TIT, 
Fig. 3.5. Relation between magnetization (o) and normalized temperature (T /7p) for various 


values of 1. The curves were calculated on the basis of (3.8). The curve with n = 2 is in good 
agreement with experimental curve for MnAs 
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Fig. 3.6. Calculated magnetization curves of MnAs at various temperatures (T/T) = t) with 
1 = 2 fixed. We see that spontaneous magnetization is o = 0.81 att = 1. The curve at t = 1.09 


shows metamagnetic behavior clearly but the curves for ¢ > 1.134 do not show hysteresis as 
well as magnetization jump 
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Fig. 3.7. Temperature dependence of Hc(up) and Hc(down) of metamagnetic phase transition 
calculated for MnAs using 7 = 2 


It is known that MnAs shows metamagnetic behavior at temperatures higher than 
Tc(up) [19]. Figure 3.6 is the result calculated from (3.8) with n = 2. The magnetiza- 
tion curves for 1 < T/Tp < 1.134 show a jump at the transition field Hc. These curves 
also show hysteresis as indicated by the vertical dotted lines for t = 1.09. Figure 3.7 
shows the temperature variation of transition fields for the process of increasing and 
decreasing field evaluated from Fig. 3.6. 
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Setting Ty) = T-(down) = 301K, the critical point (T/T = 1.13 and Hc/Tp = 
0.029) in Fig. 3.7 corresponds to T = 340K and H = 8.7 T. These values agree well 
with experimental results T = 345 K and H = 9T found by Zieba [20]. Magnetiza- 
tion curves measured by Grazhdankina [19] show that the critical point locates near 
350 K, which is also in good agreement with the calculated results mentioned above. 
We note that Fig.3.7 is also applicable to metamagnetic behavior of MnAs,_,P, 
above Jy shown in the next section. 

In summary, anomalous magnetic properties of MnAs is well explained by a 
phenomenological theory in which exchange interaction between Mn is sensitive to 
crystal volume (that is, 7) or f is large). 

Next we apply this theory to the magnetovolume effect of MnAs. We start 
from (3.7). According to experiments, the jump of magnetization at the tempera- 
ture T/T) = 1.08 is about 0.6. This value is reproduced with n = 2. By inserting 
the value of T/Tp = 1.08 and AV /V = 0.021 (obtained from experiments) shown in 
Fig. 3.3 into (3.7), we have 


NkgKToB = 0.0648. (3.10) 
Also, a relation 
N =2=2.21 x NkgKToB? (3.11) 


is known from temperature dependence of magnetization. From (3.10) and (3.11), we 
obtain B = 14 and K =3.8 x 10~!? (in CGS). Temperature dependence of magnetiza- 
tion as well as magnetovolume effect at T = Tc(up) (AV/V = 0.021) is reproduced 
using these values. However, the experimental value of AV/V is a summation of 
magnetovolume effect and volume change due to structural transformation from the 
NiAs-type to the MnP-type. To investigate pure magnetovolume effect, it is neces- 
sary to investigate MnAs;_,Sb,, which does not show the structural transformation 
for x > 0.1 [21,22]. We note that the magnetovolume effect of MnAso.9Sbo.; is well 
explained with the critical value 7 = 1. 

Finally, we investigate the effect of pressure on the magnetic transition temper- 
ature 7) of MnAs. In the procedure to derive (3.8), we set H = 0 and keep P finite 
instead of setting P = 0. Then we have a relation 


0 = 2ao[(1 — BKP) + (2b/a)no’| + (T/T) a. (3.12) 


This equation shows a relation between o and T /Tp under pressure P. We set ) = 2 
and 7) = 301K as we did for MnAs. The value of BK is evaluated to be BK = 
5.32 x 107! from (3.10) and (3.11). Using these values of parameters, we calculated 
temperature dependence of o under various pressures as shown in Fig. 3.8. From 
this figure, 7-¢(up) and 7-(down) are obtained as a function of pressure as plotted in 
Fig. 3.9. 

The calculated results shown in Fig. 3.9 agree qualitatively with the pressure de- 
pendence of 7c¢(up) and 7-(down) in the low pressure region in Fig. 3.10. We see 
from this figure that MnAs under pressure transforms easily to the MnP-type struc- 
ture. We also see that this transformation is of the first order accompanied with large 
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Fig. 3.9. Pressure dependence of first order transition temperatures Tc(up) and Tc(down) of 
MnAs evaluated from Fig. 3.8 


hysteresis for pressure. For example, if we apply pressure larger than P = 3.1 kbar 
(0.31 GP) at T = 100K, the crystal remains in the MnP-type even after the pressure 
is removed. This phase diagram is similar to that of MnAs;_,P; (Fig. 3.12), which 
we will discuss in the next section. These results suggest that NiAs-type structure of 
ferromagnetic MnAs is quite unstable; it seems to be maintained by occurrence of 
the ferromagnetic order. 

The model we have shown in this section is based on assumptions that Mn atom 
in MnAs has a local spin 3/2 and that the exchange interaction depends strongly 
and linearly on the crystal volume. This phenomenological model has explained 
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Fig. 3.10. A phase diagram (pressure-temperature) of MnAs (Reproduced from [23]) 


successfully characteristic properties of MnAs. This model, however, has defects. 
One is that it is based on the localized spin model. Another is that the structural 
transformation at Tc¢(up) is neglected. In Part II of this book, we will show that the 
anomalous magnetic properties and structural transformation of MnAs are well ex- 
plained from the viewpoint of itinerant electrons and band structures. 


3.2.2 Various Type of Phase Transition of MnAsj_,P.. 


There are various interesting mixed compounds of MnAs such as Mn;_,Cr,As 
[24-27], MnAs,_,P, [15,28-30], and MnAs,_,Sb, (see Sect. 3.2.3) [5]. As we have 
noted, MnAs is unstable against the MnP-type lattice distortion. Substituting P for 
As realizes the MnP-type structure. On the other hand, by substitution of Sb for As, 
the NiAs-type structure is stabilized. We describe MnAsj_,P, in this section. 

First we show a relation between magnetic moment and the a-axis of MnX 
(X =P, As, Sb, or Bi) in Fig.3.11. In this figure, MnP and MnAso.8gPo.12 takes 
the MnP-type structure. (We note that the a-axis of the MnP-type structure in this 
figure corresponds to that of the NiAs-type structure.) We see that magnetic moment 
is sensitive to P-substitution, but the moment does not change largely as long as the 
a-axis is longer than that of MnAs. 

Figure 3.12 shows a phase diagram of the MnAs—MnP system, which is de- 
picted using data in [31-34]. In this figure, the hatched region indicates the region 
in which the NiAs-type structure realizes. The compound in the vicinity of MnAs 
(x < 0.03) is a ferromagnet, which shows a first order transition at 7c. In the region 
of the MnP-phase, the double helical ordering realizes although thorough investiga- 
tion has not been performed yet. The first investigation of these mixed compounds 
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Fig. 3.11. Relation between magnetic moment and the a-axis of MnX compounds. The NiAs- 
type and the MnP-type structure are denoted by filled and open symbols, respectively 
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Fig. 3.12. A phase diagram (magnetic state and crystal structure) of the MnAs—MnP mixed 
crystal 


were carried out by Goodenough et al. [32]. In Fig.3.13, we show M-T curves of 
MnAso.9Po.1. The Néel temperature is about 230 K. At the temperatures T < 100K 
the double helical ordering realizes. By applying a magnetic field, the double heli- 
cal structure transforms into ferromagnetic ordering with moment 1.6Up/Mn [35]. 
In the temperature range 100K < T < Ty, magnetic structure is unknown though 
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Fig. 3.13. Magnetic moment [35] and inverse susceptibility of MnAso.9Po.; plotted against 
temperature. erp is effective Bohr magneton (Reproduced from [32]) 
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Fig. 3.14. Magnetic moment of Mn in MnAs,_,P,. The data for x = 0, 0.1, 0.25, and 1 were 
obtained from saturation moment measurement [14, 35-37]. Others were taken from results of 
neutron diffraction measurements [28, 31] 


ferromagnetic behavior of a magnetization curve has been observed. Magnetic mo- 
ment of MnAs;_,P, was measured by neutron diffraction [28,31] and saturation 
moment measurements [16, 35]. The moment slightly decreases with increasing P 
component for x > 0.03 as shown in Fig. 3.14. 
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As we see from Fig. 3.13, behavior of y~! above Ty is similar to that of MnAs 
above Tc(up). In the temperature range where y~! decreases with rising temperature, 
the MnP-type structure is changing into the NiAs-type (the parameters u, v, and 6 
approach zero) [15,36], and the NiAs-type structure realizes above 7;. Moment 
estimated from the Curie—Weiss-like behavior of y~! is Leg = 2.09UpB (pa = SUB = 
1.3p). This is close to the saturation moment observed at low temperatures shown 
in Fig. 3.14. In the range T > Ti, Ucer is 4.87 Up, indicating that magnetic moment is 
as large as the moment of MnAs. 

Figure 3.13 suggests that transition between a low spin state and a high spin 
state occurs accompanied with the structural transformation. This transition was dis- 
cussed first by Goodenough et al. [32]. There is an experiment that indicates relation 
between susceptibility and crystal structures [15,36]. Such problems have been ex- 
plained from the viewpoint of itinerant electrons by Motizuki and her co-workers as 
described in Part II of this book. 

Next we describe field-induced structural transformation of MnAso.9Po.1. We 
show temperature dependence of moment of MnAsp 9Po,; in high magnetic fields in 
Fig. 3.15 [16]. As shown in Fig. 3.16, metamagnetic behavior is seen above Ty. This 
is similar to the behavior of moment of MnAs just above 7Jc(up). The data shown in 
Figs. 3.15 and 3.16 were measured at the high magnetic facility of Osaka University. 
Behavior similar to Figs. 3.15 and 3.16 have been observed for Mno.9Cro.9As [37]. 

In a moderate magnetic field, MnAso 9Po.; has the double helical ordering with 
low spin state below Ty ~ 230K. In high-fields, metamagnetic transition to the 
state similar to MnAs occurs for T > Ty as shown in Figs.3.15 and 3.16. This 
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Fig. 3.15. Temperature change of magnetization of MnAso.9Po,1 under high magnetic fields. 
Data of MnAs are also plotted for comparison (Reproduced from [16]) 
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Fig. 3.16. Magnetization process (a) and transition fields (b) of metamagnetic behavior of 
MnAso 9Po.1. The open (filled) symbols are for increasing (decreasing) field (Reproduced 
from [16]) 
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Fig. 3.17. Schematic figure for magnetization-temperature curves of MnAs and MnAso 9Po.1 
under high magnetic fields. Temperature range for the NiAs-type and MnP-type crystal struc- 
ture is also shown 


field-induced phase transition is expected to be a transition from paramagnetic state 
(MnP-type structure, the low spin state) to ferromagnetic state (NiAs-type, the high 
spin state). This transition is similar to the metamagnetic transition of MnAs just 
above Jc, as we have noted. Figure 3.17 shows a comparison between MnAs and 
MnAspo.9Po.1- 

The fact that the magnetization of MnAsg9Po,; is almost the same as that of 
MnAs in the temperature region between Ty (=250K) and Tc (=318K), as well 
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Fig. 3.18. Magnetostriction of MnAso.93Po,97 measured using a pulse field at T = 300 K [38] 


as very large magnetostriction, which will be explored next, indicates that the 
metamagnetic transition of MnAso.9Po.1 above Ty accompanies the structural trans- 
formation. In MnAso 93P0.07, very large magnetovolume effect AV/V = 0.15 accom- 
panied with the metamagnetic transition has been observed at T = 300K, as shown 
in Fig. 3.18 [38]. 

The magnetovolume effect shown in Fig. 3.18 is extraordinarily large but rea- 
sonably understood. Temperature dependence of lattice constants of MnAsj_,P, is 
plotted in Fig. 3.19 [15]. Unit cell volume V(= abc) for x = 0.1 at 300K is 121 A?, 
while that of MnAs in the NiAs-type phase is 137 A? at 300K (Fig. 3.3). Thus, there 
is 13% difference (as also shown in Fig. 3.17), which suggests that the large magne- 
tovolume effect shown in Fig. 3.18 is reasonable. Therefore, the metamagnetic tran- 
sition shown in Fig. 3.16 is accompanied with the structural transformation (between 
the MnP-type and NiAs-type) with extraordinary volume change more than 10%. 
Metamagnetic transition accompanying large magnetostriction is possible to occur, 
because increase of elastic energy due to strain is small in the present case (this is 
because the very large thermal expansion coefficient of MnAs,_,P, results in a small 
elastic coefficient). In Fig. 3.16, we see that transition field increases with decreasing 
temperature, which is closely related to the decrease of thermal expansion coefficient 
as seen in Fig. 3.19 [15]. The thermal expansion coefficient of MnAso.9Po.1 is small 
in the temperature range T < Ty and the transition field of the metamagnetic transi- 
tion shown in Fig. 3.16 shows rapid increase. At T = 231 K, metamagnetic transition 
does not occur even in 400 kOe. Results shown in Fig. 3.16 suggests that no meta- 
magnetic transition occurs for T < Ty. These facts seems to be important hints to 
understand the field-induced metamagnetic transition, which occurs only above Ty. 
The phenomenological theory of metamagnetism of MnAs just above 7c shown in 
Sect. 3.2.1 may be applicable to MnAso 9Po,;. Divergence of transition field at Ty 
(Fig. 3.16) and metamagnetic transition of MnAso.9Po.; accompanied with structural 
change are complicated and interesting phenomena. Theoretical studies based on the 
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Fig. 3.19. (a), (b) Lattice parameters of MnAs,_,P,. a and c are the lattice constants of the 
hexagonal lattice. (c) Unit cell volume V = abc of MnAsg.9Po,; and MnAs is schematically 
plotted as a function of temperature (Reproduced from [15]) 


itinerant electron picture, as well as experiments in an ultra high field, are necessary 
to clarify the mechanism of metamagnetic transition in MnAsj_,P,. 


3.2.3 Anomalous Behavior of MnAs,_,Sb, 


The NiAs-type structure is stabilized by substituting Sb for As of MnAs. Figure 3.20 
shows a magnetic and structural phase diagram of the MnAs—MnSb system [39, 40]. 
The shaded area that indicates the MnP-phase disappears for x ~ 0.1 of MnAs;_,Sb,. 
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Fig. 3.20. A phase diagram of MnAsj_,Sb,. The shaded area denotes the MnP-type structure. 
The compound takes the NiAs-type structure in the other area. Ferromagnetism realizes for all 
composition. For x < 0.1, a first order transition occurs at Tc 


For x > 0.1, magnetic phase transition at 7c is of the second order. Jc takes the 
minimum value for x = 0.35. It is also interesting that Jc and @p are significantly 
different in composition region 0.1 <x < 0.7. Magnetic moment per Mn shows weak 
dependence on x [41] but almost constant in the range 0.15 < x < 0.7, where the 
minimum of 7c (x = 0.35) occurs. We will discuss the minimum of Tc later. 

As shown in Fig. 3.21, the a-axis increases drastically with Sb composition. The 
difference of the a-axis between MnAs and MnSb is 13%, while variation of the 
c-axis is much smaller (about 1% between MnAs and MnSb). The small jumps of 
lattice parameters at the compositions designated by the vertical dotted lines are ex- 
plained by the composition dependence of Jc shown in Fig. 3.20. If Tc is higher 
than the room temperature, the lattice parameters have an effect of spontaneous 
magnetostriction as shown in Fig. 3.22. On the other hand, when 7c is lower than 
room temperature (0.06 < x < 0.08), there is no effect of the magnetostriction. 
As shown in Fig. 3.22, magnetostriction of MnAso.7Sbo.3 at T = 80K is +1% for 
the a-axis and —0.7% for the c-axis. This corresponds to volume magnetostriction 
AV/V = 41.5 x 1072. 

The experimental result that Jc takes a minimum for the composition (Fig. 3.20) 
is phenomenologically explained from the viewpoint of localized model similar to 
the model of Bean and Rodbell shown in Sect. 3.2.1. 


34 3 Properties of the Compounds with NiAs-Type and Cu2Sb-Type Structures 
4.2 
MnAs,_-,Sb, 


4.0 
at R. T. 


3.8 emi 7 : 


3.6 


a(A) 


5.70 


0.0 0.2 0.4 0.6 0.8 1.0 
x 


Fig. 3.21. Lattice constants of MnAs,_,Sb, at room temperature [39] 
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Fig. 3.22. Temperature dependence of lattice constants of MnAso.7Sbo.3. Spontaneous mag- 
netostriction arises below 227 K (=7c) [22] 


We assume that 7c is a function of exchange parameters J, and J,, which change 
linearly with lattice parameters c and a around co and apo, respectively, as follows: 
Ja(a) = J(ao) + A(a— ao), 
J-(c) =J(co) + B(c— co), (3.13) 
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where A and B are constants. Therefore, change of 7c due to lattice change is 
written as 


dTo(Je,Ja) = (8Te/8Jq)(OJa/Aa) da + (ATe/dJ-)(AJ-/Ac) de (3.14) 
and 


Tc(a,c) — Te(ao,co) = (OTe /0Ja)A(a — ag) + (OTe /OJ-)B(c — co) 
= a(a—ay) + B(c—co), (3.15) 


where (0Tc/OJq)A and (0Tc¢/0J-)B are replaced by o and f, respectively. We insert 
the experimental lattice constants and Tc of MnSb for Tc (ao, co). By using (3.13) and 
the composition dependence of lattice constants shown in Fig. 3.21, we can calculate 
Tc of MnAs,_,Sb, as plotted in Fig. 3.23. We note that for x > 0.6 (where T¢ is 
higher than room temperature) we extrapolated the lattice constant curve (Fig. 3.21) 
in the region x < 0.5 so as to exclude the effect of ferromagnetism. 

Results shown in Fig. 3.23 were obtained by assuming linear change of exchange 
parameters J, and J, against corresponding lattice constants. In calculation we set 
the coefficients @ = 1360K/A and B = —3000K/A. With these parameters we 
successfully reproduced the minimum of 7c against composition of MnAs,_,Sb, 
qualitatively. It is apparent that anomalous composition dependence of the c-axis 
gives rise to the composition dependence of 7c, though there is no theoretical foun- 
dation for the anomalous composition dependence of the c-axis. We also see from 
Fig. 3.22 that expansion of the a-axis as well as the shrink of c-axis brings about in- 
crease of Tc. This corresponds to the sign of a and B. If elastic constants along 
the a- and c-directions are available, it is possible to evaluate quantitatively the 
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Fig. 3.23. Tc of MnAsj_,Sb, calculated using J, and J; whose temperature change was eval- 
uated from experimental values of the lattice constants. The minimum of Jc is caused by 
anomalous variation of the c-axis for the composition (Fig. 3.21 and also see text) 
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composition dependence of Jc as well as the temperature dependence of sponta- 
neous magnetostriction and magnetic moment on the basis of the Bean—Rodbell the- 
ory [21]. Lattice constants of MnSb and CrSb under high pressure up to 1.8 GPa 
have been measured by Nagasaki [42, 43]; however, a more precise measurement 
using state-of-the-art equipment is expected. 

Next we describe metamagnetic behavior of MnAs;_,Sb,. The metamagnetism 
of MnAs just above 7c arises from the fact that ferromagnetic—paramagnetic transi- 
tion at Tc is of the first order (7) > 1). In MnAs _,Sb,, similar metamagnetic tran- 
sition induced by applied magnetic field occurs only in the composition region of 
x < 0.08 (see Fig. 3.20). On the other hand, for x > 0.1, metamagnetic transition 
does not occur; however, metamagnetic-like behavior remains near x = 0.1, which 
means the parameter 7 in the argument in Sect. 3.2.1 is slightly smaller than 1. As 
an example, metamagnetic-like behavior of MnAso.7Sbo,3 just above 7c is shown in 
Fig. 3.24 [22]. 

This field-induced metamagnetic-like transition does not accompany the struc- 
tural transformation between the MnP-type and NiAs-type structures. As shown in 
Appendix, metamagnetic magnetization process is generally given by the expression 


H =A\o +A20°+A30°+..., (3.16) 


where H is a magnetic field, o magnetization, and A, etc. are expansion coefficients. 
The form of (3.16) is common to the Bean—Rodbell model and the itinerant electron 
model developed by Yamada [44], though physical meaning of the coefficients is 
different. The experimental results shown in Fig. 3.24 are well reproduced by (3.16) 


M (emv/g) 


Fig. 3.24. Magnetization curves of MnAso.7Sbo.3. The curves are measured at T = 4.2, 77, 
180, 200, 210, 220, 225, 230, 240, 250, 255, 260, 270K from larger magnetization. Inset: 
magnetic fields at inflection points are plotted as a function of temperature. From this figure 
we see that Jc = 230 K (Reproduced from [22]) 
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with temperature-dependent A;, Az, and A3. However, both the above two models 
cannot explain successfully the temperature change of Az [22]. We also note that 
there is some allowance for the experimental values of A;, Az, and A3. 


3.2.4 Effect of High Pressure on MnAs,_,Sb, 


Pressure effects on Tc of MnSb and MnAspg 8gSbo.12 [45] and changes of magnetic 
states of MnAso 7Sbo.3 [46] have been investigated under high pressure. Tc of MnSb 
decreases by 200K (from 590 to 390K) by applying pressure 5 GPa [47]. Under 
this pressure, the a- and c-axes decrease about 2% (from 4.13 to 4.03 A) and 4% 
(from 5.78 to 5.55 A), respectively [43]. This corresponds to the compressibility 
K = —(1/V)(dV/dP) = 8.6 x 107!! Pa! (= 8.6 x 107!cm* dyn7!), which is an 
ordinary value for this kind of intermetallic compounds [42]. The elastic constant of 
MnSb is anisotropic and the compressibility is much larger than that of a material 
which has closed packed structure (e.g., four times larger than that of Cu and eight 
times larger than that of Ni). Unlike MnAso.7Sbo.3 shown in Fig. 3.22, the temper- 
ature variation curve of c-axis of MnSb shows no anomaly at 7c while the a-axis 
shows small positive magnetostriction [48]. This fact indicates that the ferromag- 
netism is stabilized by the increase of a-axis. If we are based on the localized model, 
the positive magnetostriction means that increase of a-axis results in increase of an 
exchange interaction parameter J, connecting Mn spins within the c-plane. This is 
consistent to the behavior of MnAsj_,Sb, with x > 0.6 shown in Figs. 3.20 and 3.21, 
where Jc increases as a-axis increases (c-axis is almost constant) with increasing 
composition x. Therefore, the decrease of Tc of MnSb by application of pressure 
is mainly caused by the decrease of a-axis. It is apparent that there is a close cor- 
relation between magnetism and lattice parameters. To clarify the mechanism from 
microscopic electronic structure is a problem left for the future. 

Next, we discuss pressure effect on Tc of MnAsg.ggSbo.12. Unlike MnAs, this 
compound shows a second order transition at 7c. As shown in Fig.3.25, T¢ de- 
creases and temperature interval AT of hysteresis for magnetization increases with 
increasing pressure. The curve for 1.2 GPa suggests a structural change caused by 
the pressure from the NiAs- to the MnP-types. 

MnAso.ggSbo,12 has an Sb composition x just above the critical composition x = 
0.1, as seen in Fig. 3.20; therefore, shows at Jc the second order transition having 
the characteristic near the first order transition. The hysteresis generated by applied 
pressure shown in Fig. 3.25 means that pressure changes the second order transition 
to the first order transition at Tc. (The large hysteresis shown in Fig. 3.25 indicates the 
first order transition though the order is unclear especially under pressures P = 0.8 
and 1.0 GPa). 

We will show that these experimental results are explained by the Bean—Rodbell 
model. We assume 71 = | (see (3.11) and (3.12)) as this compound locates near 
the border between the first order and the second order transition. Then, by using 
parameters J) = 300K and KB = 5.3 x 10~!! (CGS), which are the same as those 
as MnAs, we can calculate by making use of (3.12) the magnetization curves plotted 
in Fig. 3.26. 
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Fig. 3.25. Temperature dependence of magnetization (in arbitrary scale) of MnAso.ggSbo,12 
under high pressure. The labels 1, 2, 3, 4, 5, 6, and 7 denote the pressure 0, 0.2, 0.4, 0.6, 0.8, 
1.0, and 1.2 GPa, respectively (Reproduced from [45]) 
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Fig. 3.26. Magnetization (o) vs. temperature curves of MnAsg.ggSbo.12 (1 = 1, To = 300 K, 
BK =5.3 x 10~!! (CGS)) under atmospheric and high pressure (P = 1 GPa). 15 K width of 
hysteresis is seen for P= 1 GPa 


Calculated results shown in Fig. 3.26 exhibits that pressure alters order of tran- 
sition from the second to the first, and the calculated values of A < 15 K is in fairly 
good agreement with the experimental results of A ~ 19K (Fig. 3.25), even with- 
out parameter adjustment. It is noted that the molecular field theory based on the 
localized picture (i.e., the Bean—Rodbell model) is successful in explaining various 
properties of itinerant electron system such as present compounds. 
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Fig. 3.27. Magnetization—temperature curves of MnAso.7Sbo,3 under atmospheric and high 
pressure (1.2 GPa) (Reproduced from [46]) 


Next, we show pressure effect on MnAso.7Sbo.3 in magnetic field up to 7T 
and in quasi-hydrostatic pressure up to 1.2 GPa [46]. We have already noted that 
this material exhibits a second order transition at Tc; however, this compound shows 
metamagnetic-like behavior just above Tc as shown in Fig. 3.24. Figure 3.27 shows 
magnetization—temperature curves under pressures 0 and 1.2 GPa. We see that Tc de- 
creases by 70 K due to 1.2 GPa pressure. (7c of MnAso.ggSbo.12 decreases by 105 K). 
It is not clear in Fig. 3.27 whether the order of transition at Tc changes by pressure 
from the second to the first. However, temperature dependence of magnetization un- 
der 8 and 10GPa shown in Fig. 3.25 suggests that phase transition under 1.2 GPa 
shown in Fig. 3.27 is of the first order. It is possible to clarify whether the transition 
is of the first order or of the second by measuring the hysteresis A or temperature 
dependence of spontaneous magnetization under pressure. 

Finally, we show pressure dependence of magnetization curves at T = 4.2K 
in Fig. 3.28. It is clearly seen that spontaneous magnetization at 4.2 K decreases 
largely by pressure; the spontaneous moment per Mn atom M(pUp) is found to be ex- 
pressed by M = 3.185 —0.126P—0.052P”, where P means the pressure in the unit of 
GPa. A relation between spontaneous magnetic moment M, and Jc under pressure, 
dinM, /dP = (2/3) d1n7¢/dP, has been obtained for a weak itinerant ferromagnet by 
Takahashi [49]; however, this relation seems to be inapplicable to a strong ferromag- 
net MnAso,7Sbo.3. Pressure dependence of spontaneous magnetization at T = 4.2K 
is explained qualitatively by the energy band calculation for changing lattice param- 
eters [46]. Note that there are many problems left for the future. 


3.2.5 Mny,_,Cr,As, Mn,_,Ti,As, etc. 


As MnAs is unstable against the MnP-type deformation, a mixed compound 
Mn,_,Cr,As (CrAs has the MnP-type structure) has the MnP-type structure in 
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Fig. 3.28. Pressure effect on the magnetization curve of MnAso,7Sbo.3 at T = 4.2 K (Repro- 
duced from [46]) 


low temperature region even when x is small, and undergoes the structural trans- 
formation to the NiAs-type when temperature rises. These structural properties of 
Mn,_;Cr,As are similar to those of MnAsj_,P, [50]; however, magnetic properties 
are somewhat different. Existence of two kind of 3d atoms in Mn,_,Cr,As makes 
the magnetic properties a little complex. 

By substituting a nonmagnetic atom Ti for Mn of MnAs and MnSb, density 
of magnetic atom in the compound decreases and the NiAs-type structure is sta- 
bilized, except the vicinity of TiAs (TiSb). 7c of these mixed compounds shows 
rapid decrease with Ti composition, approaching 0K for x = 1. Mn,_,Ti,Sb be- 
comes spin glass in low temperature for the compositions around x = 0.8 [51]. Sub- 
tracting Pauli paramagnetic contribution (¥p = 1 x 10-% emu g~!), susceptibility of 
Mn,_;Ti,As obeys the Curie-Weiss law [14]. Form this Curie—Weiss type suscepti- 
bility, the paramagnetic Curie temperature ©, and the effective Bohr magneton per 
Mn atom, Ler = gilg./S(S+ 1), were evaluated. By setting the g-factor 2 in Herr, we 
can calculate gSupg (S: spin) = Pa (atomic moment). Figure 3.29 illustrates Op, Tc, 
and Pa as a function of Ti composition. Mn;_,Ti,Sb also shows the similar behavior 
[14]. It is notable that Pa is almost independent on Ti composition. There has been 
report on temperature dependence of lattice constants and pressure effect on Tc of 
Mn,_,Ti,As (0 < x < 0.9) [21]. An itinerant electron theory is desired to explain the 
properties shown above. 


3.2.6 Magnetic Refrigeration Using MnAs and the Related Compounds 


If an applied magnetic field largely reduces entropy of a material, this material 
can be a candidate for a refrigerant of magnetic refrigeration. Enhancement of 
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Fig. 3.29. Tc, Op, and Pa (= gSip/Mn, see text) of Mnj_,TiyAs. Mnj_,TiySb shows similar 
behavior (Reproduced from [14]) 


magnetization due to applied magnetic field means increase of magnetic order or 
reduction of magnetic entropy. For a ferromagnet, large increase of magnetization 
due to applied field occurs in the temperature region around 7c; just above and be- 
low Tc for the second order transition and only just above 7c for the first order 
transition (cf. Figs. 3.4 and 3.6). A ferromagnet that shows a first order transition 
(e.g., MnAs) shows metamagnetic transition just above 7c. Metamagnetic transition 
field Hc increases with rising temperature (see Fig.3.7). In Fig. 3.30, magnetiza- 
tion curves above and below 7¢ are show schematically for a ferromagnet with the 
second or the first order transitions at Tc. 

If we measure or calculate magnetization curves at various temperatures, it is 
possible to evaluate corresponding magnetic entropy using the expansion formula 
of o (3.3). In an adiabatic condition where total entropy is constant, reduction of 


3 Temperature dependence of decrease of magnetic entropy ASy by application of external 
field H is calculated from difference of specific heats in H = 0 and H = H. ASy is also 
calculated by integrating the Maxwell relation (0S/0H)r = (0M /0OT )x of thermodynam- 
ics as 


H 
Sy ={ (OM/dT)y dH, 
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Fig. 3.30. Schematic magnetization curves above and below Tc 
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Fig. 3.31. Temperature dependence of magnetic entropy (a schematic figure for the case of a 
second order transition). See text for 7; 


magnetic entropy due to applied field enhances entropy of lattice vibration, bringing 
about increase of temperature by AT. This is called magnetocaloric effect. A material 
that has large AT is suitable for magnetic refrigerant. 

We explain briefly how AT arises below [52]. Figure 3.31 shows schematic tem- 
perature dependence of magnetic entropy Sq with and without magnetic field. 

Total entropy S; of a magnetic material is Sy + Sz (SL: lattice entropy). Assum- 
ing that the lattice specific heat Cy, is constant in the range 7; < T (Tj is a temperature 
just below Tc), Sy is calculated as 


T 
S(T) —S(T) =} C.dT/T =C,n(T/T)). (3.17) 


This equation means that the lattice entropy is a monotonically increasing function 
of temperature. Taking account of temperature dependence of Sy in Fig.3.31, we 
can express Sy schematically as shown in Fig. 3.32. 


where M is magnetization under a magnetic field. We need to measure magnetization in 
various magnetic fields to obtain ASy. By integrating 0M/0T with respect to H at respec- 
tive temperatures, we can calculate temperature variation of ASy. We note that M in this 
argument means magnetization in a state in which all the atomic moments align in the 
direction of H, that is, a state without effect of anisotropy and magnetic domain. 
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Fig. 3.32. Schematic curves of total entropy Sr vs. temperature T in H = 0 and H =H 


If we apply H adiabatically (i.e., Sp is unchanged) to a magnetic material at 
temperature T = T, (H = 0), temperature rises up to T = T) by AT. The resulting 
state at T = T, (H = H) is equivalent to a state in which H is isothermally applied 
to a material at T = T (H = 0). In the former adiabatic process, entropy of spins 
decreases (because spins align along the applied field) and lattice entropy increases, 
and then temperature of the material increases. On the other hand, dissipation of heat 
keeps the temperature unchanged in the latter isothermal process. 

The process of magnetic refrigeration is explained next. By applying H, the state 
(T),H =0), which is taken as an initial state here, changes to the state (),H = H). 
In this process, generated heat is released to outside. Adiabatic removal of magnetic 
field brings the state (71,H) to (1;,H =0). In this process, the temperature of the 
material lowers from 7> to 7;, namely the material absorbs heat from outside (i.e., air 
around). Repetition of the similar process is the magnetic refrigeration. The magnetic 
material must show large magnetocaloric effect within a finite temperature range 
because the temperature of refrigerant decreases during the refrigeration process. 
A magnetic material that has large AT in a wide temperature range is favorable. From 
Fig. 3.32 we see that dSy/dT = dS, /dT ~ ASm/AT for T > Tc. Since dS, /dT = 
C./T as seen in (3.17), we have 


AT ~ TASm/CL ~ TcASm /CL. (3.18) 


Therefore, AT is proportional to reduction of magnetic entropy due to applied field. 
For T < Tce — AT, dSp/dT = dS, /dT + dSy/dT = C./T + Cu /T is approximately 
satisfied and then 


AT ~ Tc ASm/(CL+Cu(H = 0)), (3.19) 


where Cu (H = 0) is the magnetic specific heat in H = 0. From (3.18) and (3.19) 
we see that AT for T < Tc is expected to be smaller than AT for T > Tc because of 
the difference between denominators of (3.18) and (3.19). Wada et al. investigated 
MnAs and MnAsj_,Sb, as magnetic refrigerant [53-55]. We show some of their data 
in Fig. 3.33. 
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Fig. 3.33. (a) Magnetocaloric effect AT of MnAs evaluated indirectly from temperature 
variations of total entropy Sp and magnetization. (b) Decrease of magnetic entropy —AS)4 
of MnAs;_,Sb, measured by applying magnetic field 5 T (Reproduced from [53]) 
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Fig. 3.34. Schematic figure to show the mechanism of refrigeration by pressure for 
MnAs;_,Sb, with small x. This figure is for a first order transition at Tc 


In Fig. 3.33b, AS\y at H = 5T was evaluated from temperature change of mag- 
netization in various magnetic fields, whereas in (a), AS\y was determined indirectly 
from the data shown in (b) and temperature dependence of entropy in zero field [55]. 
Method of direct measurement of the magnetocaloric effect and experimental data 
for Gd metal have been reported together with analysis based on the molecular field 
theory in [56]. 

We have discussed in Sect. 3.2.4 that Tc of MnAs;_,Sb, with small Sb content 
decreases remarkably by application of pressure. It may be possible to apply the 
compounds for refrigeration using pressure instead of magnetic field for the refrig- 
eration process. As seen from Fig. 3.34b, removal of pressure isothermally reduces 
magnetic entropy. This is equivalent to application of magnetic field isothermally. 
It is interesting to clarify which is more practical. 
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3.3 MnSb and MnBi 


We have already shown basic properties of these ferromagnetic compounds in 
Sects. 1.4 and 2.1. As MnSb has already been mentioned in Sect. 3.2.4, we explain 
MnBi in this section. As shown in Table 1.2 and Fig. 3.11, MnBi takes the NiAs- 
type structure with lattice constants larger than MnSb, and magnetic moment per Mn 
atom is 3.8Ug. MnBi shows a first order transition with temperature hysteresis 15 K 
at Tc = 628 K [57]. Figure 3.35 shows temperature dependence of spontaneous mag- 
netic moment, where data not only for low temperature phase (LTP) MnBi but also 
for Mn, .9gBi (quenched high temperature phase: QHTP) are plotted. The LTP shows 
a first order transition at Tc, which is accompanied with crystallographic transition to 
the nonstoichiometric Mn, ogBi. The QHTP has smaller moment and lower 7c than 
those of the LTP MnBi. These properties are similar to those of Mn;+,Sb; however, 
x is limited to 0.08 for Mn,,,Bi. 

The QHTP MnBi separates into Mn and Mn-Bi liquid at the peritectic tempera- 
ture of 446°C. As MnBi has large magnetic anisotropy with magnetic moment par- 
allel to the c-axis (anisotropy constant K ~ 2 x 107 ergem~* = 2 x 10°J m4) [57], 
there were attempts for applications to permanent magnet [58]. However, this ma- 
terial is not suitable for a magnet because of the relatively small magnetization and 
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Fig. 3.35. Spontaneous magnetization 477M and inverse susceptibility y—! of MnBi (T < Tc) 
and Mn; 0gBi (high temperature phase). Old data are also shown by the dashed curve for 
comparison 
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difficulty in preparation of compound. This material may be utilized for other type of 
application such as magnetic memory. There has been no theoretical study to explain 
the peculiar properties of MnBi. 


3.4 CrAs and Related Compounds 


3.4.1 Anomalous Magnetic Transition of CrAs 


In contrast with the ferromagnetic Mn-pnictides with MnX-type, the compounds 
with CrX (X = As or Sb) type take generally antiferromagnetic ordering. As mag- 
netic properties of CrAs are especially interesting, we describe mainly CrAs and its 
related compounds. 

As explained in Sect. 2.1, CrAs takes the double helical magnetic ordering below 
Ty = 272 K [59], and undergoes a first order transition at Tj to paramagnetic state in 
high temperature region above Ty. We show susceptibility of CrAs and CrAso 9Po,; in 
Fig. 3.36 (Ido, unpublished). Similar data for CrAs;_,P, have already been reported 
by Selte [60]. 

As seen in Fig. 3.36, susceptibility of CrAs increases with increasing temperature 
without anomaly at Ty, suggesting strong itinerant character. It is known that CrAs 
with the MnP-type crystal structure takes the NiAs-type structure in the temperature 
region higher than 7, = 1100K, and the magnetic susceptibility obeys the Curie— 
Weiss law in this temperature range, which is not included in Fig. 3.36 [50]. For 
CrAs}_xSb, (x < 0.6), the Curie-Weiss law in the region of T > 7; is clearly observed 
as T, decreases due to substitution of Sb for As [61]. For example, CrAso.9Sbo,; has 
transition temperature 7, = 780K and shows the Curie-Weiss behavior of suscep- 
tibility in T > 780K. From the Curie—Weiss law, moment of Cr in CrAsg.9Sbo.1 
is evaluated from the effective Bohr magneton to be 2.8p, which is equal to the 
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Fig. 3.36. Susceptibility of CrAs and CrAso.9Po.;. The magnetic order occurs in CrAs below 
Tx, while no magnetic order occurs in CrAso.9Po.1 down to 0K 
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Fig. 3.37. Temperature change of integrated intensity of (002) satellite of CrAs (Reproduced 
from [59]) 


moment of Cr in CrSb determined by neutron diffraction measurements (see also 
Sect. 3.5). It is interesting that susceptibility of CrAs in the NiAs-type phase obeys 
approximately the Curie-Weiss law expected for Cr** (i.e., 3g). This phenomenon 
will be shown and discussed later in Fig. 3.42. Ty of CrAs was determined by neu- 
tron diffraction measurements (Fig. 3.37) [59, 62]. Magnetic phase transition at Ty 
is of the first order accompanied with discontinuous change of lattice constants. We 
will discuss in detail later the magnetovolume effect of CrAs. 

It has been found from specific heat measurements that latent heat accompanied 
with magnetic phase transition is 5.3 x 10? J mol~! (=0.98 cal g~!) [63,64]. Temper- 
ature dependence of electric resistivity is similar to that of a gap-type antiferromagnet 
such as Cr and Cr;_,Mn,, (x < 0.25), whose resistivity increases first with decreas- 
ing temperature below Ty and decreases next on further cooling [63]. It is difficult to 
measure resistivity of CrAs accurately because of cracks arising due to large volume 
change at 7y. This is similar to the situation in MnAs. 

In Fig. 3.38, magnetic moment per Cr atom for some Cr pnictides is plotted as a 
function of length of the a-axis. This figure is similar to Fig. 3.11, where Mn-moment 
is plotted against a-axis length for the Mn-pnictides. We note again that the a-axis 
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Fig. 3.38. Relation between moment per a Cr atom and the a-axis at 4.2 K. Filled and open 
symbols denote the MnP-type and NiAs-type structures, respectively. This figure has been 
drawn by making use of the data in [65] and [66, 67]. 


for the orthorhombic MnP-type structure in Fig. 3.38 corresponds to the a-axis of the 
hexagonal NiAs-type structure. The reason why we take only a-axis as a parameter 
for (Cr) in Fig. 3.38 is that the c-axis variation from 5.45 A of CrSb to 5.36A of 
CrP is negligibly small compared with the large a-axis variation from 4.13 A of CrSb 
to 3.11 A of CrP. 

Figure 3.38 shows apparently that magnetism of CrAs is very sensitive to the 
change of a-axis, which will be discussed in detail in the next section. It may need to 
point out that there is little change of u(Cr) around a = 3.75 A where the structural 
change occurs. This may mean that the distortion to the MnP-type structure does not 
affect so much. This is sharply different from MnAso.9Po.1, whose moment exhibits 
large change accompanied with the structural change as shown in Fig. 3.14. 


3.4.2 Critical Lattice Constant of CrAs;_,P, and Cr;_,M,As 
(M = Mn, Ni, etc.) 


We see the anomalous magnetic properties of CrAs more clearly in the behavior of 
mixed compounds. A typical example is a CrAs-based compound in which small 
fraction of As is replaced by P. Both electron number and crystal structure remain 
unchanged by the P substitution for As. According to [60], the double helical or- 
dering with the first order transition at Ty ~ 250K disappears for CrAs;_,P, with 
x > 0.03 [60], which is also roughly expected from Fig. 3.38. This is a very interest- 
ing phenomenon whose mechanism has not been clarified yet. 

We show and discuss experimental data of magnetovolume effect of CrAsj_,P, 
in the following part of this section. Figure 3.39 shows temperature change of lattice 
constants of some of the mixed compounds derived from CrAs [68]. As we noted 
in Sect. 1.3, the A-, B-, and C-axes here correspond to the c-, a-, and b-axes of the 
hexagonal NiAs-type structure, respectively. 
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Fig. 3.39. Temperature change of lattice constants of CrAs and related compounds. Note that 
the B- and C-axes correspond to the a and b-axes of the hexagonal NiAs-type structure, re- 
spectively (Reproduced from [68]) 


In Fig. 3.39 data for the A(c)-axis, where (c) means the A-axis corresponds to the 
c-axis, are not plotted because the discontinuous change of A(c)-axis at Ty is much 
smaller than those of B(a)- and C(b)-axes. For example, in the case of CrAso 9sPo.02, 
the discontinuous change of the B- and C-axes are 5.5% and 1.1%, respectively, 
while that of the A(c)-axis is only 0.4%. Therefore, as generally seen in Fig. 3.39, 
magnetovolume effect takes place mainly in the B(a)-direction. It is noted that the 
change of the B-axis is exceptionally large. We may regard the change in the C(b)- 
direction at Ty, which has an opposite sign to that of B(a)-axis, is only induced by 
the large change of the B(a)-axis. To obtain the normal thermal expansion curve 
of the B(a)-axis, temperature dependence of the B(a)-axis of CrAso.9Po.1, which 
shows no magnetic ordering at T = 0K, was measured for the temperature from 80 
to 330 K. As shown by the dotted curves in Fig. 3.40, the part of the curve below 80 K 
was extrapolated on the basis of the specific heat (H. Ido and T. Kamimura, unpub- 
lished) and the Griineisen relation [69], which says that specific heat is proportional 
to a thermal expansion coefficient. By applying the normal thermal expansion curve 
presumed for CrAso9Po.1, magnetovolume effect at T = 0K has been estimated 
for CrAs and CrAso 9gPo.92 to be 4.6% (0.16 A) and 5.6% (0.19 A), respectively, 
as shown in Fig.3.40. These values are extremely large. As no magnetic ordering 
occurs in CrAso.9Po.1, this compound shows no magnetovolume effect at T = 0K. 
We see from Fig. 3.39 that magnetic ordering of CrAs and CrAso 9gPo.02 is stabilized 
by the large lattice distortion in the B(a)-direction. We can also say that increase of 
elastic energy due to magnetostriction is larger than decrease of magnetic energy in 
CrAso.9Po.1; therefore, this compound has no magnetic ordering. Therefore, in CrAs 
and the related compounds there seems to be a close relation between the occurrence 
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Fig. 3.40. Spontaneous magnetostriction of CrAs;—,P, and critical B-axis Bc (Reproduced 
from [68]) 


of magnetic order and the length of B(a)-axis. As indicated in Fig. 3.40, the criti- 
cal length of B(a)-axis, Bc, is about 3.37 A, namely magnetic ordering occurs in the 
compounds with B larger than Bc, and no magnetic ordering occurs if B < Bc. CrAs 
with B larger than Bc has moment 1.7Ug and Ty = 250K. Figure 3.40 indicates that 
magnetic ordering of CrAs is expected to disappear by more than 5.6% compression 
of the B(a)-axis. Figure 3.40 also suggests that, if spontaneous magnetostriction were 
fixed at the value just below the first order transition point Ty, the Néel temperature 
would be more than 350K. 

As seen from Fig.3.39, Ty (a first order transition) increases due to Mn sub- 
stitution for Cr of CrAs. In Fig.3.41a, Ty determined by thermal expansion curve 
is plotted against composition of various substitution elements. We see that Ty de- 
creases (increases) when the B(a)-axis shrinks (increases) due to the substitution. 
On the other hand, in Fig. 3.41b, Ty is plotted not for composition but for the length 
of B(a)-axis just above Ty (=250) K, which is almost close to the B-axis length at 
T = OK in the absence of magnetostriction. 

As seen from Fig. 3.40, the important point we see from Fig. 3.41b is that there 
is a critical length of the B-axis Bc = 3.38 A, which is independent of the kind of 
substituting atom. To construct a magnetic phase diagram, the crystallographic tran- 
sition temperatures T, between the MnP-type and the NiAs-type, which are taken 
from Fig. 3.42, are also shown in Fig. 3.41b. 

In Fig. 3.41b, the labels “Pauli para” and “Curie—Weiss type para” indicate the re- 
gions that magnetic susceptibility behaves as those in Fig. 3.36 and as those above 7 
in Fig. 3.42, respectively [61]. In Fig. 3.42, we show moment of a Cr atom [Uc; (evalu- 
ated from the Curie—Weiss law of susceptibility), Op, etc. As CrAsg.9Sbo.; has higher 
Ty and lower 7, compared with CrAs (see Fig. 3.36 and [50]), the Curie—Weiss-type 
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Fig. 3.41. (a) Composition dependence of Ty when Cr or As of CrAs is replaced by other 
elements, which are denoted by the following symbols: filled triangle (Ti); filled square (Mn); 
open triangle (Fe); open circle (Co); filled circle (Ni); open square (P); open inverted triangle, 
filled inverted triangle (Sb). (b) Ty of these compounds plotted as a function of length of the 
B-axis just above Ty. Data in Fig. 3.42 has also been used to make this figure. 7; is transition 
temperature of the structural transformation. See also the text (Reproduced from ref. [68]) 
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Fig. 3.42. Susceptibility of CrAs;_ Sb, (x < 0.6). Ty and T/ are the Néel temperature and tem- 
perature of the first order structural transformation (see also Fig. 3.43) (Reproduced from [61 ]) 
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behavior of susceptibility appears more clearly especially for the compound with 
x = 0.3. We will discuss again the results shown in Fig. 3.42 in Sect. 3.4.4. 


3.4.3 Phenomenological Theory for the First Order Transition of CrAs 


The phenomenological theory for the first order transition of MnAs (see Sect. 3.2.1 
and Appendix) is applicable to CrAs. We note that physical meaning of the terms of 
free energy expansion is different. When magnetization of space variation is written 
as M(r) = M, exp(igr) (q = 2/2.67b along the b-axis in the case of CrAs. See 
Fig. 2.1), we may assume molecular field AM, within the molecular field theory. 
Then we have magnetic energy Em = —(1/ 2)AM?. By setting an assumption that 
the molecular field coefficient A depends linearly on crystal volume as A = Ao[1 + 
B(V —Vo)/Vo)], we have an expression 


Em = —(1/2)Ao[1 + B(V — Vo)/Vo)] 0”, (3.20) 


where we wrote M, = o, which corresponds to the relative magnetization in the Bean 
and Rodbell theory for MnAs (see Sect. 3.2.1). As mentioned in Sect. 3.2.1, elastic 
energy due to lattice distortion is written as 


Ex = (1/2K)[(V — Vo)/(Vo)}? (3.21) 
and magnetic entropy per unit mass of the crystal is given by 
S(o) = Nkp[S(0) + ao? + bo* +co® +... (3.22) 
By setting derivative of E = E,, + Eq with respect to V zero, we have a relation 
(V —Vo)/Vo = (KBAo/2)o?. (3.23) 


Inserting (3.23) into E (= Em + Eq), we have an expression of free energy G, which 
is minimized with respect to V as 


G/NkgTy = —tS(0) — a(t— 1)0? — b(t—n)o* — cto® +..., (3.24) 


where 7p, t, and 7) are expressed by the following equations (3.25) and (3.26), re- 
spectively, 


—(Ao/2aNkg)=Ty and T/Ty=t (3.25) 
and 
n = —(a*/2Nkgb)KB7T (> 0). (3.26) 


We see that (3.26) is the same as (A8) in Appendix; therefore, a first order tran- 
sition occurs when 7) > 1, namely if the parameter expressing volume dependence 
of magnetic energy of CrAs is large enough for 7 > | to hold, the magnetization 
o (= M,) disappears discontinuously at Ty. As (3.23) shows that the spontaneous 
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magnetization is proportional to B, the experimental results shown in Fig. 3.39 
clearly show that magnetic energy is very sensitive to crystal volume (mainly to 
the B(a)-axis). Results of neutron diffraction measurements revealed that o of the 
double helical CrAs at Ty is about 0.7, which is close to the value of MnAs. Spon- 
taneous magnetostriction (V — Vo)/Vo at Tx is 0.022 for CrAs [68] and 0.021 for 
MnAs. If we express (3.23) as (V — Vo) /Vo = AO”, the coefficient A of CrAs is very 
close to A of MnAs; however, there is contribution of structural transformation to the 
magnetostriction at Tc in the case of MnAs. Therefore, we may regard the value of A 
for CrAs as twice as that of MnAs. The magnetism accompanied with the first order 
transition at Ty of CrAs, which is closely correlated to the critical length Bc of the 
B-axis, is a very interesting problem, and is expected to be clarified in the future in 
the view point of itinerant electron picture. 


3.4.4 CrAs;_,Sb, 


As noted in Sect. 2.1, CrSb takes the NiAs-type structure. In contrast to CrAs;_ Py, 
substitution of Sb for As in CrAs reduces 7; (i.e., transition temperature between 
the MnP-type and the NiAs-type structure). See also Fig. 3.42. Figure 3.43 shows a 
magnetic and crystallographic phase diagram of the CrAs—CrSb system constructed 
from temperature dependence of lattice constants obtained by X-ray diffraction mea- 
surements [61]. 

The rise of Jy in the region near CrAs seen in Fig. 3.43 is probably arising from 
the competition between magnetic energy and elastic energy due to magnetostric- 
tion; decrease of magnetic energy and increase of elastic energy is caused by the 
spontaneous magnetization in CrAs; however, increase of lattice constants due to the 
Sb-substitution reduces the increase of the elastic energy, resulting in stabilization 
of magnetic ordering and higher Ty. By this mechanism the experimental results in 
the region near CrAs is explained. A first order transition occurs when there is com- 
petition between the increase of elastic energy and the decrease of magnetic energy. 
Therefore, if there is no such competition, order of transition becomes the second, 
which is realized in the Sb-composition x between 0.1 and 0.4 as seen in Fig. 3.43. 
Tn seems to take a minimum around x ~ 0.5. This behavior is similar to that in the 
MnSb—MnAs system; however, may have more complex origin. 

Next, we explain 7 for x = 0.6 as an example in Fig. 3.43. On heating the 
lattice constants shows at T/, a large jump as Ac = 0.15A (Ac/c = 0.026), and 
Aa = —0.07 A (Aa/a = —0.018). We note that we defined 7/ as a middle point of 
temperature region about 100K wide where a high-temperature phase and a low- 
temperature phase coexist. Corresponding to the jump of the lattice constants at T/, 
susceptibility shows steep increase at 7,’ as shown in Fig. 3.42. The first order transi- 
tion at 7,’ is quite interesting but left unsolved. 

We pay our attention again to the temperature dependence of susceptibility of 
CrAso,9Sbo,1 shown in Fig. 3.42 in Sect. 3.4.2. Behavior of susceptibility in the tem- 
perature region below 7; is very similar to that of CrAs (see Fig. 3.36). On the other 
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Fig. 3.43. A magnetic and crystallographic phase diagram of CrAs—CrSb system. Thick lines 
indicate that the first order transition takes place. T/ corresponds to a transition that accom- 
panies discrete change of the c- and a-axes; however, detail of this transition has not been 
revealed yet (Reproduced from [61]) 


hand, susceptibility for T > 7; obeys the Curie-Weiss law with Cr moment 2.8 Up 
(estimated using g = 2) and paramagnetic Curie temperature 9p = —77K. These 
values agree well with the values of CrSb [70]. It is noted that CrAsg.9Sbo.; in the 
temperature region of T < 7;, where the crystal structure is with MnP-type, is similar 
to CrAs, while it is similar to CrSb in T > 7, where the crystal structure is of the 
NiAs-type. 

Finally, in Fig. 3.44, we show magnetic moment per Cr atom in the CrAs—CrSb 
system determined by Kallel et al. from their neutron diffraction measurements at 
T = 4.2 K [65]. The discrepancy around x = 0.6 is caused by the differences in mag- 
netic order and crystal structure; on the CrSb side, the compounds have an anti- 
ferromagnetic order similar to that of CrSb with the NiAs-type structure, while on 
the CrAs side the compounds take the double helical ordering similar to that of CrAs 
with the MnP-type structure. However, change of magnetic moment at around x = 0.6 
is not large, which is in contrast to the case of MnAs,_,P,., where large difference of 
moment between the NiAs-phase and the MnP-phase occurs (see Fig. 3.14). 
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Fig. 3.44. Magnetic moment of a Cr atom in CrAs;_,Sb, obtained by neutron diffraction 
measurements. For x < 0.5, the double helical ordering and the MnP-type structure realize 
(i.e., similar to CrAs). For x > 0.5, antiferromagnetism and the NiAs-type structure realize 
(similar to CrSb) [65] 


3.5 CrSb 


CrSb, which takes the NiAs-type structure, is an antiferromagnet with Ty = 718K 
and a simple spin arrangement shown in Fig. 3.45. 

Magnetic moment of a Cr atom is equal or close to Cr>+ moment (3.0ug) with 
g = 2. As shown in Fig. 3.38, the a-axis of CrSb is much longer than that of CrAs, 
while difference of the c-axis between CrSb and CrAs is small (5.65 A for CrSb and 
5.45 A for CrAs). Such a large a-axis of CrSb brings about narrow width of the 3d 
band and enhancement of intra-atomic exchange interaction, consequently the large 
moment. As we mentioned in Sect. 3.4.2, temperature dependence of susceptibility 
of CrSb behaves like an antiferromagnet with two sublattice [71]. This behavior is 
quite different from that in the temperature region Ty < T < J, for CrAs. We can 
obtain Op = —200K and Cr-moment 2.9uUp for CrSb from the Curie—Weiss type 
susceptibility together with assumption of g = 2. The value of 2.9.p thus obtained 
is in good agreement with 3.0 and 2.7uUg by the neutron diffraction measurements 
shown in Fig. 3.45. As compared with CrAs, CrSb seems to have a strong localized 
character. This kind of tendency is also seen in Mn-pnictides such as MnBi, MnSb, 
MnAs, and MnP (see Fig. 3.11). The relation between a-axis and magnetic moment 
shown in Figs. 3.38 and 3.11 seems to suggest a way of approach to understand the 
magnetism of 3d-pnictides. There is relevant discussion in Sect. 4.2 of Part IL. 

Resistivity of CrSb above Ty shows temperature dependence like a semiconduc- 
tor [72]; however, it is unsure whether such behavior of resistivity arises from an 
intrinsic property of CrSb or slight crystal decomposition at high temperatures such 
as Sb-precipitation as seen in MnSb (Sect. 1.4 and Fig. 1.6). 
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Fig. 3.45. Antiferromagnetic spin ordering of CrSb. Moment of a Cr atom has been reported 
as 3.0Up [65] or 2.7 [70] 


3.6 CrP 


CrP takes the MnP-type structure. As shown in Fig. 3.38, the B(a)-axis is shorter by 
0.5A than that of CrAs, which broadens the 3d band. Therefore, density of states 
at the Fermi level is small and no magnetic ordering appears. Susceptibility of CrP 
has minimum at about 200K [73]. Magnitude and temperature dependence of the 
susceptibility are very similar to those of CrAso.9Pp9.; as shown in Fig. 3.36. 


3.7 Properties of Cu,Sb-Type Compounds 


We have already shown magnetic structures of these compounds in Fig.2.2 in 
Sect. 2.1. Band structures are shown in Part II of this book. In this section, we mainly 
discuss the very interesting magnetic properties of Mnj_,Cr,Sb. 


3.7.1 Antiferromagnetic_Ferrimagnetic Transition of Mn2_,Cr,Sb 
and Kittel’s Model 


It has been well known that compounds Mn2Sb whose Mn or Sb is slightly replaced 
by V, Cr, Co, Cu, etc. or Ge, As etc. undergo with increase in temperature a first order 
transition from antiferromagnetic state in low temperature side to ferrimagnetic state 
in high temperature side. Among these compounds, Mnz_,Cr,Sb shows the most 
typical transition with very sharp change of magnetization. Figure 3.46 shows tem- 
perature dependence of spontaneous magnetization measured for single crystalline 
samples, and Fig. 3.47 shows the lattice constants as a function of temperature [74]. 
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Fig. 3.46. Temperature changes of spontaneous magnetization of Mnz_,Cr,Sb. MnoSb is 
ferrimagnetic. Other samples with x except for x = 0.013 show the antiferromagnetic— 


ferrimagnetic transition (Reproduced from [74]) 
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Fig. 3.47. Temperature change of lattice constants of Mnz_,,Cr,Sb (Reproduced from [74]) 


In Fig. 3.46, we see the first order antiferromagnetic—ferrimagnetic transition for 
the samples with x in the region 0.05 < x < 0.16, whereas we also see an interme- 
diate magnetic ordering with a small magnetization for the samples with x < 0.03. 
The small magnetizations of about 3emug™! is generally observed for all samples 
in the antiferromagnetic region due to small amount of impurity of ferromagnetic 
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Fig. 3.48. Schematic picture of the antiferromagnetic—ferrimagnetic transition of 
Mnz_,Cr,Sb. Mn atoms occupy two different sites Mn(I) and Mn(II). Moments at the Mn(I)- 
site and Mn(II)-site are antiparallel. M@, and M) are summation of moments in three adjacent 
layers. When M, and M2 are parallel (antiparallel), we refer it as ferrimagnetic (antiferromag- 
netic) state 


MnSb phase. As shown in Fig. 3.47, the c-axis shows a large jump at the first order 
while the small jump of the a-axis may be attributed to a kind of reaction to the large 
jump of the c-axis. 

Next we explain the change of spin alignment between the antiferromagnetic 
and ferrimagnetic states. The antiferromagnetic—ferrimagnetic transition can be ex- 
plained by the exchange inversion model developed by Kittel [75], where the sign 
of exchange parameter J, shown in Fig. 3.48 changes its sign at a critical length of 
the c-axis. The Kittel model is a phenomenological model, and the stability of the 
ferrimagnetic and the antiferromagnetic states has also been discussed from the view- 
point of itinerant electrons by Motizuki and co-workers (see Part II). We see from 
Figs. 3.47 and 3.48 that the c-axis length is closely related to J., namely it is easily 
seen that the contraction of the c-axis brings about the change of J, from positive to 
negative values. This situation is illustrated schematically in Fig. 3.49. 

We explain on the basis of the Kittel model for the antiferromagnetic- 
ferrimagnetic transition referring to Figs. 3.48 and 3.49. Among energies relevant to 
the transition, magnetic energy 


Eex = —AM,-M) (3.27) 


is directly related to the exchange parameter J... In this equation, M; and M2 are sub- 
lattice magnetizations as shown in Fig. 3.48, and A is a coefficient of the molecular 
field given by 


A = (2J-Z)/(Ngus), (3.28) 


where we note that J., as explained in Fig. 3.48, is an average exchange parameter of 
several different exchange parameters J; connecting Mn-spins belonging to M, and 
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Fig. 3.49. (a) Temperature change of the c-axis with hysteresis is schematically shown. The 
experimental data in the vicinity of 7; in the process of rising temperature are shown in 
Fig. 3.47. The dashed line (co) means the normal thermal expansion curve without the sponta- 
neous magnetostriction (c — cg). (b) Je (see Fig. 3.48) plotted schematically against the c-axis 
length. We assume that sign of J, changes at the critical length c, 


Mb, and Z the number of Mn—Mn spin pairs. We consider here only two exchange 
parameters J; and J> indicated by the dashed lines in Fig. 3.48, and we have a re- 
lation J-Z = 4J, + 2J2. For Z = 6 this relation becomes J. = (2/3)J1 + (1/3)Jo. It 
is important that the coefficient A in (3.28) is proportional to J,, that is, the summa- 
tion of J; and Jz in the form mentioned above. (Note that Kittel did not mention A 
in detail.) Considering that sign of J. or A changes at c = c,, we can expand A as 
follows, 


A(c) =A'(cc)(e— Ce), (3.29) 


where A’(c-) means the first derivative. Substituting (3.29) into (3.27) and replacing 
A’ (ce) by B, we have 


Eex = —BM? (c— ce) cos 0, (3.30) 


where 6 is a relative angle between M, and M2, and M(= M, = M2) is magnitude 
of magnetization of each sublattice. 

Next we consider elastic energy induced by the spontaneous magnetostriction 
given by the following equation: 


Ea = (1/2)[eu (e+ e) +3362.) + €12€xx€yy +613 (Cyy€ze + excl xx), 


= (1/2)[2c1,(Aa/a)? + €33(Ac/c)7] + c12(Aa/a)? 
+ 2c13(Aa/a)(Ac/c), (3.31) 


where éx; = eyy = Aa/a = (a—ao)/ag and e,, = Ac/c = (c — cg) /co are shear (see 


Fig. 3.49a), andc;;’sarestiffnessconstants [76]. By setting —(Aa/a)/(Ac/c) = R(> 0), 
the elastic energy is written as 


Ee = [€33/2 + (c11 +€12)R? — 2Re13](Ac/c)? = K(c—¢0)”, (3.32) 
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with K defined by 
K= [33/2 + (c11 + €12)R? — 2Re13] /c}. (3.33) 


From Fig. 3.47 together with Fig. 3.49a, we see that Aa/a is smaller than Ac/c, so 
that R ~ 0.3. Thus, neglecting small terms in (3.33), we find K is mainly contributed 
by the term Kc33/ oc. Therefore, it is seen from (3.32) that the elastic energy in the 
present case is approximately proportional to square of magnetostriction along the 
c-axis. We note that different origin c = c, has been used for Ex in (3.30). Magnetic 
entropy is unchanged during the transition at T,, because this is a transition between 
different magnetic orderings.* We can thus omit the entropy term in free energy. 
As a result, we can write the free energy as 


F = Ex + Ea = —BM?(c— cc) cos @ + K(c — 0)”. (3.34) 


From a condition dF /dc = 0, we have 
c— co = (B/2K)M? cos@. (3.35) 


Inserting (3.35) into (3.34), we have an expression of F which is minimum with 
respect to c as 


F(@) = —B(co — cc)M” cos @ — (B°M*/4K) cos? @ =acos@—bcos’@, (3.36) 


where b > 0 and a is negative (positive) for co > Ce (co < Ce). 
We can write as 


F /b = Dcos@ —cos’ 0 (3.37) 
with 
D=a/b = —(4K/BM’)(co — c). (3.38) 


Using (3.37) in Fig. 3.50, we calculate and plot F'/b with various values of D as a 
function of 6. Equation (3.38) shows D < 0 for cg > ce and D > 0 for co < Ce. 

We see from this figure that @ = z (antiferromagnetic state) corresponds to the 
smallest free energy when D > 0 (co < c-) as denoted by the dotted circle. As in- 
dicated by the vertical arrow in the right hand side, D decreases as cg increases 
with rising temperature (this is normal thermal expansion plotted by the dashed 
line in Fig. 3.49a). At higher temperatures, D becomes negative. However, antifer- 
romagnetic state (9 = 7) still remains stable as long as D > —2. It is apparent 
from Fig. 3.50 that a first order transition to ferrimagnetic state (@ = 0) occurs at 
D = —2. Similarly, the ferrimagnetic state in the process of decreasing temperature 
shows a transition to the antiferromagnetic state at D = 2. From (3.38), we see that 
Co = Ce + (BM? /2K) at D = —2 and cy = ce — (BM* /2K) at D = 2. This situation is 
summarized in Fig. 3.51. 


+ There is little change of entropy due to change of the molecular field. 
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Fig. 3.50. F/b (F: free energy) plotted as a function of 6 (an angle between the sublattice 
moments) for various values of D(= —(4K/BM7”) (co — cc)) 
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Fig. 3.51. Normal thermal expansion (cg) and magnetostriction (c — cy) = +(BM?/2K) in the 
vicinity of transition temperature. AT denotes the temperature hysteresis. Transition occurs 
when cy = Ce + (BM? /2K) is satisfied (see text). The dotted lines indicates length of the c-axis 
to be measured 


From a relation c, + (BM?/2K) = c-[1 + a(AT/2)], AT is written as 
AT = (B/c-K)M*/a, (3.39) 


where « is a thermal expansion coefficient of co. We see from (3.39) that width of 
hysteresis is proportional to M?. 

We have explained qualitatively the first order antiferromagnetic—ferrimagnetic 
transition of Mn2_,Cr,Sb in the vicinity of the first order transition temperature. 
Next, we compare the result mentioned above with experimental data. From (3.35), 
the jump of the c-axis at T = T, (up or down), Ac/c, is approximately written as 


A[(c — ¢0)/co] & Al(e — co) /¢e] = (B/Kec)M?. (3.40) 
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Fig. 3.52. Relation between the square of sublattice moment M? and the jump of the c-axis at 
the first order transition temperature 7;. (Drawn by using experimental data shown in Figs. 3.46 
and 3.47) 
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Fig. 3.53. Temperature hysteresis of magnetization for Mno_,Cr,Sb (x = 0.07). AT ~ 15K 


M? and A(c—co)/¢c, which is simply expressed by Ac/c, at T = T; can be determined 
from experimental results shown in Figs. 3.46 and 3.47. We note that magnetization 
of a sublattice M has been assumed to be a half of the jump of magnetization at T). 

Figure 3.52 shows experimental result corresponding to (3.40). The linear rela- 
tionship between Ac/c and M* shown in Fig. 3.52 indicates that (3.35) is satisfied. 
From the slope of the dotted line in Fig. 3.52 and (3.35), we can evaluate (B/Kc-) as 
1.9 x 10-°. Applying this value to (3.39), we obtain 


AT = 1.9 x 10-°M?/a. (3.41) 


Now we show temperature variation of magnetization in Fig. 3.53 to show an ex- 
ample of the hysteresis AT. Note that it is somewhat difficult to determine the cor- 
rect value of AT from experiments, because the hysteresis depends on the rate of 
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temperature change and the sample quality, etc. In addition, a large coexistence re- 
gion of low temperature antiferromagnetic phase and high temperature ferrimagnetic 
phase has been found by X-ray diffraction measurements. 

From Fig.3.53, we see that AT ~ 15K. Inserting AT = 15K and M? ~ 
150(emu/g)* (see Fig.3.46) in (3.41), we can estimate as @ ~ 2.5 x 10-4 K7! 
though direct measurement of & is not available. This value of & is larger than @ of 
a similar compound Mn2Sbo.375ASo.125 [77], however, usual one in the intermetallic 
compounds. 

We note that there is no other model to explain the first order antiferromagnetic— 
ferrimagnetic transition of Mnz_,Cr,Sb. The model described above is based on the 
assumption that the exchange parameter J, connecting different sublattices is sensi- 
tive to length of the c-axis and changes its sign at the critical length c,. As shown in 
Fig. 3.49, J. is the summation of J; and J2. If J; and J2 have different sign and com- 
parable magnitude, sign of J, may change due to even small change of the c-axis. It 
is noted that J; and Jz mean exchange interactions between rather distant Mn atoms; 
the nearest distance is 5.05 A for J; and 3.94A for J> as shown in Fig. 3.48. It is 
expected to investigate J; and Jz on the basis of band structures. 


3.7.2 Magnetic Transition of Fe,_,Mn,.As (a ~ 2) 


As for compounds written as Feg_,Mn,As, both Fe2As and Mn2As show antiferro- 
magnetic ordering. For a = 2.1 and 1.25 <x < 1.5, however, ferrimagnetic structure 
occurs in low temperature region. This ferrimagnetism undergoes a first order tran- 
sition to antiferromagnetic ordering with rising temperature [78]. 


3.7.3, Layered Ferromagnets MnAlGe and MnGaGe 


These compounds have an easy axis of magnetization along the c-axis. In these com- 
pounds Mn atoms occupy only the Mn(I)-site shown in Fig. 3.48, and Al (Ga) and 
Ge occupy the Mn(II)-site and the Sb-site, respectively. Mn atoms thus align within 
a two-dimensional plane (001), and the distance between the neighboring planes 
is as long as the c-axis (5.933 A for MnAlGe). These compounds therefore have 
a magnetically two-dimensional character. Probably the RKKY (Ruderman-Kittel— 
Kasuya-—Yoshida) interaction and super-exchange interaction via Al, Ga, and Ge are 
responsible for the magnetic interaction between Mn atoms in the adjacent layers. 
Curie temperature 7c and magnetic moment per 3d atom in Mn;_,M,AJGe are plot- 
ted against x, respectively, in Figs. 3.54 and 3.55 (M = 3d metals) [79,80]. By making 
use of the band of MnAIGe [81], we can explain the reduction of moment for the case 
of M = Fe, which has one more electron within the rigid band picture. The case of 
Cr substitution is also explained in the similar way. The prominent increase of Tc by 
80K due to 15% Cr substitution, which is shown in Fig. 3.54, contrasts to that the 
Cr-substitution decreases the averaged moment of 3d atom. Such increase of Tc may 
be caused by change of the inter-layer interaction as well as intra-layer interaction. 
MnGaGe shows properties similar to MnAlGe [82]. Band structures and magnetic 
properties of these compounds are also discussed in Part II of this book. 
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Fig. 3.54. Tc of Mnj_,M,AlGe (M = 3d atom) plotted as functions of composition (Repro- 
duced from [79]) 


3.7.4 Application of the First Order Transition of Mnj_,,Cr,Sb 


As shown in Fig. 3.53, the hysteresis of the antiferromagnetic—ferrimagnetic transi- 
tion is rather large. A composite material made of compounds Mn2_,Cr,Sb with two 
slightly different Cr compositions is expected to show a magnetization—temperature 
curve as shown schematically by solid bold lines in Fig. 3.56. 

To explain this figure, let us start from an antiferromagnetic state at the initial 
temperature 7), which is not shown in the figure, however, between es and T;. When 
the temperature is raised, the initial magnetization Mp is zero at temperatures T < 7), 
and then the nonmagnetic state changes discontinuously to the magnetic state with 
residual magnetization M; when temperature is raised up to 7}. M, remains un- 
changed in the temperature range of 7; < T < 7,. With further heating, residual 
magnetization M, similarly increases to M. 

It is thus possible to store three values Mo, M,, and Mp) at one position (bit) 
by different three heating temperatures T’s with fy <T<T,1<T <t, and 
Th < T. Thus, a composite material made of the three compounds with different Cr- 
composition can, in principle, store four values. It may be possible to store more 
data using a composite material made of compounds with continuously different 
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Fig. 3.55. Magnetic moment per 3d atom in Mn;_,.M,AIGe (M = 3d atom). Pare is effective 
moment estimated from the Curie—Weiss constant, and Ua is magnetic moment evaluated from 
spontaneous moment (Reproduced from [79]) 


Magnetization 


es Us Ts 


M, 


composition | 


composition 2 
— 1+2 


T (K) 


Fig. 3.56. Magnetization—temperature curves of composite material made of two compounds 
having different Cr compositions. T; and Ty (or Tj and T;) are the first order transition tem- 
peratures of the compound in heating and cooling processes, respectively. (7; — T/) etc. is 


temperature width of hysteresis 


Cr compositions. This mechanism is superior to an usual magnetic memory, which 
stores two values at a position (bit) corresponding to plus/minus magnetization. We 
need to pay attention that the hysteresis curves in Fig.3.56 are idealized. Actual 
magnetization curves show jump of magnetization within a finite temperature range 
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as shown in Figs. 3.53 and 3.46. It is necessary to reduce this temperature range 
by, for example, homogenization of Cr composition [83]. This compound system 
Mn>_,Cr,Sb may be applicable easily for temperature controller because large vari- 
ety of transition temperatures can be obtained by changing the Cr composition. 
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Appendix: Magnetic Transition and Free Energy 


Free energy G of a magnetic material can be written in the Landau type expansion 
form as 


G=G(0)+a07 + ayo" +.a60° + ..., (Al) 


where o is a magnetic moment of a ferromagnet (or a sublattice moment for an an- 
tiferromagnet) normalized by the value at T = 0 K. When an external magnetic field 
is applied, Zeeman energy —O# is added. This expansion form is irrelevant to the 
nature of electrons, localized or itinerant; however, physical meaning of the coeffi- 
cients depends on the character of electrons and theoretical model. In this Appendix, 
we show an expansion form of magnetic entropy, which appears in the Bean—Rodbell 
theory of magnetic phase transition of MnAs, up to the sixth order o°. Further, we 
show expression of G(o) and its temperature dependence. 

Let us consider a paramagnetic material consisting of N particles with spin j in 
a magnetic field H. In this case free energy is written as G = —gjUpoHN —TS(o). 
Replacing S(o) by NkgS(o), we have G/NkgT = —(gjUpo/kgT)o — S(o) = 
—ao —S(o). Applying the minimum condition dG/do = 0 to the final equality, 
we obtain 


a =—dS(o)/do. (A2) 


On the other hand, @ is related to o through the Brillouin function with spin j, 
Bj (a), as 


o =B,(a) =Aa+Bo3+Da>+..., (A3) 
where A, B,D,... are given by 


A= (1/3)[(24+1)?—1]/(2/)?, B= —(1/45)[(27 + 1)* - 1/(2i)%, 
D= (1/945)[(2j +1)®—1]/(2))*.... (A4) 
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By expanding S(o) in (A2) with respect to o, we have 


S(o) = S(0) +407 + bo* +co®+--- 
=In(2j+1)+a0*+bo*+co®+--. (AS) 


As the number of microscopic states, in the case of o = 0, is given by Z = (27+ 1)" 
and S = kg InZ, it is reasonably understood that S(0) becomes In(2'+ 1). By inserting 
(A5) into (A2), we have a = —(2a0 + 4bo0? + 6co* + -:-). By inserting this relation 
into (A3) and setting the coefficient of each o” term zero, we have the expansion 
coefficients a,b,c,... as 


a=-1/A<0, b=B/4A* <0, 
c = —(1/2)(B’/A’) + D/6A® <0, (A6) 


where A, B, and D are given by (A4). We insert (A5) into the expression of magnetic 
entropy (3.2) of Sect. 3.2.1, and then we put a condition of a minimum energy with 
respect to volume 0G/dV = 0 into (3.2) again. Thus, we obtain the free energy G 
also minimized by crystal volume. By setting T/T) = t (the relative temperature), 
P =O, and H = 0, we have 


G(o) /NkgTo = —tln(2j+ 1) —a(t—1)0? —b(t—)o* —cto®+---, (AT) 
with 
n = —(a*/2b)NkgKToB" > 0. (A8) 


Figure 3.57 displays [G(o) — G(0)|/Nkg7o calculated by (A7) with parameters n = 2 
and j = 3/2, which were chosen so as to be applicable to MnAs. 


(GCa) —G(0))/NKT, 


Fig. 3.57. Free energy with parameters j = 3/2 and n = 2 vs. normalized moment o drawn 
for various relative temperatures 
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On heating, the minimum point of free energy disappears at temperature tf = 
1.083 as seen in Fig. 3.57, which means that magnetization o disappears discontinu- 
ously due to the absence of energy barrier. Therefore, this is of the first order and also 
illustrated in the o-T curve with 7 = 2 in Fig. 3.5. While in the process of decreas- 
ing temperature from about t = 1.083, the peak of free energy disappears at ¢t = 1. 
This correspond to sudden appearance of magnetization. We note that the curves in 
Fig. 3.57 are somewhat inaccurate in the region o ~ | because of truncation of the 
expansion formula of entropy. (Numerical calculation is necessary to evaluate en- 
tropy without the expansion form.) The coefficient of the o+-term of free energy has 
an important meaning. When 7 > 1, this coefficient is negative in the temperature 
range | <t < 17. As coefficients of o” and o° terms are positive for t > 1, the free 
energy does not have minimum at finite o unless the coefficient of the o*+ term is 
negative. Therefore, the first order transition takes place only for n > 1. 


Part II 


Itinerant Electron Theory 


4 


Electronic Band Structure and Magnetism 
of NiAs-Type Compounds 


4.1 Band Calculation 


Extensive studies on band structures of NiAs-type transition-metal compounds 
started in the middle of 1980s by the group of Motizuki, who is one of the present 
authors, of Osaka University and the group of Haas of the University of Groningen. 

For the band calculations, the group of Osaka University employed the self- 
consistent augmented plane wave (APW) method, while the group of the University 
of Groningen used the self-consistent augmented spherical wave (ASW) method. In 
spite of the difference in the method used for calculations, the general features of 
calculated bands are similar to each other. In the APW calculations carried out by 
the group of Osaka University, muffin-tin approximation was used for atomic poten- 
tial. For exchange and correlation interaction, the local (spin) density approximation 
(LDA or LSDA) was adopted within the formula of Gunnarsson and Lundqvist (GL). 
The equation of von Barth and Hedin (HB) for the LDA was also used in some cal- 
culations. In this book, we show mainly the results calculated by the group of Osaka 
University. 

In recent years, we have applied more sophisticated methods. For example, we 
carried out band calculations by using the self-consistent full-potential linearized 
augmented plane wave (FLAPW) method to obtain more accurate band structures. 
Kulatov and his co-workers carried out band calculations by using the self-consistent 
linear-muffin-tin orbital (LMTO) method. Nakada and Yamada recently calculated 
bands by using the LMTO method with atomic sphere approximation (ASA); this 
method which employs localized basis functions is convenient for the calculations 
for a system including many atoms in the unit cell. 

Scalar relativistic effects are usually included, however, spin—orbit interaction 
(SOD) is neglected in many band calculations. The SOI removes the degeneracy in 
the dispersion curves and also induces mixture between up-bands and down-bands, 
resulting in reduction of magnetic moment of ferromagnetic state. Furthermore, the 
SOI in the band calculations is needed to investigate optical processes. Recently, the 
SOI has been included in the calculations by Kulatov et al. 
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Generalized gradient approximation (GGA) has been used to improve the LDA. 
The GGA has enabled us to see difference of total energies between different struc- 
tures. For example, total energy calculation for CoAs having the NiAs- and the MnP- 
type structures has been carried out by Nakada, and his results by the full-potential 
LMTO method with the GGA show good agreement with experiments. 


4.2 Band Structures and Optical Properties 


4.2.1 Pnictides: MnAs and MnSb 
Nonmagnetic Band of MnAs and MnSb with NiAs-Type Structure 


Figure 4.1a, b shows the crystal structure and the first Brillouin zone of the NiAs- 
type structure, respectively. Figure 4.2a shows the dispersion curves of nonmagnetic 
MnAs calculated by the self-consistent APW method with the muffin-tin approxi- 
mation and the LDA [1]. The two bands in the low-energy region consist of As-4s 
orbitals. Above the gap, there are 16 mixing bands consisting of As-4p and Mn-3d 
orbitals. Figure 4.2b shows the density of states (DOS) of the p—d mixing bands. 
It is seen from this figure that there are two remarkable peaks: the low-energy re- 
gion is contributed by bonding orbitals of As-4p and Mn-3d, and the peaks in high- 
energy region by antibonding orbitals mainly consisting of 3d states of Mn. As seen 
in Fig. 4.2b, the width of the antibonding p—d mixing bands is about 5 eV (or 0.4 Ry), 
which suggests that d states of MnAs are not localized but are of itinerant character. 
Density of states of MnSb has features similar to those of MnAs, however the band 
width is slightly smaller than that of MnAs. 


O anion 


Fig. 4.1. (a) The NiAs-type crystal structure and (b) the first Brillouin zone 
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Fig. 4.2. (a) Dispersion curves and (b) density of states of the nonmagnetic MnAs calculated 
by the APW method 


As shown in Fig. 4.2b, the Fermi level is near the large peak of the density of 
states, which is contributed mainly by the d orbitals of Mn. The density of states at 
Ep, p (Er), of MnAs is larger than that of NiAs-type compounds with 3d elements 
other than Mn, which is consistent with the fact that both MnAs and MnSb are fer- 
romagnetic below the Curie temperatures Tc = 318 K and 537K, respectively. 

These densities of states calculated by the APW method agree well with the 
bands of MnAs and MnSb calculated by Sandratskii by using the Green’s function 
method [2]. Recently, Nakada and Yamada have calculated nonmagnetic bands of 
MnAs having the NiAs-type structure by the FLAPW method with the GL formula 
and the LMTO-ASA method with the BH formula. Their results shown in Fig. 4.3a, 
b are rich in fine structure, however the general features of the bands are similar to 
the bands obtained by using the APW method shown in Fig. 4.2b. 


Ferromagnetic Bands of MnAs and MnSb with the NiAs-Type 


Among the NiAs-type compounds, only MnAs, MnSb, and MnBi are ferromagnetic 
as shown in Part I. We have calculated ferromagnetic bands using the self-consistent 
APW method. The muffin-tin approximation for atomic potential and the LSDA 
within the GL formula for the exchange and correlation interaction were used. 

The density of states of the p—d mixing bands of ferromagnetic MnAs is shown 
in Fig. 4.4. Due to spin-polarization, the mixing bands of As-4p and Mn-3d orbitals 
split in the up-spin and the down-spin subbands, whereas no polarization occurs 
in the 4s bands in the low-energy region. In Fig. 4.4, the broken and dot-dashed 
curves show components of Mn-3d and As-4p orbitals in the muffin-tin spheres, 
respectively. Comparing Fig. 4.4 with the nonmagnetic bands shown in Fig. 4.2, we 
see that the p—d mixing bands are not shifted rigidly, but deformed significantly due 
to the spin-polarization. 
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Fig. 4.3. Density of states of nonmagnetic MnAs calculated using the (a) FLAPW method and 
(b) LMTO-ASA method (K. Nakada, Dr. thesis) 


We evaluated total magnetic moment from the calculated band as 3.1 Up/formula. 
This value is somewhat smaller than the observed magnetic moment 3.4{p/formula 
(see Part I). Magnetic moment in the muffin-tin spheres of a Mn-site and an As-site 
is 3.12Up/formula and —0.15uUp/formula, respectively, hence most of the magnetic 
moment is attributed to Mn atoms. In addition to the magnitude of the total magnetic 
moment, the calculated values of magnetic moments at Mn- and As-sites correspond 
well to the results of neutron diffraction measurements, therefore we see that mag- 
netic moments in MnAs are well understood from the band theory. 

Ferromagnetic band of MnSb with the NiAs-type structure is similar to the 
band of MnAs except that the width of p—d mixing band is slightly smaller than 
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Fig. 4.4. Density of states of ferromagnetic MnAs calculated by the APW method. The broken 
and dash-dotted lines are Mn-3d and As-4p components, respectively 


that of MnAs. According to the band calculation, the total magnetic moment is 
3.2Up/formula; magnetic moment inside the muffin-tin sphere at Mn- and Sb-sites 
are 3.41 {p/formula and —0.12pp/formula, respectively. These values of magnetic 
moment agree fairly well with those measured by Yamaguchi et al. [3]. 

Results of our band calculation for MnSb agree well with the bands calculated 
by Coehoorn et al. by using the self-consistent ASW method [4], namely both dis- 
persion curves and density of states in these two works agree well with each other, 
and the magnetic moments by Coehoorn et al. are 3.31p at a Mn-site and —0.06 Lup 
at an Sb-site. Density of states at the Fermi level p (Ep) is obtained by us to be 37.0 
(states/Ry unit cell), which agrees well with 35.8 (states/Ry unit cell) calculated by 
Coehoorn et al. These values are also in good agreement with 32.6 (states/Ry unit 
cell) evaluated from measurements of specific heat. 

X-ray photoelectron spectroscopy (XPS) measurements for MnSb, CoSb, and 
NiSb were carried out by Liang and Chen [6]. Coehoorn et al. compared the XPS 
spectrum with their ferromagnetic bands (the 5s band and the p—d mixing band) of 
MnSb. As shown later in Fig. 4.17 in Sect. 4.2.5, there is good correspondence for 
the peaks of the XPS spectrum and density of states (see also Sect. 4.2.5). 
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Fig. 4.5. Density of state of ferromagnetic MnAs. (a) FLAPW method, (b) LMTO-ASA 
method [5] 


Nakada and Yamada calculated ferromagnetic band of MnAs with the NiAs- 
type structure by using both the FLAPW method and the LMTO-ASA method [5]. 
Calculated densities of states are shown in Fig. 4.5a, b. Compared with the band 
obtained by the APW method, these densities of states are rich in fine structures, 
however, their general features are similar to those calculated by the APW method. 
Nakada and Yamada also evaluated contributions of five Mn-d orbitals and three 
As-p orbitals to the APW wavefunctions. Their results are shown in Fig. 4.6, where 
five d-components for density of states and their values (or numbers of electron be- 
longing to the respective components) are depicted as a function of electron energy. 
Contributions of five d-components to charge and magnetic moment are also shown 
in Table 4.1, where the values are normalized so that the summation of contributions 
from the five d orbitals is 100%. These calculations revealed that Mn-d orbitals are 
not extended in the c-plane, but have rather spherical distribution. They also plot- 
ted charge distribution in the (0001) plane and (1100) plane which include Mn—Mn 
bonds. The results also confirmed that charge density of Mn-d orbitals has spherical 
character. Furthermore, they clarified by plotting charge density in the (1120) plane 
including Mn and As atoms that the Mn—As bond is the strongest bond. The charge 
distributions are shown in Fig. 4.7. 
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Fig. 4.6. Density of states of ferromagnetic MnAs calculated by the FLAPW method. (a) 


Integrated DOS and (b) components of each d orbitals (K. Nakada, private communication) 


Table 4.1. Contribution rates of d orbitals for charge and magnetic moment, which are inte- 
grated up to the Fermi level, of MnAs 


327-1? xx zy x y xy 
Charge (%) 22 19 19 20 20 
Moment (%) 20 18 18 22 22 


They also investigated the effect of crystal volume change on magnitude of mag- 
netic moment. They carried out band calculations in two cases: one is the case of 
fixed-axial ratio c/a (i.e., the bond angle Mn—As-Mn is kept constant) and another 
is the case of unfixed-axial ratio (i.e., the bond angle is variable). They found that 
magnitude of magnetic moment is significantly affected by the change of bond an- 
gle. This result indicates that magnetic moment is sensitive to lattice distortion in the 
direction of the a-axis. On the other hand, effect of lattice distortion along the c-axis 
on magnetic moment is small, though Mn—Mn distance in this direction is smaller 
than the distance along the a-axis. These results of band calculation are consistent 
with the experimental data seen, for example, in Figs. 3.3 and 3.11. 


Nonmagnetic Band of MnAs with MnP-Type Structure 


MnAs shows a structural transformation to the MnP-type structure at 7; = 318 K (see 
also Fig. 3.2). We calculated band structure of MnAs with the MnP-type structure 
by using the self-consistent APW method. The band structure deforms due to the 
lattice distortion accompanied with the structural transformation, which results in 
decrease of p (Ep). In Sect. 4.3, we will show in detail that anomalous behaviors of 
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Fig. 4.7. Charge distribution of ferromagnetic MnAs (a) (0001) plane, (b) (1100) plane, (c) 
(1120) plane. In each figure, (i), (ii), and (iii) show majority spin, minority spin, and moment, 
respectively [5] 
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paramagnetic susceptibility and thermal expansion observed for MnAs (see Fig. 3.2) 
are well explained by considering spin fluctuation effect and the deformation of the 
band structure. 


4.2.2 FeAs, CoAs, and NiAs 


It is known that CoAs shows structural transformation at 7, = 1,250 K from the MnP- 
type structure (low-temperature phase) to the NiAs-type structure (high-temperature 
phase). FeAs has the MnP-type structure, and NiAs has the NiAs-type structure.! 
FeAs has the double-helical magnetic ordering lower than Ty = 77 K, whereas CoAs 
and NiAs are nonmagnetic at all temperatures (see Table 2.7). 


Nonmagnetic Bands of CoAs and NiAs with the NiAs-Type Structure 


We calculated band structures of CoAs and NiAs with the NiAs-type structure by 
using the self-consistent APW method within the muffin-tin potential approxima- 
tion and the LDA. Density of states of the p—d mixing band of CoAs is shown in 
Fig. 4.8. Density of states at the Fermi level is p (Ep) = 46 (states/Ry unit cell) for 
CoAs and 31 (states/Ry unit cell) for NiAs. These values are much smaller than 156 
(states/Ry unit cell) of MnAs, which is consistent with the fact that CoAs and NiAs 
are nonmagnetic but MnAs is ferromagnetic. 

We also depicted Fermi surfaces of CoAs and NiAs, and then found that CoAs 
has two hole surfaces around the A-axis and these surfaces are expected to have good 
nesting for a displacement of wavevector I'M. Since this wavevector I'M is relevant 
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Fig. 4.8. Density of states of nonmagnetic CoAs with the NiAs-type structure calculated by 
the APW method 


' We have related argument in Sect. 4.4. 
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to deformation into the MnP-type structure, the NiAs-type structure of CoAs has an 
instability to the MnP-type structure. This fact corresponds well to the high 7, of 
CoAs. On the other hand, the structural transformation to the MnP-type structure 
will not occur for NiAs because no nesting of Fermi surfaces is expected [7]. We 
have related discussion in Sects. 4.4 and 4.7. 


Nonmagnetic Bands of FeAs and CoAs with the MnP-Type Structure 


Figure 4.9 shows DOS curves for the p—d mixing band of nonmagnetic CoAs 
and FeAs with the MnP-type structure. These compounds have characteristic DOS 
curves, namely, the densities of states show abrupt increase just below the Fermi level 
[7]. We will show in Sect. 4.3 that anomalous temperature dependence of paramag- 
netic susceptibility of CoAs and FeAs can be well explained by the spin fluctuation 
theory including such feature of the densities of states. 

It is known that total energies evaluated from band structures make it possible 
to know which crystal structure is most stable at 0K. Recently, Nakada and co- 
workers calculated nonmagnetic bands of CoAs with the NiAs-type and the MnP- 
type structures by using the full-potential LMTO method with the GL formula for 
the exchange and correlation interaction. From the calculated total energies they 
found that the MnP-type structure is more stable than the NiAs-type structure for 
CoAs. This result is consistent with observation that CoAs remains in the MnP- 
type structure up to 1,250 K. They obtained also similar result by calculation of full- 
potential LMTO method with the GGA term. 
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Fig. 4.9. DOS curves calculated by the APW method for nonmagnetic CoAs and FeAs having 
MnP-type structure 
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4.2.3 CrSb, CrAs, and CrP 


These compounds take the NiAs-type or the MnP-type structure [8]. As shown in 
Table 4.2, the lattice constants of these compounds generally decrease as the an- 
ion changes in the order Sb, As, and P. In particular, change of the a-axis is re- 
markable. The CrX compounds show various magnetic properties with variation of 
anion atom. CrSb is an antiferromagnet with Cr moment of 3.0up [9, 10]. CrAs 
takes a double-helical magnetic ordering, and the magnetic moment of Cr is 1.7Up, 

which is only about a half of Cr-moment of CrSb [11]. CrP is a paramagnetic com- 
pound whose magnetic susceptibility is almost independent on temperature [12]. The 
double-helical magnetic ordering of CrAs has been known to disappear by a few per- 
cent substitution of P atoms for As [13]. The experimental results mentioned earlier 
suggest that there is a close correlation between the crystal structure (or lattice pa- 
rameters) and magnetic properties in the CrX compounds. 

To investigate the correlation between crystal structure and magnetic properties 
of the CrX compounds, we carried out first-principle band calculations by using the 
full potential linear muffin-tin orbital (FP—LMTO) method [14]. The crystallographic 
parameters in Table 4.2 were used in the calculations. As for the local density approx- 
imation (LDA) for exchange and correlation interaction, we have used the formula 
of Vosko, Wilk, and Nussair [15]. We carried out iterative calculations using 125 
k-points in the 1/24 of the first Brillouin zone so that muffin-tin potential converges 
within 10-7 mRy. 

Figure 4.10 shows densities of states of nonmagnetic CrP, CrAs, and CrSb. From 
these densities of states we see following three features: 


1. Total band width increases in the order CrSb, CrAs, and CrP. 

2. Hybridization between Cr 3d-orbitals and anion orbitals takes place in both high- 
and low-energy regions, but not in the intermediate region. 

3. Width of the large DOS peak near the Fermi level in CrSb, where there is little 
Cr—X mixing broadens and splits roughly in two groups as X goes from Sb to P. 


Table 4.2. Crystallographic parameters of CrX (X = P, As, or Sb) at the room tempera- 
ture. Lattice constants a, b, and c are common to both NiAs- and MnP-type structures (see 
Sect. 1.3). u, v, w, and x denote displacement parameters [8] 


CrP CrAs CrSb 
Crystal structure MnP-type MnP-type NiAs-type 
a(A) 3.114 3.463 4.127 
b (A) 6.018 6.212 7.148 
c (A) 5.360 5.649 5.451 
u 0.05 0.05 0 


v 0.065 0.05 0 
w 0.018 0.006 0 
x 0.0073 0.0065 0 
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Fig. 4.10. Local densities of states for CrP, CrAs, and CrSb (/Ry atom) 


As shown in Table 4.2, the smallest distance between Cr atoms in the c-direction, 
which is equal to c/2 (~2.7A) in the case of the NiAs-type structure, is almost com- 
mon to all the CrX compounds. On the other hand, length of the a-axis which cor- 
responds to the distance between Cr atoms within the c-plane, depends considerably 
on X of the CrX compounds. This difference in the length of the a-axis influences 
the bonding character in the c-plane, resulting in the characteristics of the density of 
states especially in the vicinity of Ep as mentioned earlier. 

We evaluated contribution of five d-orbitals to the total density of state. The re- 
sults for CrSb and CrAs are shown in Fig. 4.11. As shown in Table 4.2, CrSb and 
CrAs take different crystal structure. The length of the a-axis of CrAs is much smaller 
than that of CrSb. We discuss relation between such difference of the a-axis length 
and magnetic properties. As shown in Fig. 4.11, the density of states of CrSb has a 
sharp peak near the Fermi level. This part consists of mainly d(xy) and d(x? — y’) 
components which extend in the c-plane and do not mix with p states of Sb as shown 
in Fig. 4.10. On the other hand, states near the Fermi level of CrAs consist of not 
only d(xy) and d(x” — y*) but also d(xz) and d(yz) orbitals. Because of the short 
a-axis of CrAs compared with that of CrSb, width of d(xy)- and d(x* — y*)-bands 
become larger for CrAs, and then the Fermi level is pushed down to the low-energy 
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Fig. 4.11. Contribution of 3d orbitals on density of states of CrSb and CrAs 


side of the peak of the density of states. The atomic displacement of As which is 
denoted by v in Table 4.2 due to the MnP-type distortion gives rise to lowering the 
energy of d(xz) and d(yz) components, however, these states remain in the energy 
region higher than the Fermi level. The low-energy (high-energy) peak is bonding 
(antibonding) band consisting of d(xz), d(yz), d(3z? — r7) states and p state of anion 
as seen in Fig. 4.10. 

It is important that the change of the a-axis, which induces change of d(xy) and 
d(x* —y*) bands, plays an important role to explain the variety of magnetic properties 
seen in CrSb, CrAs, and CrP. 
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Fig. 4.12. Density of states of Cr for one of the two magnetic sublattices of antiferromagnetic 
CrSb (in unit of states/Ry atom) 


From the density of states of CrP shown in Fig.4.10, p(E£p) is evaluated as 
27.3 (states/Ry f.u.), which is 1/2 and 1/4 of the value of p(Epf) estimated from 
measurements of specific heat [16] and magnetic susceptibility measurement [17], 
respectively. 

Figure 4.12 shows density of states for one of the two magnetic sublattices of an- 
tiferromagnetic CrSb, which gives Cr magnetic moment of 2.7{1p. This value agrees 
well with experimental values 2.7 Up [10] and 3.0 Up [9]. Since there is little p—d 
mixing in the region near the Fermi level as shown in Fig. 4.10, there is strong po- 
larization of Cr-d states; most of the up-spin states lies above the Fermi level. This 
fact indicates that Cr-d states in CrSb has a considerable localized character probably 
because of the large length of a-axis. It is also notable that the down-spin band has a 
gap-like shape near the Fermi level. 

Calculated total energies for paramagnetic, ferromagnetic, and antiferromagnetic 
states of CrSb are plotted in Fig. 4.13 as a function of crystal volume. As seen in this 
figure, the band calculation indicates that the antiferromagnetic state is most stable. 
This is consistent with observation that CrSb is an antiferromagnet. 

To understand magnetic properties of CrAs (double helix, Uc; = 1.7Up, Tn = 
250K, see Sect. 3.4 for detail), it is expected to investigate wavevector-dependent 
susceptibility y(q). This is an important problem left for the future. 


4.2.4 Chalcogenides: CrTe, CrSe, and CrS 


Compounds between Cr and a chalcogen (S, Se, or Te) also take the NiAs-type 
structure. CrTe is a ferromagnet with 7c = 340K, whereas CrSe and CrS take an- 
tiferromagnetic ordering below Ty = 285 K and 460K, respectively. Effects of high 
pressure on these compounds have attracted much interest. From electron-spin res- 
onance (ESR) measurements, Shanditsev pointed out that ferromagnetism of CrTe 
disappears at the critical pressure 2.8 GP [18]. Kanomata and his co-workers inves- 
tigated pressure dependence of Curie temperature 7c and magnetic moment [19]. 
They pointed out that the ferromagnetic state is not the ground state under high 
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Fig. 4.13. Total energy of CrSb as a function of unit cell volume. The total energies for para- 
magnetic state, ferromagnetic state, and antiferromagnetic state are plotted. Vexper denotes 
experimental value of unit cell volume 


pressure. According to observation by Ishiduka and Eto, the ferromagnetism of CrTe 
disappears at 7 GP and a structural transition to the MnP-type structure occurs at 
13 GP [20]. 

To explain the results of these high-pressure experiments, we have theoretically 
studied the possible magnetic transitions by investigating pressure dependence of 
band structures. We will discuss the pressure effect in detail in Sect.4.5, and in 
this section we describe the results of band structures calculated by Nakada and 
co-workers; they investigated dependence of bond length (Cr—Cr) and bond angle 
(Cr-Te—Cr) on charge distribution and magnetic moment in ferromagnetic CrTe. 

On the basis of the FLAPW band, calculations for ferromagnetic CrTe, Nakada 
decomposed the APW wavefunction into the components from three p orbitals of Te 
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and five d orbitals of Cr. The results are shown in Fig. 4.14a, b and Table 4.3. The 
density of states and integrated density of states are shown in Fig. 4.14b and 4. 14a, 
respectively. 

Contribution from each d orbital to charge and magnetic moment in CrTe is 
shown in Table 4.3, where the values are normalized so that the summation of con- 
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Fig. 4.14. (b) Density of state of ferromagnetic CrTe calculated by the FLAPW method. (a) 
Integrated DOS component from each d orbital (K. Nakada, private communication) 


Table 4.3. Contribution rates of respective d orbitals to charge and magnetic moment in CrTe 


Bear ge ee oP ey oy 
Charge (%) 22 19 19 20 20 
Moment(%) 21 18 18 21.5 21.5 
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tributions from five d orbitals is 100%. We see that all d orbitals have almost equal 
contribution to charge and magnetic moment, which means that charge of the d states 
of a Cr atom distributes spherically as also seen in MnAs. This result that charge dis- 
tribution of Cr is spherical has been confirmed by depicting charge distribution in the 
(0001) and (1100) planes which include Cr—Cr bond. They also found strong Cr-Te 
bond from charge distribution in the (1120) plane which include Cr and Te. They 
have clarified that the strongest bond between Cr and Te is the bond in the (1120) 
plane. 

They also investigated the effect of volume change on magnetic moment. They 
calculated magnetic moment as a function of volume for the two cases: (i) bond angle 
of Cr-Te-Cr is fixed (i.e., the axial ratio c/a is fixed) and (ii) bond angle of Cr-Te— 
Cr is unfixed (i.e., c/a is variable). The calculated results indicate that magnetic 
moment is sensitive to lattice change along the a-axis where Cr—Cr distance is large, 
and insensible to lattice change along the c-axis where Cr—Cr distance is small. This 
is a situation similar to that of MnAs. 


4.2.5 Optical Properties 


We can discuss optical properties of compounds on the basis of band calculation 
including the spin-orbit interaction. 

Kulatov and co-workers calculated ferromagnetic bands of MnAs and MnSb by 
using the LMTO method including the spin-orbit interaction [21]. In the calculations, 
they used the muffin-tin approximation and the expression of von Barth and Hedin 
(BH) for the exchange and correlation interactions, and they also took the spin-orbit 
interaction into account. 

Figure 4.15a, b shows the calculated densities of states of MnAs and MnSb, 
respectively. Except the split of the band due to the spin-orbit interaction, overall 
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Fig. 4.15. Densities of states of (a) ferromagnetic MnAs and (b) MnSb. (A) shows total DOS. 
(B) and (C) show partial DOS of Mn-d and As(Sb)-p, respectively (Reproduced from [21]) 
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features of the bands are similar to the bands calculated by Motizuki and her 
co-workers with the self-consistent APW method. Magnetic moment evaluated here 
by the LMTO method agrees well with experimental value. 

On the basis of these bands, Kulatov investigated optical properties of MnAs and 
MnSb. He obtained a dielectric function by making use of the Kubo formula for the 
linear response theory. Optical absorption coefficient is proportional to the imaginary 
part €)(@) of the dielectric function, and €)(@) consists of two terms: the Drude term 
e? (@) arising from intra-band transition, and the €(@) term arising from interband 
transition. In the dipole approximation, eb (@) is given by 


872 2 : 
€3() = Base, NAAR) Psa — fea) 9 (Exat— Ex, —fiw), (4.1) 


Sona eu sh he eh ; 
where f,, is the Fermi—Dirac distribution function, j = —i—V, and A is a band 
m 


suffix. The matrix elements for the dipole transition (kA’|j|kA) were calculated from 

the LMTO wavefunctions. Interband optical photoconductivity o(@) is related to 
b 

€;(@) as 


eb(w) = > 0\®): (4.2) 


Figure 4.16a, b shows the photoconductivity o(@) calculated for ferromagnetic 
MnAs and MnSb, where magnetic moments are assumed to align along the c-axis 
(z-axis). In these figures, Oy,(@) and 6;,(@) are plotted by the solid and broken 
curves, respectively. Due to anisotropy arising from the hexagonal symmetry of the 
crystal structure, shape and peak positions are different for O,, and 0;z. 

The features of o(@) of MnAs and MnSb in Fig. 4.16a, b are similar even in 
the structures. The main peaks of o(@) are attributed to transition of electrons as 
follows: 


e The peak at around f@ ~ 0.3 eV is attributed to direct interband transition which 
has become allowed owing to mixing between the up-spin and down-spin bands 
due to the spin-orbit interaction. 

e The peaks at around f#@ ~ | eV correspond to interband transition between down- 
spin bands of Mn. 

e The broad and strong peak in the energy region of hw = 2.5 ~ 4.5 eV is mainly 
due to transition between p-band of Sb (As) and d-band of Mn. 


In Fig. 4.17, we show XPS spectrum of MnSb measured by Liang and Chen (by 
the dashed curve) [4], and the calculated density of states of ferromagnetic MnSb 
below Fermi level (the solid curve). Since the measurements were carried out with a 
poly-crystal sample, the experimental spectrum corresponds to 0;, + 20;,. There is 
rough correspondence between the experiment and the calculation, however, experi- 
mental data on single-crystal are necessary for detailed comparison. 
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Fig. 4.16. Optical conductivity for (a) ferromagnetic MnAs and (b) ferromagnetic MnSb cal- 
culated by Kulatov. The solid and dashed curves show 0, and 0;z, respectively 


4.3 Spin Fluctuations and Anomalous Magnetic 
and Elastic Properties 


4.3.1 Paramagnetic Susceptibility and Anomalous Thermal Effect of MnAs 
and MnAs,_ ,P,. 


With increasing temperature, MnAs shows a first-order transition from ferromag- 
netic state to paramagnetic state at Tc = 318K. ? At this temperature, the structural 
transformation from the NiAs-type to the MnP-type structure also occurs with 2% of 


2 See Sect. 3.2 for experimental detail of MnAs, MnP and their mixed compounds. 
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Fig. 4.17. Density of states for ferromagnetic MnSb (solid curve) and XPS spectrum (dashed 
curve). The dotted curve shows DOS broadened with a Gaussian of full with at half height = 
0.6eV (Reproduced with permission from [4], Copyright (2009) by the American Physical 
Society) 


volume decrease. In the temperature region above Jc, reverse change to the NiAs- 
type structure takes place at 7; = 398 K (see Fig. 3.3. The NiAs-type and MnP-type 
structures are also shown in Fig. 4.38). In the temperature region, T > T;, paramag- 
netic susceptibility obeys approximately the Curie—Weiss law. According to obser- 
vation carried out by Ido et al., the y~!-T curve shows a slight upward convex in this 
temperature range. Magnitude of magnetic moment evaluated from the y~!-T curve 
in the high-temperature region is 3.8 Ug. This is larger than saturate moment in the 
ferromagnetic state 3.4 Up. 

In the intermediate temperature range (Jc < T < 7T;) where the crystal structure 
is of the MnP-type, the y~!-7 curve shows anomalous behavior that the curve shows 
a broad peak, and crystal volume increases rapidly with increasing temperature (see 
Figs. 3.2 and 3.3, etc.). Crystal volume at J; is almost equal to the volume just below 
Tc. Such anomalous behavior of susceptibility and thermal expansion in the MnP- 
type phase are more remarkable in MnAs j_,P, (0 < x < 0.275). CrSbj_,As, and 
Cro.4Mno.6As in the MnP-type phase also show such anomalous behavior as shown 
in Part I. According to neutron diffraction measurements by Schwartz [22], magnetic 
moment of MnAs and MnAs;_,P, in the MnP-type phase is smaller than the moment 
in the NiAs-type phase. They also found that magnitude of magnetic moment shows 
smooth change with temperature. 

According to measurements of magnetization process carried out by the group 
of Osaka University (Fig.3.4 and Fig.3.15-3.17, etc.), magnetic moment of 
MnAs\_,P,, depends little on magnetic fields at temperatures below 7c, whereas 
it depends considerably on magnetic fields at around 7;. In the MnP-type phase in 
the temperature region just above 7c, field-induced metamagnetic magnetization 
curve with hysteresis was observed. As shown in Fig. 3.4, temperature dependences 
of magnetization measured in 20 and 400 kOe revealed that magnetic moment in 
20kOe shows a discontinuous change at 7c. This change is considered due to the 
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structural phase transition at 7c. On the other hand, magnetization in 400 kOe con- 
tinuously changes with temperature, which suggests that MnAs remains to be of 
the NiAs-type structure even in the temperature region above 7c. These results 
mentioned earlier indicate that the metamagnetic behavior is accompanied with the 
structural transition from the MnP-type phase (low-fields) to the NiAs-type phase 
(high-fields). 

Yamaguchi et al. carried out neutron diffraction measurements for the ferromag- 
netic state of MnAs [3]. They reported that a magnetic moment 0.23Up exists at an 
As-site and it is antiparallel to magnetic moment at Mn-site. Their neutron diffrac- 
tion measurements also revealed that an Sb atom in MnSb has moment 0.3 Ug in the 
opposite direction of the moment of Mn atom. 


4.3.2 Spin Fluctuation and Magnetism 


In the itinerant model, correlation between electrons is a very important factor to 
discuss magnetic properties in finite temperature. Moriya and his co-workers de- 
veloped the theory in which the correlation is treated in terms of fluctuation of spin 
density. They succeeded in the unified theory in describing both the extreme case 
of weak ferromagnetism, where the spin fluctuation has a small amplitude and local 
character in the g-space, and another extreme the case of local moment, where all 
q-modes of spin fluctuation and their mode—mode coupling are important [23-25]. 

Motizuki and Katoh investigated magnetism of MnAs [26, 27] by applying the 
spin-fluctuation theory within the formalism of Moriya—Takahashi [24] and Usami-— 
Moriya [25]. 

We start from the following single-band Hubbard Hamiltonian: 


FE = K+ KH 


U 
= Y tea act, Fe - 1s} ; (4.3) 
oj. ri 


where the symbol f;¢ denotes a transfer integral, a\, a creation operator which creates 
an electron with spin o ata site j. nj and S; in #% are operators of charge density and 
spin density, respectively. U and J are an effective intra-atomic Coulomb integral 
and an effective intra-atomic exchange integral, respectively. Assuming that five d 
orbitals degenerate, U and J are related to U’ and J’ (Coulomb and exchange integrals 
for a single orbital Hubbard Hamiltonian, respectively) as 


1 


U = =(9U' —4J') 


Nn 


1 
J = —(U'+4J'). 
By using the Stratonovich—Hubbard transformation, we can transform this many- 
body problem into a single-particle problem under a magnetic field () and an elec- 
tric field (7)) fluctuating temporally and spatially. Within the static approximation, 
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a partition function Z (or free energy F’) can be written in the functional integral 
form as 


Zeit ee) ete (4.4) 
and 


e Pa — [TY agan exp(—B¥(E.n)). (4.5) 


Here, LU is the chemical potential, N a total electron number operator, and B = 1 /kpT. 
The functional ‘¥ [6,1] = ‘¥o[€,n]+ ¥i[&, n], is given by 


¥ols l= 5 > [8F +0)| (4.6) 


and 


B 
e PF ilS.n] — (row - { at, {a8 -$)(2) eames) ) (4.7) 
where 


2, = \/4nJ/B, > = /—nU/B, (4.8) 


and €; and 7; denote local magnetic and electric fields, respectively, operating at the 
jth site. In addition, we introduce a quantity 


1 f 
Xa = 5B 2ibea5-ao (@—=x9,2); (4.9) 


where &, is a Fourier component of the fluctuating magnetic field €; and No is the 
number of lattice points in the system. We note that pat means that the term g = 0 is 
excluded. xq in (4.9) is a square average of local amplitude of the §-field, and related 
to a square average of local amplitude of spin fluctuation as 


( (Sia - (Sja))) = “ (x. a aa) (4.10) 


By assuming the mode—mode coupling of spin fluctuation to be local and applying 
a saddle-point approximation for the charge field, Moriya and the co-workers intro- 
duced the following form of the functional: 


W,[é] = ~221 YY Kao Slab st NoE be Eo. (4.11) 
q @ 
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The first term of (4.11) represents the nonlocal part, where terms independent 
on relative direction of spin fluctuation are neglected. An appropriate function of g is 
assumed for Xgq[x, 0]. The second term of (4.11) represents the local part, for which 
Usami and Moriya obtained an expression for L[x, 9] by using the coherent potential 
approximation (CPA). 

From the functional ‘¥; thus obtained, the partition function of the system is 
calculated and then various quantities such as paramagnetic susceptibility and local 
magnetic moment can be evaluated successfully. 

We show here the results of calculations carried out for paramagnetic MnAs and 
MnAs1_xPx [26,27]. First, we simplified the density of states calculated by the self- 
consistent APW method by replacing the curves with a set of straight lines and then 
obtained model density of states as shown in the lower part of Fig. 4.18. We use the 
model density of states for the present calculations. For example, the curve denoted 


by (1) in Figs. 4.18 and 4.19 show calculated local magnetic moment (S5) and 


inverse paramagnetic susceptibility y~', respectively. They are calculated with the 
model density of states (1) and in the case that electron number per orbital n is 
0.65 and U and J are equal to 2.8 (in unit of 5/6eV). As shown in Fig. 4.19 the 
calculated y~'-T curves show the Curie-Weiss behavior with weak upward convex. 
Magnetic moment estimated from the y~!-T curve is 2.58, which is slightly larger 
than the value 2.30 tug evaluated from the model density of state and the Hartree— 
Fock approximation at 0K. As seen from Figs. 4.19 and 4.20, there is thus good 
correspondence between the present calculations and the experimental results except 
for the numerical values. For example, Curie temperature estimated from an equation 
x~! =0 is only 1/10 of the observed value. It is needed to point out that the relation 
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Fig. 4.18. Calculated local magnetic moments myo, using the model density of states (1)-(4). 
The model density of states (1) corresponds to the NiAs-type MnAs, (2)-(4) correspond to 
MnAs distorted to the MnP-type structure (Reproduced with permission from [26]) 
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Fig. 4.19. Inverse paramagnetic susceptibility calculated for the model densities of states (1)— 
(4). The dotted and dash-dotted curves represent y~! in the intermediate temperature range 
To <T < T with T, = 2Tc and T, = 37c, respectively (Reproduced with permission from [26]) 
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Fig. 4.20. Temperature dependence of (Ap — ARP) calculated for the model densities of states 
(1)+(4) shown in Fig. 4.18. The dotted and dash-dotted curves represent (Ap — AHF ) in the 
intermediate temperature range Tc < T < T, with 7; = 27¢ and T, = 37¢, respectively (Repro- 
duced with permission from [27]) 
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yt = 0 gives Jc for the case of a second-order transition and the transition at Tc 
of MnAs is of a first order. We may expect a better estimation of Tc by considering 
effect of degeneracy of orbitals and quantum effects which are neglected in the spin- 
fluctuation theory. 

Next, we investigate paramagnetic susceptibility of MnAs in the intermediate 
temperature range of Tc < T < J where the crystal structure is of the MnP-type. 
As noted in Sect. 4.2, the APW band calculations have revealed that the height of 
the sharp peak of density of state near the Fermi level in the NiAs-type structure 
becomes smaller with lattice distortion to the MnP-type structure. In order to take 
such change of band structure into account, we modified the model density of states. 
The model density of states in the MnP-type structure is shown in the lower part 
of Fig.4.18 for the three cases denoted by (2), (3), and (4), which correspond to 
the cases where p(Ep) is 90%, 85%, and 80% of p(Ep) of (1) in the NiAs-type 
phase, respectively. The solid curves denoted by (2), (3), and (4) in Figs. 4.18 and 
4.19 display the temperature variations of the local magnetic moment and inverse 
paramagnetic susceptibility calculated for the model densities of states (2), (3), and 
(4), respectively. It is seen from Figs. 4.18 and 4.19 that the local magnetic moment 
decreases in the order (1), (2), (3), and (4), resulting in steeper inclination of the 
x—'-T curves. In order to explain the observed y~!-T curve of MnAs, in Fig. 4.19 
we start from a point on the curve (1) in the high-temperature region and consider the 
x_—'-T curve in the process of lowering temperature. With decreasing temperature, 
x! decreases first along the curve (1), and in the MnP-phase below the transition 
temperature 7, ye changes to the values of points on the curves (2), (3), and (4) 
successively because the lattice distortion increases as temperature decreases. The 
dotted- and the dashed-curves show temperature dependences of y~! for two cases 
of G) T = 27¢ and ii) T; = 37¢, respectively. For the case (i), we assumed that the 
lattice distortion changes monotonically with temperature in the temperature range 
between 7; and Jc. For the case (ii), the monotonical change of distortion in the 
temperature region between 7; and 1.5 7c, and constant distortion in the temperature 
region between 1.5 Jc and Tc are assumed. Measurements of lattice distortion for 
MnAs and MnAso,8Po.2 has revealed that MnAso.gPo.2 takes the MnP-type structure 
in wider temperature range than MnAs (see, for example, Fig. 3.19 and [28]). It also 
has been found that temperature dependence of the lattice distortion of MnAso gPo.2 
is large only in a narrow temperature region just below 7, but almost constant in 
other temperature regions below 7, [28]. Thus, the y~!-7 curve calculated for the 
case (i) corresponds to MnAs and that for the case (ii) to MnAso.gPo.2. The ab 
curves of these compounds show a peak in the temperature region of the MnP-type 
phase. The peak of MnAsp gPo 2 is especially remarkable (see, for example, Figs. 3.2 
and 3.13). 

As mentioned earlier, we have shown that the anomalous temperature depen- 
dences of susceptibility observed for MnAs and MnAs;_,P; compounds are rea- 
sonably understood by the spin-fluctuation theory including the effect of lattice 
distortion of the MnP-type structure. 
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4.3.3 Spin Fluctuation and Elastic Properties 


We describe in this section how to explain anomalous elastic properties of MnAs and 
MnAs\_,P, in the MnP-type phase from the viewpoint of itinerant electrons taking 
the effect of spin fluctuation into account [27]. 


OF 
Cohesive pressure of electron system is given by P = —=— with the free en- 


ergy F’. For simplicity we assume that electronic energy & is proportional to band 
width W and that only W depends on volume. The pressure P in this approximation 
is written as 


P=AA/V. (4.12) 
In this equation, A is defined by A = — ve and A is given by 
W oV 
OF 
A=) &——. 4.13 
2 * Oe ( ) 


From free energy evaluated by using the spin-fluctuation theory described in 
Sect. 4.3.2, we calculated the quantity A for paramagnetic NiAs-phase and para- 
magnetic MnP-phase. The free energy per atomic site is written as a summation 
of Fo and AF, which are terms irrelevant and relevant, respectively, to the spin 
fluctuation. The quantities Fo and AF are given by 


1 
Fo = ——~ Yin] 1 + cP H-&) 4.14 
0 nop =|! +e (4.14) 
and 
_ nl? 
AF =n +o + EG), (4.15) 


where the nonlocal term has been omitted. By using an expression of the local term 
L(x) obtained within the coherent potential approximation (CPA), we can calculated 
the quantity A from (4.13). The parameter Ap, which means A of the paramagnetic 
state, is given by 


2 
Ay =—= fdef(e +n) xIm{(e-+ 2, %) Fo}, (4.16) 
where f is the Fermi—Dirac distribution function and € energy of electrons, and other 


symbols in (4.16) are explained next. We note that |7; |? and x in the expression of AF 
are assumed not to depend on & explicitly. The chemical potential is determined by 


2 
-= | defle +n) ImF, =n, (4.17) 
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where n is electron number. The symbol ¥% is nU/2, Fp is a diagonal element of 
single-particle Green’s function, and 2, denotes site-independent self-energy. Within 
the Hartree-Fock approximation in which the spin fluctuation is neglected, the local 
magnetic moment is zero and Xp = 0. In this case, Ap becomes AY given by 


An = 2 fae fle +Y%)ep(é), (4.18) 


where p(é€) is the density of states. Contribution of the spin fluctuation to A is given 
by (Ap =A," ). Change of pressure AP is related by the compressibility « to change 
of volume AV as 

AV 

— = KAP (4.19) 

Vv 

We see from (4.12) that AP is written by AA, therefore volume change in (4.19) is 
given by 


AV 
a A KA(Ap—ABP)/V +AKAATT/V 


AV, AVo 
=F + v~ (4.20) 
The first and the second terms denote volume changes arising from the spin fluctua- 
tion and that from change of electronic energy, respectively. 

Calculated (Ap ao ) is plotted against temperature in Fig. 4.20. The curves 
(1)-(4) are results calculated for the model densities of states (1)—(4) shown in 
Fig. 4.18, respectively. Since (Ap =A ) arises from the spin fluctuation, temper- 
ature dependence of this quantity reflects temperature dependence of local ampli- 
tude of spin fluctuation plotted in Fig. 4.18. As seen in Fig. 4.20 that (Ap =a ) 
at each temperature decreases in the order (1), (2), (3), and (4), which means that 
spin fluctuation is suppressed as the lattice distortion to the MnP-type increases. To 
explain anomalous thermal expansion observed for MnAs, etc., (see Figs. 3.3, 3.19, 
etc., in Part I) we consider temperature dependence of (Ap Ar ) in the process of 
decreasing temperature; (Ap — AHF) decreases first along the curve (1), and in the 
temperature region of 7c < T < T, below 7j, it changes successively to the points 
on the curves (1) (2) —(3) —(4) as the lattice distortion increases. The dotted- 
and dashed-dotted-curves denote temperature dependence of (Ap =A for the case 
T, = Tc and T, = 37¢, respectively. These results clearly indicate that the spin fluctu- 
ation plays an important role on the remarkable volume change with temperature for 
MnP-type MnAs and MnAsj_,P,, because temperature dependence of (Ap — AS ) 
gives rise to the volume change. 

Next, we investigate the volume change at 7c. Using values k = 4.5 x 
10-!'m2/N, V = 34 x 1077°m?, and A = 5/3 (estimated by Heine [29]), we eval- 
uated at Tc the quantity AV; /V from (Ap =A) to be 3.4%, 4.8%, and 5.8% for 
the changes of (Ap =A5) between curves (1) and (2), (1) and (3), and (1) and (4), 
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respectively. We also estimated AV>/V comparable to AV, /V from the tight-binding 
calculation for the MnP-type structure, which means that decrease of local magnetic 
moment due to the structural transition and energy change of electron system have 
comparable contributions to the discontinuous volume change at 7c in MnAs (see 
Fig. 3.3 in Part I). Calculated values of AV /V is 7.1% for the change (1) (2), 9.5% 
for (1) (3), and 12% for (1) (4), respectively. Experimental volume change at 
Tc is 2% for MnAs and 8 ~ 10% for MnAs)_,Py. 

Finally, we note that we estimated A ~ 1.1 from our band calculation for MnAs. 
If we use this value of A instead of Heine’s result A = 5/3, the calculated values of 
AV /V become smaller by about 2/3 than the value mentioned earlier. 


4.3.4 Paramagnetic Susceptibility of CoAs and FeAs 


The compounds FeAs and NiAs crystallize in the MnP-type and the NiAs-type 
structures, respectively, while CoAs shows the structural phase transition at 1,250 K 
from the MnP-type to the NiAs-type structures as temperature rises. Below this tem- 
perature, the crystal structure of CoAs becomes the MnP-type. FeAs takes the MnP- 
type structure at all temperatures, and NiAs remains in the NiAs-type structure ex- 
cept at very temperatures. 

Morifuji and Motizuki calculated paramagnetic susceptibility vy of these com- 
pounds by using the spin-fluctuation theory and nonmagnetic bands calculated by 
the APW method [30]. It has been revealed that 7 of NiAs is almost temperature- 
independent, while 7’s of CoAs and FeAs, which have the MnP-type structure, 
shows anomalous temperature dependence. According to the experimental data [31], 
inverse susceptibilities of CoAs and FeAs obey the Curie—Weiss law in high temper- 
ature region and increase for further decreasing of temperature. 

The temperature dependences of paramagnetic susceptibility of CoAs and FeAs 
are similar to that of a weak ferromagnetic semiconductor FeSi. Takahashi and 
Moriya explained anomalous temperature dependence of susceptibility of FeSi on 
the basis of the spin-fluctuation theory and characteristic shape of density of states 
[32]. According to their theory, negative mode—mode coupling is thermally induced 
when the Fermi level lies at a dip of density of state and the density of states just 
below the Fermi level has a large value. In this situation, the negative mode—mode 
coupling gives rise to a dip of the y~'-T curve and a rapid increase of y~! in the tem- 
perature region below the dip. Furthermore, the Curie-Weiss behavior of susceptibil- 
ity occurs in the temperature region above temperature of the dip. The Curie—Weiss 
behavior is brought about by the saturation of local amplitude of spin fluctuation. 

Morifuji and Motizuki calculated v —!_T curves for CoAs and FeAs by taking 
account of the spin-fluctuation effect, as well as the fact that the densities of states 
for FeAs and CoAs shown in Fig. 4.9 are similar to that of FeSi [33, 34]. 

We show in Fig. 4.21 the calculated results for CoAs. We used the model density 
of states, shown in the inset, in which the calculated curve is replaced by fragments of 
straight lines. We see that the temperature dependence of inverse susceptibility y~! 
agrees qualitatively with that of measurement plotted by x [31]. This result indicates 
that the anomalous temperature dependence of paramagnetic susceptibility y~! of 
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0 500 T(K) 


Fig. 4.21. Inverse paramagnetic susceptibility of CoAs calculated by using the model density 
of states inserted in the figure is plotted by the solid curve. The observed susceptibilities at 
several temperatures are shown by x 


CoAs is well understood in terms of the spin-fluctuation effect except discrepancy in 
magnitude of y~! and the temperature of the dip. We have obtained similar results 
for FeAs. 


4.4 Fermi Surface of NiAs and the de Haas—van Alfen Effect 


We can investigate Fermi surface by observing the de Haas—van Alfen effect. Since 
the cross section normal to an applied magnetic field is related to period of the 
de Haas—van Alfen oscillation, it is possible to investigate shape of Fermi surfaces 
from the de Haas—van Alfen measurements carried out in various directions of ap- 
plied magnetic fields. 

Recently, various techniques and theories for band calculation have also been 
developed to calculate more accurate band structures. As for the correlation between 
electrons, the local density approximation, which is usually used, may be insufficient 
because average electron density is used. Recently, improvement of the LDA such 
as GGA has been carried out. It is thus important to justify calculated bands by 
examining the shape of Fermi surfaces. 

The Group of Kamimura measured the de Haas—van Alfen oscillation using sin- 
gle crystals of NiAs and NiSb with the NiAs-type structure [35]. Since the NiAs-type 
structure has two metal atoms and two anions in the unit cell, the dispersion curves 
near the Fermi level have a complicated dependence on wavevector due to mixing 
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between d-orbitals of Ni and p-orbitals of As or Sb. Many branches are observed in 
the experiments, and therefore it is difficult to reproduce the shape of Fermi surfaces 
from de Haas—van Alfen measurements. 

Recently, however, Harima developed a computer program to calculate de Haas— 
van Alfen frequencies from Fermi surfaces obtained by band calculations [36]. It 
has enabled us to easily make comparison between calculated band structures and 
de Haas—van Alfen measurements. To make direct comparison, we first constructed 
Fermi surfaces of NiAs from results of band calculation and determined dependences 
of the de Haas—van Alfen frequencies in the direction of magnetic field, and then 
compared them with the experimental data [35,37]. 

We calculated initially band structure of nonmagnetic NiAs and NiSb by us- 
ing the APW method, and later the FLAPW method to obtain more accurate band 
structures. We applied the LDA with the GL formula for the exchange and corre- 
lation interaction. Scalar relativistic effects were taken into account and the spin— 
orbit interaction was omitted in the calculations. The program codes TSPACE and 
KANSAI-92 were used for these calculations. 

Dispersion curves of NiAs and NiSb are shown in Fig. 4.224, b, respectively. The 
bands in the low-energy region consist of As-4s or Sb-5s orbitals. The 16 bands 
above the gap are p-—d mixing bands consisting of Ni-3d and As-4p (or Sb-5p) 
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Fig. 4.22. Dispersion curves of nonmagnetic band of (a) NiAs and (b) NiSb (Fig. 4.22a is 
reproduced form [35]) 
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Fig. 4.23. Density of states of nonmagnetic band of (a) NiAs and (b) NiSb (Fig. 4.23a is 
reproduced form [35]) 


orbitals. We show the total densities of states of the p—d mixing bands p(£) with 
the components of Ni-d orbitals and As- (or Sb-) p orbitals for NiAs and NiSb in 
Fig. 4.23a, b, respectively. The density of states in both compounds can be divided 
into three main parts: p—d bonding band, nonbonding band consisting mainly of d or- 
bitals in the middle region of energy, and p—d antibonding band in the highest region 
of energy. The gross feature of the density of states is similar to that of MnAs. Total 
width of the p—d mixing band is larger than that of MnAs shown in Figs. 4.2 and 4.3. 
The Fermi level is an energy far above the energy region with the large peak of the 
density of states, p(Ep) = 29.12 (states/Ry unit cell) for NiAs and 14.53 (states/Ry 
unit cell) for NiSb are obtained. These values are much smaller than 156 (states/Ry 
unit cell) of MnAs, which is consistent with the fact that MnAs is a ferromagnet but 
NiAs and NiSb are paramagnetic. 

The value of p (Ep) of NiAs evaluated from measurement of temperature depen- 
dence of specific heat is 34.54 (states/Ry unit cell) [37], which is somewhat larger 
than the calculated value of 29.12, however, the difference between experimental 
and calculated values in NiAs is smaller than that in MnAs. This fact suggests that 
electron—electron interaction and electron—phonon interaction in NiAs is weaker than 
those in MnAs. Width of the density of states below the Fermi level for both NiAs 
and NiSb is about 6eV (~0.44 Ry) and there is a sharp peak at 2 eV (~0.15 Ry) en- 
ergy lower than the Fermi level. Such character of the density of states corresponds 
well to the results of photoemission measurements [38]. 

Figure 4.24 shows calculated Fermi surfaces of NiAs. NiAs has three Fermi sur- 
faces: two hole surfaces with cylindrical shape around the [’A-axis and an electron 
surface whose shape is somewhat complicated. 
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Fig. 4.24. Fermi surfaces of nonmagnetic state of NiAs (Reproduced form [35]) 


The de Haas—van Alphen frequencies calculated from the Fermi surfaces are plot- 
ted as a function of angle of direction of magnetic field in Fig. 4.25. The branches 
labeled e, m, and d in Fig. 4.25 seem to correspond generally well to the experimental 
branches labeled €, 1, and 6 as shown in Fig. 4.26 [35]. 

However, the a and f branches in Fig. 4.26 cannot be explained by the present 
calculations, and the c- and n-branches which are seen in the theoretical curves are 
not measured. The absence of observations of the c- and n-branches may be ex- 
plained by a possible structural transformation due to nesting effect between the 
hole surface around the ’A-axis and the electron surface around the KH-axis. Recent 
X-ray diffraction measurements have revealed that NiAs shows a structural transfor- 
mation at a very low temperature [39]. Due to this structural transformation, it is 
considered that the hole surface around the I’A-axis and the electron surface around 
the KH-axis disappear. 

Figures 4.27 and 4.28 show the Fermi surfaces and the calculated de Haas— 
van Alfen frequency of NiSb, respectively [37,40]. The results of measurements 
are shown in Fig. 4.29. As seen in NiAs, there is good correspondence between the 
e branch in Fig. 4.28 and the € branch in Fig. 4.29. Unlike the case of NiAs, the 
a- and the B-branches which correspond to the cj- and the ;-branches in Fig. 4.28 
are observed in NiSb. This result corresponds well to the fact that NiSb does not 
show structural transformation, but remains in the NiAs-type structure even in low 
temperature. 
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Fig. 4.25. Calculated de Haas—van Alfen frequencies of nonmagnetic NiAs (Reproduced 
form [35]) 


4.5 Pressure Effect on Magnetic State of CrTe, CrSe, and CrS 


CrTe becomes a ferromagnet below 7c = 340K. By application of high pressure, 
Tc decreases remarkably and the magnetic moment also decreases, which suggests a 
kind of pressure-induced phase transition. 

To clarify the effects of high pressure on magnetic properties of CrTe, we carried 
out band calculations for nonmagnetic, ferromagnetic, and antiferromagnetic states 
for various values of the unit cell volume. In these calculations, the axial ratio c/a 
is fixed to the value at room temperature. For the antiferromagnetic calculation, the 
sublattice magnetic moments are assumed to be in the c-plane and are antiparallel 
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Fig. 4.26. Experimental de Haas—van Alfen frequencies of NiAs (Reproduced form [35]) 


a_ NiSb, hole surface c NiSb, electron surface 


Fig. 4.27. Fermi surfaces of nonmagnetic NiSb 
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Fig. 4.28. Calculated de Haas—van Alfen frequencies of NiSb 


between adjacent c-planes. We employed the APW method for the first calculations, 
later we used the FLAPW method to obtain more accurate band structures [41-43]. 
In this section, we show the results obtained by the FLAPW method. The density of 
states of ferromagnetic CrTe calculated by the FLAPW method is shown in Fig. 4.14. 

In Fig. 4.30, we show the calculated total energy as a function of the lattice con- 
stant a for nonmagnetic, ferromagnetic, and antiferromagnetic states in CrTe [41]. 
The total energy for the ferromagnetic state takes the minimum value at a = 3.970 A, 
which agrees well with the experimental value a = 3.981 A. We can express the total 
energy curve of ferromagnetic state by a function 


FEtotal = Eo + bo(a — ag)? + b3(a — ag)* + ba(a — ao)’, (4.21) 


where Ey = —31357.17 (Ry/unit cell), ag = 3.970, bz = 0.471, b3 = —0.575, and 
b4 = 1.686. 
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Fig. 4.29. Experimental de Haas—van Alfen frequencies observed for NiSb 
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Fig. 4.30. Total energy curves for ferromagnetic, antiferromagnetic, and nonmagnetic CrTe as 
a function of lattice constant a (Reproduced from [41]) 


4.5 Pressure Effect on Magnetic State of CrTe, CrSe, and CrS 111 


Since the APW method gives the minimum energy at the lattice constant a = 
4.180 A, we see that the FLAPW method brings about improvement of calculated 
results. Generally, a full potential calculation reproduces more correctly the exper- 
imental values of lattice constant than the APW method. As shown in Fig. 4.30, the 
ferromagnetic state is most stable at a = 3.981 A, which is consistent with the fact 
that CrTe is a ferromagnet under atmospheric pressure. 

It is seen in Fig. 4.30 that the antiferromagneic state has smaller energy than that 
of the ferromagnetic state for the lattice constant smaller than a = 3.58 A, which 
indicates that pressure-induced phase transition from the ferromagnetic state to the 
antiferromagnetic state can occur. The inset of Fig. 4.31 shows relation between cal- 
culated electronic pressure and the lattice constant a. From this curve, we see that the 
pressure-induced magnetic transition may occur at a critical pressure of about 40 GP 
which corresponds to the lattice constant 3.58 A. The APW calculation gives critical 
pressure 20 GP. In Fig. 4.31 total magnetic moment of Cr and Te in the ferromagnetic 
CrTe are plotted against the lattice constant a. Large part of the magnetic moment 
arises from Cr atom, however, Te atom has a small negative magnetic moment in the 
opposite direction. (see also Sect. 4.2.4). 

The total magnetic moment of CrTe is calculated to be 3.1up/formula at a = 
3.970A. This value is not so different from observed moment 2.29uUp/formula. The 
APW calculation gives magnetic moment 3.9,p/formula. 

We can evaluate electron pressure P from the equation P = —dEqotq,/dV. Since 
the unit cell volume is given by V = 1.35a°, the pressure P is expressed as P = 
—0.247a~* dEtqotqi/da. The bulk modulus, which is given by B = —dP/dlogV 
(V = Vig), is calculated to be B = 404kbar, or the compressibility (B-!) 2.43 
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Fig. 4.31. Magnetic moments in ferromagnetic CrTe plotted against the lattice constant a. 
Inset: Electronic pressure as a function of lattice constant (Reproduced from [41]) 
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Fig. 4.32. Total energy for ferromagnetic, antiferromagnetic, and nonmagnetic CrSe (Repro- 
duced from [41]) 


x 107? kbar~!. The calculated compressibility evaluated by the FLAPW calculation 
is slightly smaller than the value 2.6x 10 ~* kbar—! evaluated by the APW calcula- 
tion. Pressure coefficient of magnetic moment is calculated as OM/dP = —1.33x 
10~* Up kbar~! at a = 3.981 A, and then the magnetic moment decreases at the rate 
k =—0InM/0P = —4.21 x 10-3 kbar~!, which is small but in the similar order in 
comparison with the observed value —12 x 10-3 kbar~! [42]. 

We also carried out the total energy calculations for CrSe and CrS. As an exam- 
ple, we show the relation between Eyota) and the a-axis of CrSe in Fig. 4.32. The total 
energy calculations show that the antiferromagnetic state is most stable even under 
high pressure both for CrSe and CrS. 


4.6 Magnetic Ordering and Instability of Paramagnetic State 
4.6.1 Double-Helical Magnetic Ordering of MnP-Type Compounds 


The total energy calculation, as described in Sect. 4.5, makes it possible to clarify 
which magnetic state occurs in a material with itinerant electrons. It is, however, a 
hard task to find out the most stable magnetic state from the calculations of total en- 
ergies, especially in a material with a lot of magnetic atoms in a unit cell because we 
need to carry out total energy calculations for many possible magnetic orderings. In 
addition, some of transition metal compounds such as CrAs explained in Part I have 
helical magnetic ordering. If wavevector of the helical ordering and lattice vectors 
are incommensurate, it is impossible to carry out band calculation. 
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An alternative way to search for a possible magnetic state is to study instability of 
paramagnetic state. In this section, we describe how we can investigate the instability 
of paramagnetic state on the basis of paramagnetic band calculations as well as the 
Hubbard model [44]. 

The transition metal pnictides such as FeAs and CrAs with the MnP-type crys- 
tal structure show the double-helical magnetic ordering. Magnetic transition tem- 
perature is 77 K for FeAs and 250K for CrAs. We describe later the procedure to 
investigate instability of the paramagnetic state for the case of FeAs. 

As shown in Fig. 4.33, the unit cell of FeAs (orthorhombic MnP-type structure) 
contains four Fe atoms and four As atoms. Paying attention only to Fe atoms and 
neglecting the difference in As atoms surrounding the respective Fe atoms, we can 
take a reduced unit cell which contains only two Fe atoms. The thick lines in Fig. 4.33 
shows the reduced unit cell containing two Fe atoms labeled as “a” and “b”. 

Experimental studies have revealed that the helical ordering of FeAs is described 
by a wavevector Q = (0,0,0.375 x 27/c) (See also Tables 2.4 and 2.6). We note that 
Q of 0.375 x 27/c corresponds to wavelength 2.67b. We also note that the lattice 
vector c here corresponds to b given in Table 2.6. It also has been known that the 
phase difference between magnetic moments at the a-site and the b-site in the same 
unit cell is @ = 140°, and magnitude of the magnetic moment of the Fe atom is 
0.5up [31]. 


Fig. 4.33. The magnetic ordering of FeAs with the MnP-type structure. The thick lines denote 
the reduced unit cell (see the text). f; and f) denote transfer integrals between Fe atoms at 
different sublattices. f3 is a transfer integral between Fe atoms at the same sublattices. See 
also Fig. 2.1 and Table 2.6 


3 @ is denoted by 69 in Part I. 
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4.6.2 Instability of Paramagnetic State 


We start from a Hubbard Hamiltonian with transfer integrals and intra-atomic 
Coulomb interaction. We assume that each Fe atom has a single orbital. For the 
transfer integral, we consider three types of interaction: tz between the same sublat- 
tice (a—a or b—b) and ft, and f) between the different sublattices (a—b) as shown in 
Fig. 4.33. 

By Fourier transformation and the Hartree-Fock approximation for the Coulomb 
interaction, we have the following expression of the Hubbard Hamiltonian: 


w= 2 In (k)a} ko + Tj (k)b} ako + T3(k) (Gace ns biabko) 
oO 


+UYY |Ag-c4h goto + Bq-o4j sgobko] — NUD, [Ag+Ag- + Be+Bq-| 
kq © q 


-uy, [Apap tq Mee + BEDE, bi +ATgGtt, deg +Bi,bf Pesa-| 
kq 


+NU [IAG +1877], (4.22) 
q 
where 
1 + 1 } 
Ago = N Yi %o%+9o ) Bao = N ie bk+a0 )» 
k k 
ba eg 1 + + 1 + 
xq — N DiMth thea »P+q — N Ze? 1 Okge 4 (4.23) 
and 


T, (k) =r Se + tr Yel i—ri-9) 
j j 

T3(k) = ts Yel ir), (4.24) 
j 


In (4.24), r; denotes position of an a-atom in the ith unit cell, and r; + T position 
of a b-atom in the jth unit cell. 

Using this Hamiltonian, we investigate instability of paramagnetic state at OK 
against an external magnetic field modulated by a wavevector q. In order to do it, 
first we add Zeeman terms to the Hamiltonian. Next, by solving equations of mo- 
tion of the quantities Ago, etc., we have the following relation between q-dependent 
magnetic fields and magnetic moments: 
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ee 7 ee | x = 70) ea mes 
H’(q) x(q) x(a) ) \MPg) M’(q) 

In this equation, M“(q) = Ug(Ag— — Ag+) and M?(q) = Up(Bg— — Bg+) are Fourier 
components of spin density at the a- and b-sites, respectively. H“(q) and H(q) 
denote g-dependent magnetic fields operating on M“(q) and M?(q), respectively. 


x"“(q), etc., are elements of the inverse magnetic susceptibility tensor. These quan- 
tities are given by 


1 G(g)+U [Gi(@)* -|B(@)/"] 


aa ~bb 
= = F (4.26a) 
0 =O) =F @?- Be 
x? (a) = x"**(q) = = Big 7 (4.26) 
2up Ti(q)* — |B(@)| 
where 
nwa 010% (fe BBR OB CK ok 4.27 
1) =a Xo ( W+X ( DM+X ( W+Xo ( Q) ( . a) 
k 
and 
Bq 2 Ty (k+q) T,(k) 
~ aN IT (k+ 91/10) 
x [x8 (ka) + 16° (ka) — x5" (ka) — x8 Cka)]. (4.270) 


In these expressions, @ and B denote bands, and 5 (kK. 4); etc., are bare electronic 
susceptibilities which are given by 


f(Exu) — f(Ex+qv) 
Exu = Ex+qv 


wv 
XxX = 


where E. i denotes the band energy obtained by diagonalizing the Hamiltonian. When 
paramagnetic state is unstable against an infinitesimal q-dependent magnetic field, a 
relation 


det[Z(q)] = X(a) — 1x°"(@) |’ = 0 (4.28) 


must be satisfied. Substituting (4.26a)—-(4.27b) into (4.28), we can evaluate a critical 
value of the Coulomb interaction U, as a function of the wavevector q. If U > Us, 
paramagnetic state is unstable, resulting in occurrence of a magnetic order. When U, 
takes the minimum value for a wavevector g = Q, we may say that a magnetic order 
described by the wavevector Q is most likely to realize. 

Substituting H“(q) = H’(q) = 0 into (4.25), we have a relation 7““(q)M“(q) + 
x (q)M°(q) = 0. When q = Q, we see that |7““(q)| is equal to |7“(q)|. Therefore, 
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absolute value of M“(Q) is equal to the absolute value of M?(Q). In the real space, 
magnetic moments at the a- and the b-sites are written as M“(r) = M“(Q) e~'27 and 
M?(r) = M°(Q) ei2", respectively. 

From these expressions, the phase difference @ between the moments of the a-site 
and the b-site is defined by the following equation: 


bi, b(O) | 
Mitt) _ M'Q) 19-2 = gio, (4.29) 


M“(r;) — M*(Q) 


Since 7““(Q) is a real number and the phase factor of 7“?(Q) is derived from 
Ih(Q)*, @ the phase difference between moments at a-site and adjacent b-site is 
given by 


@ = m-+arg[I3(Q)] + arg [e12"] (4.30) 
We carried out calculations using parameters t, = —0.07Ry, f2/|t1| = —1.0, 
t3/|t,| = —0.2, and electron number n = 1.1. These values were chosen so that they 


reproduce as nicely as possible the features of Fermi surfaces and density of states 
of FeAs which were obtained by the APW band calculation. 

The calculated results are shown in Figs. 4.34 and 4.35. As shown in Fig. 4.34, 
U, takes the minimum value for Q = 0.4 x 277/c, which means that the paramagnetic 
state is most unstable against magnetic ordering modulated with this wavevector. We 
see that this value of Q agrees well with the experimental value Q = 0.375 x 27/c. 
Furthermore, as shown in Fig. 4.35, the phase difference @ takes 158° at Q = 0.4 x 
21/c. This value of @ is also close to the experimental value @ = 140°. 


2n r 22 1 Qn 
c b a 


Fig. 4.34. Critical value of Coulomb interaction U, evaluated as a function of wavevector q. 
As shown by the arrow, Uc takes the minimum for Q = 0.4 x 2m7/c 


4.6 Magnetic Ordering and Instability of Paramagnetic State 117 


Fig. 4.35. Phase difference @ between spins at the a-site and b-site calculated as a function of 
q. For Q = 0.4 x 277/c, @ = 158°. This value of @ is close to the observed value @ = 140° 


4.6.3 Energy of Double-Helical Spin Density Wave State 


Next, we investigate the double-helical spin density wave (DHSDW) state and eval- 
uate electronic energy with the Hubbard Hamiltonian used in the instability study. In 
order to describe a DHSDW state modulated by a wavevector Q, we neglect g-terms 
except the one g = Q in (4.22), and have 


2 
H =F ¥(k,Q)20 (k,Q)¥(k,Q) + = +NU (l4gP x 135") . (431) 
k 


The #(k,Q) and ‘¥(k,Q) in this equation are expressed by the following equations: 


T3(k) Ti(k) —UA_p 0 
ae Ti\(k) T3(k) 0  —UBRy ia 
1 -UAS 0 (kK +Q) T(k+Q) 


0 —-UBS Ti (k+Q) Ts(k+Q) 
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and 
i (k,Q) = fe sl it op. cP (4.33) 
Total energy of the DHSDW state is given by 
4 2 
NUn 
Espw =) Y) Exuf(Emu) + — +NU (lap)? + 18g I") (4.34) 
k u=1 


where Ex, is eigen energies evaluated by diagonalizing the Hamiltonian, f(E,,) the 
Fermi—Dirac distribution function, Ao and Bo Fourier component of spin density at 
the a- and b-atoms, respectively. The phase difference @ between spins at the a-site 
and the b-site is given by 


BS. 
2 it — gid, (4.35) 
AO 


We obtained Ao and Bo by self-consistent calculations, and evaluated the total en- 
ergy Espw at 0 K as a function of Q and U. 

In Fig. 4.36, Espw is plotted as a function of wavevector Q for the two values 
of U/|t;|. The total energy Espw takes the smallest value at Q = 0.4 x (27/c). This 
result corresponds to the result of the instability study described earlier. 


Espw 


N 
0.3 
Ul t:|=8 

0.2 
—0.2 


—0.3 


—0.4 


UN t1=6 


—0.5 


r 21 
@ “e- 
Fig. 4.36. Energy of the DHSDW state calculated as a function of wavevector Q for two values 
of U/|t,|. Corresponding to the result of the instability study, the energy takes the minimum 
value at Q = 0.4 x 2m/c 
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Fig. 4.37. Upper part: Espw — Epara plotted as a function of U/|t;|. For the curves of each 
value of Q, in the region of U /|t;| larger than the value denoted by the arrow, Espw is smaller 
than Epara- Lower part: lA plotted as a function of U/|t,| calculated for Q = 0.4 x 27/c 


We show energy difference between the DHSDW state and paramagnetic state, 
Espw — Epara, at 0 K in Fig. 4.37. The vertical arrows indicate the values of U, where 
Espw is equal to Epara for the respective values of Q shown in the figure. These 
values of U corresponds to U, shown in (4.28). Since U, takes the smallest value 
for Q = 0.4 x (2z/c) as shown in Fig. 4.37, the calculated energies of the DHSDW 
state is consistent with the result of instability study. The phase difference @ between 
spins in the same unit cell is evaluated from (4.35) as @ = 153. The value of lAg| for 
Q =0.4 x (27/c) is shown in lower part of Fig. 4.37 as a function of U. Since the 
magnetic moment per atom is given by 2\Apl- we see that a parameter U /|t,| = 4.8 
reproduces the experimental value of magnetic moment 0.5,p/Fe, therefore we can 
estimate the magnitude of the Coulomb interaction to be U = 4.6eV by using ft) = 
—0.07 Ry. 

From the discussion described in this section, we see that the double-helical mag- 
netic ordering of FeAs is well explained by the itinerant picture based on the multi- 
band Hubbard model whose parameters (transfer integrals and electron number) are 
obtained so as to reproduce results of the APW band calculations. This theory is 
widely applicable not only to other MnP-type compounds which has the double- 
helical ordering but also to compounds with several magnetic atoms in a unit cell. 


4.7 Phase Transition from the NiAs-Type 
to the MnP-Type Structure 


The structural phase transition between the NiAs-type and the MnP-type struc- 
tures, which is easily induced by temperature change or by composition change of 
compound, is one of the characteristic properties of the transition metal pnictides 
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and chalcogenides with the NiAs-type structure. Correlation between magnetic prop- 
erties and the crystal structure is an interesting problem which has not been fully 
understood yet. In transition metal chalcogenides, VS takes the MnP-type and the 
NiAs-type structures below and above 7; = 850K, respectively. Tig.95Se takes the 
NiAs-type structure, whereas Ti; 95Se takes the MnP-type structure at room temper- 
ature. TiyV;_,S takes the NiAs-type for x > 0.66 and the MnP-type for x < 0.66. As 
for the pnictides, MnAs, CrAs, and CoAs show the structural phase transition be- 
tween the NiAs-type structure above 7; and the MnP-type crystal structure below T, 
with 7, = 398K, 1,100K, 1,250K, respectively. Only MnAs, with decreasing tem- 
perature from 7; = 398 K, shows the second transformation from the MnP-type to 
the NiAs-type structures at 7; = 318 K, where magnetic phase transition to the fer- 
romagnetic order also occurs discontinuously (see also Fig. 3.17, etc.). In contrast to 
MnAs, etc., NiAs does not show the transformation to the MnP-type structure, how- 
ever the de Haas—van Alfen effect and the Fermi surface suggest an existence of a 
structural transformation at very low temperatures as noted in Sect. 4.4. 

Various measurements carried out for various mixed compounds revealed that the 
structural transformation between the NiAs-type and the MnP-type structures occurs 
in the following compounds: Cr;_,Mn,As (all x), Crj_;Co,As (all x) , Fe;_,,Co,As 
(0.8 <x < 1.0), Mnj_,Ni,As (0 < x < 0.58), Mnj_,Co,As (O< x <0.2,0.5<x< 
1.0), MnAs;_,P; (0 < x < 0.28) See also Sect. 3.2. In contrast to these arsenides, all 
the antimonides investigated so far do not show the transformation to the MnP-type 
structure. 

Taking account of the fact that the phase transition occurs in MnAs, CrAs, CoAs, 
and VS but not in MnSb, CrSb, NiAs, and TiS, we investigated mechanism of the 
structural phase transition by studying instability of the NiAs-type structure. This is 
a microscopic theory based on band structures obtained by the first principle calcu- 
lations, that is, no fitting parameters are used. 


4.7.1 Electron—Lattice Interaction Coefficient 


Change of energy of an electronic system AF due to lattice deformation is described 
by a phonon normal mode Q,,, (q is a wavevector, A denotes mode) is given by 


1 2 
AF =—5x (qa) |Q,a| - (4.36) 
7(q/) in this equation is a generalized susceptibility, which is given by 


ae %(ga,v)eP"(qd,v') x8 (vv',g) (4.37) 
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In these equations, My is mass of the v atom, €(q/,V) is the polarization vector of 
phonon, Ee is the band energy in the absence of lattice distortion, and f (E%) is the 
Fermi—Dirac distribution function. 

Displacement of a v atom in the o@-direction brings about change of crystalline 
potential AV, resulting in coupling between eigen states |nk) and |n’k +q) of the 
system. The coefficient Lilt dy in (4.38) denotes magnitude of the coupling between 
states ‘V,, and ‘Y’,,49 due to the potential change AV, and is called an electron—lattice 
interaction coefficient, which is given by 


Tenlk-tg = (Pn |AV Pweg) (4.39) 
An expression of the electron—lattice interaction coefficient based on the APW for- 
malism is given in [45]. 
4.7.2 Tendency of Structural Transformation from the NiAs-Type 


to the MnP-Type 


The structural transformation to the MnP-type is caused by freezing of a phonon 
mode M, at the M point of the first Brillouin zone for the NiAs-type structure. The 
normal coordinate of the M, mode is expressed as 


1 (Z1 +22) + ¢2(x, — x2) +.€3(z3 — 4), (4.40) 


where the subscripts 1 and 2 denote metal atoms, and 3 and 4 denote anion atoms. 
The symbol z is the displacement along the c-axis, and x (y) the displacement in the 
c-plane (see Fig. 4.38). As indicated in (4.40), metal atoms shift in the c-plane as well 
as along the c-axis, while anions only displace along the c-axis (see also Chap. 1). 


O metal @ anion 


Fig. 4.38. (a) The NiAs-type structure. The thick lines show the unit cell corresponding to 
the MnP-type structure. (b) Unit cell of the MnP-type structure. The arrows indicate displace- 
ments of atoms which vanish for the NiAs-type structure. We note that magnitude of the 
displacement is exaggerated. The lattice parameters a, b, and c in this figure correspond to A, 
B, C in Fig. 1.4, respectively 
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We calculated matrix elements of the electron—lattice interaction hae ig for g = 


ITM as a function of k. We note that when a wavevector k is on the Y-line (IM), 
k + q is also on the 2-line. There are four irreducible representations 2), 22, 3, and 
24 in the group of symmetry on the X-line. Among product representations, only 
x x X4 and Ly x 3 satisfy compatibility relation with the Mj representation. 
Therefore, matrix elements /’s vanish except for the combination of states 2j- 
and X4-bands or 2- and 2%3-bands. Furthermore, it is clear from (4.38) that the 
electron—lattice interaction reduces the free energy only when one of the states |nk) 
and |n'k’) is below the Fermi level and the other is above the Fermi level, and the 


reduction of free energy increases as the energy difference E . —E°_, decreases. 
ni n'k 


MnAs, MnSb, CrAs, and CrSb 


Dispersion curves along the X-line near the Fermi level calculated by using the APW 
method for MnAs, MnSb, CrAs, and CrSb are shown in Fig. 4.39. We see that the 
2}-band and 24-band lie near the Fermi level, so we calculated J rhe eta for the states 
k of the 2|-band and k +q of the X4-band. Calculated electron—lattice interaction co- 
efficients are plotted as a function of k along the X-line in Fig. 4.40. The symbols e, 
x, and o denote the results for the displacements in the x-direction of a metal atom, 
in the z-direction of a metal atom, and in z-direction of an anion, respectively. In 
the hatched regions, on the k-axis, the electron—lattice interaction contributes to the 
reduction of energy of electron system. Figure 4.41 shows the reduction of free en- 
ergy 2% sykrgl / Exe — Ey,x+9|. When this quantity is large in the hatched region 
on the k-axis shown in Fig. 4.41, large reduction of energy of the electron system 
due to the lattice distortions arises, therefore we can say that phase transition to the 
MnP-type structure is likely to occur. The results shown in Fig. 4.41 suggests that the 
phase transition to the MnP-type structure is expected to occur in MnAs and CrAs, 
but not in MnSb and CrSb. 


E (Ry) CrAs E (Ry) MnAs E (Ry) CrSb E(Ry) MnSb 
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Fig. 4.39. Dispersion curves of the nonmagnetic bands of MnAs, MnSb, CrAs, CrSb along 
the X-line. The broken line denotes the Fermi level. The numbers on the curves 1, 2, 3, and 4 
denote the irreducible representations 21, 22,23, and 24 of the bands (Reproduced form [45] 
with permission) 
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Fig. 4.40. Elecron-lattice matrix elements, |/¥% |, for MnAs, MnSb, CrAs, CrSb as func- 
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tions of the wavevector k. The wavevector q is fixed to I'M. The metal x-, metal z-, and anion 
z-components of the electron—lattice matrix elements represented by filled circles, times, and 
open circles, respectively. The hatched regions on the k-axis denote, where the coupling of 
electronic states due to lattice distortion reduces electronic energy (Reproduced form [45] 
with permission) 
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Fig. 4.41. Energy reduction |/ Se. Syke sl /\Ex\xk — Ex,k+q| (q =T-M) calculated as functions of 


k. The symbols filled circles, times, and open circles denote the metal x-, metal z-, and anion 
z-components of the MnP-type distortion (Reproduced with permission from [45]) 


We also see from the figure that the energy reduction due to the displacement of 
metal atoms in the x-direction is larger than that due to other two displacements. Ob- 
servation for MnAs in the MnP-type structure has revealed that the displacement of 
metal atoms in the x-direction is the largest in the three displacements. 
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From the discussion mentioned earlier, we see that wavevector dependence of the 
electron—lattice interaction coefficients and dispersion curves near the Fermi level 
have a decisive role in the phase transition form the NiAs-type to the MnP-type 
structure. 


CoAs and NiAs 


Using the APW method, we calculated the band structure of CoAs and NiAs near 
the Fermi level for k-points in the planes k, = 0, k; = 1/3IA, and k, = 2/3IA [46]. 
As an example of calculations, we show dispersion curves along the lines parallel 
to the X-line (IM) on the plane of k, =0 in Fig. 4.42. As we have noted earlier, 
the combination of 2;-band and 24-band (or the 22-band and 23-band) contribute 
to the reduction of electronic energy due to the MnP-type lattice distortion. Only 
when one of the states is above the Fermi level and the other is below the Fermi 
level, the electron—lattice interaction reduces electronic energy. The hatched region 
in Fig. 4.42 indicates where the reduction of electronic energy occurs. 


CoAs 


Fig. 4.42. Dispersion curves of the band of CoAs and NiAs near the Fermi level along the di- 
rection parallel to the I'M. the hatched regions on the k-axis mean where the energy reduction 
due to the electron—lattice interaction occurs (Reproduced with permission from [46]) 
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By calculating 25 ukeg! / Exe = Exjk+al> we found that reduction of elec- 
tronic energy due to the MnP-type lattice distortion is large for CoAs but small for 
NiAs, which indicates that structural transformation to the MnP-type structure is 
likely to occur in CoAs but not in NiAs. This is consistent with observation. Calcu- 
lated electron—lattice interaction coefficients as well as |JX7 5444 ?/|Exin — Exyetgl 
are depicted in [46]. 


VS and TiSe 


We carried out similar calculations for VS and TiSe. Dispersion curves near the 
Fermi level of these materials show that the electron—lattice interaction between the 
states k of the X3-band and k +q of the X-band is important. Calculated electron— 
lattice interaction coefficients show that reduction of free energy due to MnP-type 
lattice distortion is large for VS in comparison with that for TiSe, which suggests that 
VS has a possibility to deform to the MnP-type structure, however TiSe does not. 
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5 


Itinerant Electron Theory of Magnetism 
of Cu,Sb-Type Compounds 


5.1 Crystal Structure and Magnetic Properties 


The Cu2Sb-type compounds that consist of Cr, Mn, Fe, and As or Sb take tetragonal 
crystal structure. Figure 5.1 shows the crystal structure and the first Brillouin zone 
of the CuSb-type structure. The unit cell contains four magnetic atoms denoted by 
(1, 2, 3, 4) and two As or Sb atom denoted by (5, 6). Magnetic atoms occupy the 
sites M(I) and M(II), which have different symmetry. The Cu2Sb-type compounds 
have attracted much interest because they show various magnetic ordering such as 
ferromagnetic, ferrimagnetic, and antiferromagnetic orderings [1]. 

Magnetic structure with many magnetic atoms in a unit cell is described by two 
parameters Q and @;;. The vector Q (propagation vector) is a wavevector that spec- 
ifies angle between total moment in a unit cell and that in an adjacent unit cell, and 
;; is a relative angle between moments at the ith and the jth sites in a unit cell. The 
structural parameters and magnetic structure (Q, @;;, and magnetic moments at the 
respective sites) are shown in Tables 5.1 and 5.2, respectively. 

In CrMnAs and FeMnAs, which are mixed compounds between Mn2As and 
Cr2As or Fe2As, a Cr or Fe atom occupies mainly the M(I)-site and a Mn atom occu- 
pies mainly the M(ID-site. In MnAlGe and MnGaGe, which also take the Cu2Sb-type 
structure, Mn atoms occupy only the M(1)-site, and Al and Ga occupy the M(II)- 
site. The compounds MnAIGe and MnGaGe are ferromagnetic below T¢ = 503 and 
453 K, respectively. As the Mn atomic layers, as seen in Fig.5.1, are separated by 
the layers of Al (or Ga) and Ge, these compounds are of interest as a quasi two- 
dimensional ferromagnet. 


5.2 Band Structures of Cu2Sb-Type Compounds and Magnetic 
and Optical Properties 


Since early 1990s, Motizuki and her co-workers have carried out band calculations of 
the Cu2Sb-type compounds for magnetic and nonmagnetic phases to investigate var- 
ious magnetic properties as shown in Table 5.2 [2]. The self-consistent APW method 
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a 2 


Fig. 5.1. (a) Crystal structure and (b) the first Brillouin zone of the Cu2Sb-type structure 


Table 5.1. Lattice parameters of Cu2Sb-type compounds 


Lattice parameters 

a(A) c(A) ue Ve 
Cr2As 3.60 6.34 0.325 0.275 
CrMnAs 3.88 6.28 0.327 0.266 
Mn2As 3.78 6.25 0.33 0.265 
MnFeAs 3.735 6.035 (0.33 0.25 
Fe2As 3.627 5.973 0.33 0.265 
Mn2Sb 4.078 6.557 0.295 0.280 
Cu2Sb 3.99 6.09 0.27 0.30 
MnAlGe 3.914 5.933 0.273 0.280 
MnGaGe 3.963 5.895 0.29 0.29 


Data are taken from [1] 


and the LAPW method are employed [3] together with the muffin-tin approximation 
for crystal potential and the local density approximation for the exchange and corre- 
lation interaction. The formula of Gunnarsson and Lundqvist (GL) was also used for 
the LDA term [4]. As for relativistic effects, we neglected the spin-orbit interaction 
but scalar relativistic effects are taken into account. Densities of states calculated by 
the APW and LAPW methods for the nonmagnetic bands of the Cu2Sb-type com- 
pounds were found to agree well with each other. APW band calculations including 
the spin-orbit interaction were carried out by Kulatov and his co-workers for the 
Cu2Sb-type compounds. They also calculated optical conductivity from the band 
structure [5]. 
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Table 5.2. Magnetic structure of Cu2Sb-type compounds 


Magnetic 
moments 
Magnetic structure (Up/atom) 
Q diz $23 34. = M() = Md) 
Cro As AF (0,0,2/c) @ 1 0 0.40 1.34 
CrMnAs_ AF __ (0,0,0) T 1 1 0.41 3.14 
Mn, As AF (0,0,z/c) 0 1 1 2.2 41 
MnFeAs AF (0,0,a/c) 0 0 1 0.2 3.6 
Fe2As AF (0,0,z/c) 0 0 1 1.28 2.05 
Mn2Sb FI (0,0,0) 0 1 0 2.13 3.87 
Cu2Sb P - 
MnAlGe F (0,0,0) 0 0 0 1.70 - 
MnGaGe F (0,0,0) 0 0 0 166 —- 
Data are taken from [1] 
Cr.As Non-magnetic state 
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Fig. 5.2. Dispersion curves of nonmagnetic band of Cr2As 


5.2.1 Nonmagnetic State of Cr2As, Mn2As, Fe.As, Mn,Sb, CrMnAs, 
and FeMnAs 


We calculated dispersion curves and densities of states for these compounds in non- 
magnetic state [6]. As an example, we show the dispersion curves of Cr2As in 
Fig. 5.2. Two bands in the low energy region consist of As-4s orbitals. Above the gap, 
there are mixing bands between Cr-3d orbitals and As-4p orbitals. Total width of the 
p-d mixing bands is about 3-4 eV, which is large enough to regard the d electrons 
as itinerant electrons. Figure 5.3 shows the total and the partial densities of states 
of M2As (M = Cr, Mn, and Fe) and Mn2Sb. When we divide the p—d mixing band 
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Fig. 5.3. Densities of states of nonmagnetic (a) Cr2As, (b) Mn2As, (c) Fe2As and (d) Mn2Sb 


shown in Fig. 5.3 into three parts specified by (1), (2), and (3) from low energy side, 
it was revealed by our calculations that (1) and (3) have bonding and antibonding 
characters between M and As (or Sb), respectively, and (2) is regarded as nonbond- 
ing bands consisting almost of d obitals. In the nonbonding bands, d-orbitals at the 
M(1)-sites and M(ID-sites show strong hybridization as seen in Fig. 5.3. 

It is generally said from Fig. 5.3 that the d-band arising from M(II)-sites is nar- 
rower than that from M(I)-sites, which indicates that d electrons at the M(ID)-sites 
have weak itinerancy compared with d electrons at the M(I)-sites. Therefore, we ex- 
pect that magnitude of magnetic moment at the M(ID-site is larger than that at the 
M(1)-site. Both calculations and experiments actually indicate that magnetic moment 
of M(ID is larger than that of M(I) as shown in Table 5.3. 
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Table 5.3. Magnetic moments of Cu2Sb-type compounds evaluated by calculations 


Magnetic moments (Up) 


Total Md) Md) Anion 
Cro As = 0.34 (040) 137 (1.34) 0.03 
MnAs 172 (23) 3.50 (4.1) 0.04 
Fe) As = 1.01 (1.28) 211 (2.05) 0.03 
MnSb 1.76 (1.74) —2.11. (—2.13) 3.67. ~— (3.87) «0.01 
MnAlGe 1.81 (1.70) 1.90 —0.02 —0.06 
MnGaGe 2.22 (1.66) 2.31 —0.02 O07 


Observed moments are shown in the parenthesis 
Ey 
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Fig. 5.4. Density of states of nonmagnetic Cu2Sb 


Fermi level lies in the nonbonding band. As seen in Fig. 5.3, Mn2Sb has larger 
p(Ef), the density of states at Ep, while compounds other than Mn2Sb have smaller 
p(Ef) than that of Mn2Sb, which is consistent with the fact that Mn2Sb is ferrimag- 
netic and others are antiferromagnetic. 

Next, we show in Fig. 5.4 the density of states of Cu2Sb which shows no magnetic 
ordering [7]. We see that hybridization between d-orbitals of M() and M(ID) is weak, 
and density of states at Fermi level, p (Ep), is very small. 
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Fig. 5.5. Density of states of ferrimagnetic Mn Sb 


5.2.2 Ferrimagnetic Band of Mn2Sb 


We show calculated density of states of ferrimagnetic Mn2Sb in Fig. 5.5 [8]. Com- 
ponents of 3d-orbitals of Mn(I) and Mn(ID) and 5p-orbitals of Sb inside muffin-tin 
spheres are also plotted. It is noted that the shape of the densities of states for the 
up-spin and down-spin bands are quite different from each other in the respective 
components, and also different from the densities of states in the nonmagnetic states. 
The band component of Mn(II) shows a large splitting between up- and down-spin 
bands, while the splitting of the band of Mn(I) is small, which results in that the mag- 
netic moment of Mn(II) is about two times larger than that of Mn(I) (see Table 5.3). 
We also see that small magnetic moment is induced at the Sb-site in the direction of 
moment of Mn(II). Calculated total magnetic moment as well as those of Mn(I) and 
Mn(ID) agrees well with experimental vales as seen in Table 5.3. 

The ferrimagnetic band of Mn2Sb is also calculated by Haas and his co-workers 
using the ASW method [9, 10]. Kulatov also carried out band calculation using the 
LMTO method. Gross features of density of states obtained by them are similar to 
the result shown in Fig. 5.5. 
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5.2.3 Ferromagnetic State of MnAlGe and MnGaGe 


In these compounds, atomic planes of Mn at M(I)-site are separated by two planes 
consisting of Al (or Ga) and Ge, which suggests a two-dimensional character in mag- 
netic properties. Observed magnetic moment is much smaller than moment expected 
for localized Mn orbitals. Therefore, we can regard to some extent these compounds 
as an itinerant electron ferromagnets with two-dimensional character. 

Figure 5.6 shows calculated density of states for MnAlGe [11]. MnAlGe has 
density of states similar to that of MnGaGe. Hybridization between d-orbitals of Mn 
and p-orbitals of Al (Ga) and Ge is small. As shown in Table 5.3, calculated moment 
of MnAIJGe is in good agreement with the experimental value, while the agreement 
is not good for MnGaGe. According to the band calculations, very small magnetic 
moment is induced at Al (Ga) and Ge atomic site in the opposite direction of the 
Mn-moment. 

We found, as shown in Fig.5.7, three Fermi surfaces of the down-spin bands 
consisting of two electron surfaces and a cylindrical hole surface around the Z-axis, 
and four Fermi surfaces of the up-spin bands consisting of a small hole pocket around 
the R-point, a big hole pocket, and two electron pockets near the M-point. The Fermi 
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Fig. 5.6. Density of states of ferromagnetic bands of MnAlGe 
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Fig. 5.7. Fermi surfaces of ferromagnetic bands of MnAlGe 


surfaces of the down-spin bands has large contribution to p (Ef). We expect that the 
strong two-dimensional character of the Fermi surfaces of the down-spin bands will 
be observed as anisotropy of transport properties such as resistivity. 


5.2.4 Antiferromagnetic Bands of Cr2As, Mn2As, and Fe2As 


We carried out band calculations for antiferromagnetic state of these compounds. As 
seen in Fig. 2.2 of Part I and Fig. 5.1 of Part II, magnetic moments of Cr(1) and Cr(2) 
of CroAs are in the c-plane and antiparallel, and so Cr(1) and Cr(2) are regarded 
as different sites. On the other hand, in Mn2As and Fe2As, the sites 1 and 2 are 
equivalent because momentat these sites are parallel. Magnetic moments at the M(1)- 
site and adjacent M(ID-site are antiparallel in Mn2As and parallel in Fe2As. The spin 
arrangements for the compounds with Cu2Sb-type structure are illustrated in Fig. 2.2 
of Part I. 
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Fig. 5.8. Total and local densities of states for (a) Cr2As, (b) Mn2As, and (c) Fe2As in anti- 
ferromagnetic state 


Figure 5.8 shows total and local densities of states for MAs with M = Cr, Mn, and 
Fe. It is noteworthy that in Cr2 As the splitting of Cr(II)-bands is about | eV, whereas 
there is little splitting of Cr(I)-bands. The very small splitting in Cr(I)-local density 
of states results in the very small magnetic moment of Cr(1), which is related to the 
antiferromagnetic order in the c-plane. As for FezAs and Mn2As, mixing between 
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3d-states at the M(I)- and the M(II)-sites in Mn2As is smaller than that in Fe2As. 
Calculated moments are compared with the experimental data in Table 5.3. As shown 
in Table 5.3, magnetic moment of M(ID) is larger than that of M(D) in the respective 
compounds. Agreement between calculation and experiment is good for Cr2As and 
Fe,As but not so good for Mn2As. 


5.2.5 Magnetic Ordering of Cu,Sb-Type Compounds 


The Cu2Sb-type compounds except for MnAlGe and MnGaGe have four magnetic 
atoms in a unit cell, and therefore various types of magnetic structures have been 
observed. We have studied such various magnetic structures from the viewpoint of 
itinerant electrons. For this picture we investigated instability of paramagnetic phase 
within the theoretical framework described in Sect. 4.6. In this section we explain 
only the calculated results. 

Figure 5.9 shows the magnetic phase diagram calculated for Q = (0,0,0) and 
Q = (0,0,2/c), respectively. The dashed-curves express the relation between the 
critical values of intra-atomic Coulomb interaction divided by a transfer integral, 
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Fig. 5.9. A magnetic phase diagram of the Cu2Sb-type compounds. The symbols open circle, 
open triangle, open square, etc. denote critical values of intra-atomic Coulomb interaction 
calculated for each magnetic ordering (see also the text). The magnetic orders expressed by 
CrMnAs(II)-type and Cr2As-type are hypothetical ones for the theoretical calculation [14] 
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Ue /|ti|. Uc/|ti| is plotted as a function of average electron number n. As we have 
noted in Sect. 4.6, if U is larger than U, paramagnetic phase is unstable. The magnetic 
structure that corresponds to the smallest value of U; is most likely to realize. 

From the calculations, it is indicated that Mn2Sb-type, CrMnAs(II)-type, and fer- 
romagnetic ordering are theoretically possible for Q = (0,0,0). For Q = (0,0,2/c), 
the Cr2 As(ID)-type, Mn2As-type, and Fe2As-type can theoretically realize. The mag- 
netic orderings actually observed are included in Fig. 5.9. 


5.2.6 Photoemission and Inverse Photoemission 


It is possible to obtain a rough estimation of density of states below and above the 
Fermi level by photoemission and inverse photoemission measurements, respec- 
tively. Suga and Kimura measured photoemission and inverse photoemission spec- 
tra of Cu2As, MnoAs, FezAs, MnAlGe, and MnGaGe. Their experimental data are 
shown in Fig. 5.10 [12, 13] together with the calculated density of states summed 
over the up-spin and down-spin bands [14]. Overall features of the spectra corre- 
spond well to the calculations. For Mn2Sb, the spectrum has remarkable peaks at 
about 3 eV below the Fermi level and 2eV above the Fermi level, while in MnAlGe 
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Fig. 5.10. Densities of states determined by photoemission and inverse photoemission spectra 
(dots) and calculated densities of states (solid lines) for the Cu2Sb-type compounds 
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there is no clear peak but broad hump with the width about 4.5 eV lies below Fermi 
level. Such difference between the spectra of Mn2Sb and MnAI]Ge can be attributed 
to the absence of Mn at M(ID) site in the case of MnAIGe. 

From the spectra of Cr2As, Mn2As, and Fe2As, we see that energy difference 
between the peaks above and below the Fermi level in Cr2As and Fe2As is smaller 
than the energy difference in Mn As. 
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